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Abstract. We consider the 3-dimensional gravitational n-body 
problem, n > 2, in spaces of constant Gaussian curvature k ^ 0, 
i.e. on spheres for n > 0, and on hyperbolic manifolds H^, 
for k < 0. Our goal is to define and study relative equilibria, 
which are orbits whose mutual distances remain constant in time. 
We also briefly discuss the issue of singularities in order to avoid 
impossible configurations. We derive the equations of motion and 
define six classes of relative equilibria, which follow naturally from 
the geometric properties of and Hj. Then we prove several 
criteria, each expressing the conditions for the existence of a certain 
class of relative equilibria, some of which have a simple rotation, 
whereas others perform a double rotation, and we describe their 
qualitative behaviour. In particular, we show that in the bodies 
move either on circles or on Clifford tori, whereas in they move 
either on circles or on hyperbolic cylinders. Then we construct 
concrete examples for each class of relative equilibria previously 
described, thus proving that these classes are not empty. We put 
into the evidence some surprising orbits, such as those for which 
a group of bodies stays fixed on a great circle of a great sphere of 
S^, while the other bodies rotate uniformly on a complementary 
great circle of another great sphere, as well as a large class of 
quasiperiodic relative equilibria, the first such non-periodic orbits 
ever found in a 3-dimcnsional n-body problem. Finally, we briefly 
discuss other research directions and the future perspectives in the 
light of the results we present here. 
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Part 1. INTRODUCTION 

In this introductory part, we will define the problem, explain its 
importance, tell its history, and outline the main results we obtained. 

1. The problem 

We consider in this paper a natural extension of the Newtonian n- 
body problem to spaces of non-zero constant curvature. Since there is 
no unique way of generalizing the classical equations of motion such 
that they can be recovered when the space in which the bodies move 
flattens out, we seek a potential that satisfies the same basic properties 
as the Newtonian potential in its simplest possible setting, that of one 
body moving around a fixed centre, which is a basic problem in celestial 
mechanics known as the Kepler problem. Two basic properties stick 
out in this case: the Newtonian potential is a harmonic function in 3- 
dimensional space, i.e. it satisfies Laplace's equation, and it generates 
a central field in which all bounded orbits are closed, a result proved 
by Joseph Louis Bertrand in 1873, [5]. 

On one hand, the curved potential we define in Section [6] approaches 
the classical Newtonian potential when the curvature tends to zero, 
whether through positive or negative values. In the case of the Kepler 
problem, on the other hand, this potential satisfies Bertrand's prop- 
erty and is a solution of the Laplace-Beltrami equation, [39], the nat- 
ural generalization of Laplace's equation to Riemannian and pseudo- 
Riemannian manifolds, which include the spaces of constant curvature 
k ^ we are interested in, namely the spheres S^, for k > 0, and the 
hyperbolic manifolds H^, for k < 0. For simplicity, from now on we 
will refer to the n-body problem defined in these spaces of constant 
curvature as the curved n-body problem. 

In the Euclidean case, the Kepler problem and the 2-body problem 
are equivalent. The reason for this overlap is the existence of the in- 
tegrals of the centre of mass and of the linear momentum. It can be 
shown with their help that the equations of motion appear identical, 
whether the origin of the coordinate system is fixed at the centre of 
mass or fixed at one of the two bodies. For non-zero curvature, how- 
ever, things change. As we will later see, the equations of motion of the 
curved n-body problem lack the integrals of the centre of mass and of 
the linear momentum, which prove to characterize only the Euclidean 
case. Consequently the curved Kepler problem and the curved 2-body 
problem are not equivalent anymore. It turns out that, as in the Eu- 
clidean case, the curved Kepler problem is integrable, but, unlike in the 
Euclidean case, the curved 2-body problem is not, [TS]. As expected, 
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the curved n-body problem is not integrable for n > 3, a property that 
is well-known to be true in the Euclidean case. 

Our main goal is to find solutions that are specific to each of the 
spaces corresponding to n < 0, k = 0, and n > 0. We already succeeded 
to do that in some previous papers, but only in the 2-dimensional 
case, [22], [25], [26], [27], [28]. In this paper, we will see that the 3- 
dimensional curved n-body problem puts into the evidence even more 
differences between the qualitative behaviour of orbits in each of these 
spaces. 

2. Importance 

The curved n-body problem lays a bridge between the theory of 
dynamical systems and the geometry and topology of 3-dimensional 
manifolds of constant curvature, as well as with the theory of regu- 
lar polytopes. As the results obtained in this paper show, many of 
which are surprising and nonintuitive, the geometry of the configura- 
tion space (i.e. the space in which the bodies move) strongly influences 
their dynamics, including some important qualitative properties, stabil- 
ity among them. Some topological concepts, such as the Hopf fibration 
and the Hopf link, some geometric objects, such as the pentatope, the 
Clifford torus, and the hyperbolic cylinder, or some geometric proper- 
ties, such as the Hegaard splitting of genus 1, become essential tools 
for understanding the gravitational motion of the bodies. Since we 
provide here only a first study in this direction of research, reduced to 
the simplest orbits the equations of motion have, namely the relative 
equilibria, we expect that many more connections between dynamics, 
geometry, and topology will be discovered in the near future through 
a deeper exploration of the 3-dimensional case. 

3. History 

The first researchers that took the idea of gravitation beyond the Eu- 
clidean space were Nikolai Lobachevsky and Janos Bolyai, who also laid 
the foundations of hyperbolic geometry independently of each other. In 
1835, Lobachevsky proposed a Kepler problem in the 3-dimensional hy- 
perbolic space, H 3 , by defining an attractive force proportional to the 
inverse area of the 2-dimensional sphere of radius equal to the distance 
between bodies, [58]. Independently of him, and at about the same 
time, Bolyai came up with a similar idea, [BJ. Both of them under- 
stood the intimate connection between geometry and physical laws, a 
relationship that proved very prolific ever since. 
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These co- discoverers of the first non-Euclidean geometry had no fol- 
lowers in their pre-relativistic attempts until 1860, when Paul Joseph 
SerrelQ extended the gravitational force to the sphere § 2 and solved 
the corresponding Kepler problem, [71]. Ten years later, Ernst Scher- 
ing revisited Lobachevsky's gravitational law for which he obtained an 
analytic expression given by the curved cotangent potential we study 
in this paper, [7TJ. Schering also wrote that Lejeune Dirichlet had told 
some friends to have dealt with the same problem during his last years 
in Berlin^ [72], but Dirichlet never published anything in this direc- 
tion, and we found no evidence of any manuscripts he would have left 
behind in which he dealt with this particular topic. In 1873, Rudolph 
Lipschitz considered the problem in S 3 , but defined a potential propor- 
tional to 1/ sin J^, where r denotes the distance between bodies and R 
is the curvature radius, [57]. He obtained the general solution of this 
problem only in terms of elliptic functions. His failure to provide an 
explicit formula, which could have him helped draw some conclusions 
about the motion of the bodies, showed the advantage of Schering's 
approach. 

In 1885, Wilhelm Killing adapted Lobachevsky's gravitational law 
to § 3 and defined an extension of the Newtonian force given by the 
inverse area of a 2-dimensional sphere (in the spirit of Schering), for 
which he proved a generalization of Kepler's three laws, [UJ. Then 
an important step forward took place. In 1902, Heinrich Liebmann^] 
showed that the orbits of the Kepler problem are conies in § 3 and H 3 
and further generalized Kepler's three laws to k ^ 0, [53]. One year 
later, Liebmann proved S> 2 - and H 2 -analogues of Bertrand's theorem, 
[5], [82], which states that for the Kepler problem there exist only two 
analytic central potentials in the Euclidean space for which all bounded 
orbits are closed, [M]- He thus made the most important contributions 
of the early times to the model, by showing that the approach of Bolyai 



Paul Joseph Serret (1827-1898) should not be confused with another French 
mathematician, Joseph Alfred Serret (1819-1885), known for the Frenet-Serret for- 
mulas of vector calculus. 

2 This must have happened around 1852, as claimed by Rudolph Lipschitz, |56j . 

3 Although he signed his papers and books as Heinrich Liebmann, his full name 
was Karl Otto Heinrich Liebmann (1874-1939). He did most of his work in Munich 
and Heidelberg, where he became briefly the university's president, before the Nazis 
forced him to retire. A remembrance colloquium was held in his honour in Hei- 
delberg in 2008. On this occasion, Liebmann's son, Karl-Otto Liebman, donated 
to the University of Heidelberg an oil portrait of his father, painted by Adelheid 
Liebmann, [75] 
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and Lobachevsky led to a natural extension of Newton's gravitational 
law to spaces of constant curvature. 

Liebmann also summed up his results in the last chapter of a book 
on hyperbolic geometry published in 1905, [55], which saw two more 
editions, one in 1912 and the other in 1923. Intriguing enough, in the 
third edition of his book he replaced the constant-curvature approach 
with relativistic considerations. 

Liebmann's change of mind about the importance of the constant- 
curvature approach may explain why this direction of research was 
ignored in the decades immediately following the birth of special and 
general relativity. The reason for this neglect was probably connected 
to the idea that general relativity could allow the study of 2-body prob- 
lems on manifolds with variable Gaussian curvature, so the constant- 
curvature case appeared to be outdated. Indeed, although the most 
important subsequent success of relativity was in cosmology and re- 
lated fields, there were attempts to discretize Einstein's equations and 
define a gravitational n-body problem. Remarkable in this direction 
were the contributions of Jean Chazy, [TT], Tullio Levi-Civita, |50j . 
[5T] . Arthur Eddington, [33], Albert Einstein, Leopold InfelcQ Banesh 
Hoffmann, [31], and Vladimir Fock, [3E]. Subsequent efforts led to re- 
fined post-Newtonian approximations (see, e.g., [16], [17], [IE]), which 
prove very useful in practice, from understanding the motion of artifi- 
cial satellites — a field with applications in geodesy and geophysics — to 
using the Global Positioning System (GPS), [T9"] . 

But the equations of the n-body problem derived from relativity are 
highly complicated even for n = 2, and they are not prone to analytical 
studies similar to the ones done in the classical case. This is probably 
the reason why the need for some simpler equations revived the research 
on the motion of 2 bodies in spaces of constant curvature. 

Starting with 1940, Erwin Schrodinger developed a quantum-mecha- 
nical analogue of the Kepler problem in § 2 , [73]. Schrodinger used the 
same cotangent potential of Schering and Liebmann, which he deemed 
to be the natural extension of Newton's law to the sphere^} Further 
results in this direction were obtained by Leopold Infeld, [43], [SI]- In 
1945, Infeld and his student Alfred S child extended this problem to 
spaces of constant negative curvature using a potential given by the 
hyperbolic cotangent of the distance, [15] . A comprehensive list of the 



A vivid description of the collaboration between Einstein and Infeld appears in 
Infcld's autobiographical book [44] . 

5 "The correct form of [the] potential (corresponding to 1/r of the fiat space) is 
known to be cotx," [73], p. 14. 
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above-mentioned works also appears in [76], except for Serret's book, 
|74j . which is not mentioned. A bibliography of works about mechanical 
problems in spaces of constant curvature is given in [79], [78] . 

The Russian school of celestial mechanics led by Valeri Kozlov also 
studied the curved 2-body problem given by the cotangent potential 
and considered related problems in spaces of constant curvature start- 
ing with the 1990s, [49J. An important contribution to the case n = 2 
and the Kepler problem belongs to Jose Carinena, Manuel Ranada, and 
Mariano Santander, who provided a unified approach in the framework 
of differential geometry with the help of intrinsic coordinates, empha- 
sizing the dynamics of the cotangent potential in S> 2 and H 2 , [9] (see 
also [10], [UJ). They also proved that, in this unified context, the conic 
orbits known in Euclidean space extend naturally to spaces of constant 
curvature, in agreement with the results obtained by Liebmann, |76j . 
Moreover, the authors used the rich geometry of the hyperbolic plane 
to search for new orbits, whose existence they either proved or conjec- 
tured. 

Inspired by the work of Carinena, Ranada, and Santander, we pro- 
posed a new setting for the problem, which allowed us an easy deriva- 
tion of the equations of motion for any n > 2 in terms of extrinsic 
coordinates, [26], [27]. The combination of two main ideas helped 
us achieve this goal: the use of Weierstrass's hyperboloid model of 
hyperbolic geometry and the application of the variational approach 
of constrained Lagrangian dynamics, |39| . Although we obtained the 
equations of motion for any dimension, we explored so far only the 
2-dimensional case. In [26] and [27], we studied relative equilibria and 
solved Saari's conjecture in the collinear case (see also [29], [30]), in 
[28] and [2T] we studied the singularities of the curved n-body prob- 
lem, in [23] we gave a complete classification of the homographic solu- 
tions in the 3-body case, and in [22] we obtained some results about 
polygonal homographic orbits, including a generalization of the Perko- 
Walter-Elmabsout theorem, [35], [68]. Knowing the correct form of the 
equations of motion, Ernesto Perez-Chavela and J. Guadalupe Reyes 
Victoria succeeded to derive them in intrinsic coordinates in the case of 
positive curvature and showed them to be equivalent with the extrin- 
sic equations obtained with the help of variational methods, [67J. We 
are currently developing an intrinsic approach for negative curvature, 
whose analysis appears to be more complicated than the study of the 
positive curvature case, [25] . 
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4. Results 

This paper is a first attempt to study the 3-dimensional curved n- 
body problem in the general context described in the previous section, 
with the help of the equations of motion written in extrinsic coordi- 
nates, as they were obtained in [26] and [27]. We are mainly concerned 
with understanding the motion of the simplest possible orbits, the rel- 
ative equilibria, which move like rigid bodies, by maintaining constant 
mutual distances for all time. 

We structured this paper in 6 parts. This is Part 1, which defines 
the problem, explains its importance, gives its historical background, 
describes the structure of the paper, and presents the main results. In 
Part 2 we set the geometric background and derive the equations of 
motion in the 3-dimensional case. In Part 3 we discuss the isometric 
rotations of the spaces of constant curvature in which the bodies move 
and, based on this development, define several natural classes of rela- 
tive equilibria. We also analyze the occurrence of fixed points. In Part 
4 we prove criteria for the existence of relative equilibria and a give 
qualitative description of how these orbits behave. In Part 5 we pro- 
vide concrete examples of relative equilibria for all the classes of orbits 
previously discussed. Part 6 concludes the paper with a short analysis 
of the recent achievements obtained in this direction of research and a 
presentation of the future perspectives our results offer. 

Each part contains several sections, which are numbered increasingly 
throughout the paper. Sections 1, 2, 3, and 4 form Part 1. In Section 5, 
which starts Part 2, we lay the background for the subsequent develop- 
ments that appear in this work. We thus introduce in 5.1 Weierstrass's 
model of hyperbolic geometry, outline in 5.2 its history, define in 5.3 
some basic concepts of geometric topology, describe in 5.4 the metric 
in curved space, and introduce in 5.5 some functions that unify cir- 
cular and hyperbolic trigonometry. In Section 6, we start with 6.1 in 
which we give the analytic expression of the potential, derive in 6.2 
Euler's formula, present in 6.3 the theory of constrained Lagrangian 
dynamics and then use it in 6.4 to obtain the equations of motion of 
the n-body problem on 3-dimensional manifolds of constant Gaussian 
curvature, k^O, i.e. spheres, for k > 0, and hyperbolic manifolds, 
H^, for k < 0. Then we prove in 6.5 that the equations of motion are 
Hamiltonian and show in 6.6, Proposition 1, that the 2-dimensional 
manifolds of the same curvature, and H^, are invariant for these 
equations. In Section 7 we prove the existence and derive the expres- 
sions of the 7 basic integrals this Hamiltonian system has, namely we 
come up in 7.1 with 1 integral of energy and in 7.2 with the 6 integrals 
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of the total angular momentum. Unlike in the Euclidean case, there 
are no integrals of the centre of mass and linear momentum. Although 
we don't aim to study the singularities of the equations of motion here, 
we describe them in Section 8, Proposition 2, in order to avoid singular 
configurations when dealing with relative equilibria. 

Section 9, which starts Part 3, describes the isometric rotations of 
E> 3 and H^, a task that helps us understand the various types of trans- 
formations and how to use them to define natural classes of relative 
equilibria. It turns out that we can have: (1) simple elliptic rotations 
in §^; (2) double elliptic rotations in E> 3 ; (3) simple elliptic rotations in 
M 3 K ; (4) simple hyperbolic rotations in M 3 ; (5) double elliptic-hyperbolic 
rotations in M. 3 ; and (6) simple parabolic rotations in M 3 . Section 10 
continues this task by exploring when 2-dimensional spheres, in 10.1, 
and hypeboloids, in 10.2, of curvature k or different from k are pre- 
served by the isometric rotations of §^ and M 3 K , respectively. Based on 
the results we obtain in these sections, we define 6 natural types of 
relative equilibria in Section 11: /t-positive elliptic, in 11.1; /t-positive 
elliptic-elliptic, in 11.2; /t-negative elliptic, in 11.3; ^-negative hyper- 
bolic, in 11.4; /t-negative elliptic-hyperbolic, in 11.5; and /t-negative 
parabolic, in 11.6; then we summarize them all in 11.7. Since relative 
equilibria can also be generated from fixed-point configurations, we 
study fixed-point solutions in Section 12, which closes Part 3. It turns 
out that fixed-point solutions occur in E> 3 , as we prove in 12.1, but not 
in M 3 (Proposition 3 of 12.3) or hemispheres of E> 3 (Proposition 4 of 
12.3). Moreover, some of these fixed points are specific to E> 3 , such as 
when 5 equal masses lie at the vertices of a regular pentatope, because 
there is no 2-dimensional sphere on which this configuration occurs, as 
we show in 12.2. 

Section 13 opens Part 4 of this paper with proving 7 criteria for the 
existence of relative equilibria as well as the nonexistence of /t-negative 
parabolic relative equilibria, the latter in 13.6 (Proposition 5). Criteria 
1 and 2 of 13.1 are for /t-positive elliptic relative equilibria, both in the 
general case and when they can be generated from fixed-point configu- 
rations, respectively. Criteria 3 and 4 of 13.2 repeat the same pattern, 
but for /t-positive elliptic-elliptic relative equilibria. Finally, Criteria 5, 
6, and 7, proved in 13.3, 13.4, and 13.5, respectively, are for ^-negative 
elliptic, hyperbolic, and elliptic-hyperbolic relative equilibria, respec- 
tively. In Section 14 we prove the main results of this paper after laying 
some background of geometric topology. In 14.1 we introduce the con- 
cept of Clifford torus. In 14.2 we prove Theorem 1, which describes 
the general qualitative behaviour of relative equilibria in E> 3 , assuming 
that they exist, and shows that the bodies could move on circles or on 
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Clifford tori. In 14.3 we prove Theorem 2, which becomes more specific 
about relative equilibria generated from fixed-point configurations; it 
shows that we could have orbits for which some points are fixed while 
others rotate on circles, as well as a large class of quasiperiodic orbits, 
which are the first kind of non-periodic relative equilibria ever hinted 
at in celestial mechanics. In Section 15, which closes Part 4, we study 
what happens in H^. In 15.1 we introduce the concept of hyperbolic 
cylinder and in 15.2 we prove Theorem 3, which describes the possible 
motion of relative equilibria in H^; it shows that the bodies could move 
only on circles, on geodesies, or on hyperbolic cylinders. 

Since none of the above results show that such relative equilibria 
actually exist, it is necessary to construct concrete examples with the 
help of Criteria 1 to 7. This is the goal of Part 5, which starts with 
Section 16, in which we construct 4 classes of relative equilibria with 
a simple rotation in S^. Example 16.1 provides a Lagrangian solution: 
the bodies, of equal masses, are at the vertices of an equilateral trian- 
gle) in the curved 3-body problem in S 3 K and rotate on a non-necessarily 
geodesic circle of a great sphere, the frequency of rotation depending 
on the masses and on the circle's position. Example 16.2 generalizes 
the previous construction in the particular case of the curved 3-body 
problem when the bodies rotate on a great circle of a great sphere. 
Then, for any acute scalene triangle, we can find masses that lie at 
its vertices while the triangle rotates with any non-zero frequency. In 
Example 16.3, the concept of complementary great circles, meaning 
great circles lying in the planes wx and yz, plays the essential role. 
We construct a relative equilibrium in the curved 6-body problem in 
for which 3 bodies of equal masses are fixed at the vertices of an 
equilateral triangle on a great circle of a great sphere, while the other 
3 bodies, of the same masses as the previous 3, rotate at the vertices of 
an equilateral triangle along a complementary great circle of another 
great sphere. We show that these orbits cannot be contained on any 
2-dimensional sphere. Example 16.4 generalizes the previous construc- 
tion when the triangles are acute and scalene and the masses are not 
necessarily equal. 

Section 17 provides examples of relative equilibria with double ellip- 
tic rotations in In Example 17.1 an equilateral triangle with equal 
masses at its vertices has 2 elliptic rotations generated from a fixed- 
point configuration in the curved 3-body problem. Example 17.2 is in 
the curved 4-body problem: a regular tetrahedron with equal masses 
at its vertices has 2 elliptic rotations of equal frequencies. In Example 
17.3 we construct a solution of the curved 5-body problem in which 5 
equal masses are at the vertices of a regular pentatope that has two 



Relative equilibria in the 3-dimensional curved n-body problem 



13 



elliptic rotations of equal-size frequencies. As in the previous example, 
the orbit is generated from a fixed point, the frequencies have equal 
size, and the motions cannot take place on any 2-dimensional sphere. 
Example 17.4 is in the curved 6-body problem, for which 3 bodies of 
equal masses are at the vertices of an equilateral triangle that rotates 
on a great circle of a great sphere, while the other 3 bodies, of the same 
masses as the others, are at the vertices of another equilateral triangle 
that moves along a complementary great circle of another great sphere. 
The frequencies are distinct, in general, so, except for a negligible set 
of solutions that are periodic, the orbits are quasiperiodic. This seems 
to be the first class of solutions in n-body problems given by various 
potentials for which relative equilibria are not periodic orbits. Exam- 
ple 17.5 generalizes the previous construction to acute scalene triangles 
and non-equal masses. Sections 18, 19, and 20 provide, each, an exam- 
ple of a ^-negative elliptic relative equilibrium, a /t-negative hyperbolic 
relative equilibrium, and a /t-negative elliptic-hyperbolic relative equi- 
librium, respectively, in H^. These 3 examples are in the curved 3-body 
problem. So Part 5 proves that all the criteria we proved in Part 4 can 
produce relative equilibria, which behave qualitatively exactly as The- 
orems 1, 2, and 3 predict. 

Finally, Part 6 concludes the paper with a short discussion about 
the current research of solutions of the curved n-body problem and 
their stability (Section 21), and about future perspectives for relative 
equilibria in general, in 22.1, and their stability, in 22.2. 
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Part 2. BACKGROUND AND EQUATIONS OF MOTION 

The goal of this part of the paper is to lay the mathematical back- 
ground for future developments and results and to obtain the equations 
of motion of the curved n-body problem together with their first inte- 
grals. We will also identify the singularities of the equations of motion 
in order to avoid impossible configurations for the relative equilibria 
we are going to construct in Part 5. 

5. Preliminary developments 

In this section we introduce some concepts that will be needed for the 
derivation of the equations of motion of the n-body problem in spaces 
of constant curvature as well as for the study of the relative equilibria 
we aim to investigate in this paper. 

The standard models of 2-dimensional hyperbolic geometry are the 
Poincare disk and the Poincare upper-half plane, which are conformal, 
i.e. maintain the angles existing in the hyperbolic plane, as well as 
the Klein-Beltrami disk, which is not conformal. But none of these 
models will be used in this paper. In the first part of the section, we 
introduce the less known Weierstrass model, which physicists usually 
call the Lorentz model, and we extend it to the 3-dimensional case. 
This model is more natural than the ones previously mentioned in the 
sense that it is dual to the sphere, and will thus be essential in our 
endeavours to develop a unified n-body problem in spaces of constant 
positive and negative curvature. We then provide a short history of 
Weierstrass's model, to justify its name, and introduce some geometry 
concepts that will be useful later. In the last part of this section, we 
define the metric that will be used in the model and unify circular and 
hyperbolic trigonometry, such that we can introduce a single potential 
for both the positive and the negative curvature case. 

5.1. The Weierstrass model of hyperbolic geometry. Since the 
Weierstrass model of hyperbolic geometry is not widely known among 
nonlinear analysts or experts in differential equations and dynamical 
systems, we present it briefly here. We first discuss the 2-dimensional 
case, which we will then extend to 3 dimensions. In its 2-dimensional 
form, this model appeals for at least two reasons: it allows an obvious 
comparison with the 2-dimensional sphere, both from the geometric 
and from the algebraic point of view; it emphasizes the difference be- 
tween the hyperbolic (Bolyai-Lobachevsky) plane and the Euclidean 
plane as clearly as the well-known difference between the Euclidean 
plane and the sphere. From the dynamical point of view, the equations 
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of motion of the curved n-body problem in S 3 . resemble the equations 
of motion in H 3 ., with just a few changes of sign, as we will see later. 
The dynamical consequences will be, nevertheless, significant, but we 
will still be able to use the resemblances between the sphere and the 
hyperboloid in order to understand the dynamics of the problem. 

The 2-dimensional Weierstrass model is built on one of the sheets of 
the hyperboloid of two sheets, 



2,2 2 -1 

x + y — z — k 



where k < 
Minkowski space 



represents the curvature of the surface in the 3-dimensional 
Dace IR 2 - 1 := (IR 3 , □), in which 



aHb 



with a = (a x ,a y ,a z ) and b = (b x ,b y ,b z 



defines the Lorentz inner 
product. We choose for our model the z > sheet of the hyperboloid 
of two sheets, which we identify with the hyperbolic plane H 2 . We can 
think of this surface as being a pseudosphere of imaginary radius iR, 
a case in which the relationship between radius and curvature is given 
by (iR) 2 = kT 1 . 

A linear transformation T: IR 2,1 — y IR 2 ' 1 is called orthogonal if 



T(a) HT(a) = a Ha 

for any a G IR 2,1 . The set of these transformations, together with the 
Lorentz inner product, forms the orthogonal group C^IR 2,1 ), given by 
matrices of determinant ±1. Therefore the group .SC^IR 2,1 ) of orthogo- 
nal transformations of determinant 1 is a subgroup of (^(IR 2 ' 1 ). Another 
subgroup of ©(IR 2,1 ) is (^(IR 2 ' 1 ), which is formed by the transformations 
T that leave H 2 invariant. Furthermore, (^(IR 2 ' 1 ) has the closed Lorentz 
subgroup, Lor^ 2 ' 1 ) := GiM 2 ' 1 ) n SOiM 2 ' 1 ). 

An important result, with consequences in our paper, is the Principal 
Axis Theorem for Lor(IR 2 ' 1 ), [31], [66J. Let us define the Lorentzian 
rotations about an axis as the 1-parameter subgroups of Lor(IR 2 ' 1 ) that 
leave the axis pointwise fixed. Then the Principal Axis Theorem states 
that every Lorentzian transformation has one of the forms: 



A = P 



cos 9 — sin 9 
sin 9 cos 9 
1 



B = P 



1 
cosh s 
sinh s 





sinh s 
cosh s 
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or 



C = P 



1 

t 



-t t 

1 - t 2 /2 t 2 /2 P 
-t 2 /2 1 + 1 2 /2 



-i 



where 9 G [0, 2ir), s,t el, and P e Lor(IR 2 ' 1 ). These transformations 
are called elliptic, hyperbolic, and parabolic, respectively. The elliptic 
transformations are rotations about a timelike axis — the z axis in our 
case — and act along a circle, like in the spherical case; the hyperbolic 
rotations are about a spacelike axis — the x axis in our context — and 
act along a hyperbola; and the parabolic transformations are rotations 
about a lightlike (or null) axis — represented here by the line x — 0, 
y = z — and act along a parabola. This result is the analogue of Euler's 
Principal Axis Theorem for the sphere, which states that any element 
of 5*0(3) can be written, in some orthonormal basis, as a rotation about 
the z axis. 

The geodesies of H 2 are the hyperbolas obtained by intersecting the 
hyperboloid with planes passing through the origin of the coordinate 
system. For any two distinct points a and b of H 2 , there is a unique 
geodesic that connects them, and the distance between these points is 
given by 



In the framework of Weierstrass's model, the parallels' postulate of 
hyperbolic geometry can be translated as follows. Take a geodesic 7, 
i.e. a hyperbola obtained by intersecting a plane through the origin, 
O, of the coordinate system with the upper sheet, z > 0, of the hyper- 
boloid. This hyperbola has two asymptotes in its plane: the straight 
lines a and 6, intersecting at O. Take a point, P, on the upper sheet 
of the hyperboloid but not on the chosen hyperbola. The plane aP 
produces the geodesic hyperbola a, whereas bP produces (5. These two 
hyperbolas intersect at P. Then a and 7 are parallel geodesies meet- 
ing at infinity along a, while (3 and 7 are parallel geodesies meeting at 
infinity along b. All the hyperbolas between a and (3 (also obtained 
from planes through O) are non-secant with 7. 

Like the Euclidean plane, the abstract Bolyai-Lobachevsky plane 
has no privileged points or geodesies. But the Weierstrass model has 
some convenient points and geodesies, such as the point (0,0, |/t| -1 / 2 ), 
namely the vertex of the sheet z > of the hyperboloid, and the 
geodesies passing through it. The elements of Lor(M 2,1 ) allow us to 
move the geodesies of H 2 to convenient positions, a property that can 
be used to simplify certain arguments. 



(1) 



d(a,b) = (-K^cosh-^/caHb). 
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More detailed introductions to the 2-dimensional Weierstrass model 
can be found in [37] and [69]. The Lorentz group is treated in some 
detail in [4] and [69], but the Principal Axis Theorems for the Lorentz 
group contained in [4] fails to include parabolic rotations, and is there- 
fore incomplete. 

The generalization of the 2-dimensional Weierstrass model to 3 di- 
mensions is straightforward. We consider first the 4-dimensional Min- 
kowski space M 3 ' 1 = (M 4 ,H), where □ is now defined as the Lorentz 
inner product 

a □ b = a w b w + a x b x + a y b y - a z b z , 

with a = (a w ,a x ,a y ,a z ) and b = (b w ,b x ,b y ,b z ). In the Minkowski 
space we embed the z > connected component of the 3-dimensional 
hyperbolic manifold given by the equation 

/r>\ 2,2,2 2 -1 

(2) w + x + y — z — K , 

which models the 3-dimensional hyperbolic space H 3 of constant cur- 
vature k < 0. The distance is given by the same formula Q, where 
a and b are now points in M 4 that lie in the 3-dimensional hyperbolic 
manifold (|2]) with z > 0. 

The next issue to discuss would be that of Lorentzian transformations 
in H 3 . But we postpone this subject, to present it in Section [9] together 
with the isometries of the 3-dimensional sphere. 

5.2. History of the Weirstrass model. The first mathematician 
who mentioned Karl Weierstrass in connection with the hyperboloid 
model of the hyperbolic plane was Wilhelm Killing. In a paper pub- 
lished in 1880, |46| . he used what he called Weierstrass's coordinates 
to describe the "exterior hyperbolic plane" as an "ideal region" of the 
Bolyai-Lobachevsky plane. In 1885, he added that Weierstrass had 
introduced these coordinates, in combination with "numerous appli- 
cations," during a seminar held in 1872, [48], pp. 258-259. We found 
no evidence of any written account of the hyperboloid model for the 
Bolyai-Lobachevsky plane prior to the one Killing gave in a paragraph 
of [IB], p. 260. His remarks might have inspired Richard Faber to 
name this model after Weierstrass and to dedicate a chapter to it in 
[37], pp. 247-278. 

5.3. More geometric background. Since we are interested in the 
motion of point particles in 3-dimensional manifolds, the natural frame- 
work for the study of the 3-dimensional curved n-body problem is the 
Euclidean ambient space, 1R 4 , endowed with a specific inner-product 
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structure, which depends on whether the curvature is positive or nega- 
tive. We therefore aim to continue to set here the problem's geometric 
background in IR 4 . For positive constant curvature, K > 0, the motion 
takes place in a 3-dimensional sphere embedded in the Euclidean space 
IR 4 endowed with the standard dot product, • , i.e. on the manifold 



For negative constant curvature, k < 0, the motion takes place on the 
manifold introduced in the previous subsection, the upper connected 
component of a 3-dimensional hyperboloid of two connected compo- 
nents embedded in the Minkowski space IR 3,1 , i.e. on the manifold 

M 3 K = {(w,x,y,z)\w 2 + x 2 + y 2 - z 2 = n~ x ,z > 0}, 

where IR 3 ' 1 is IR 4 endowed with the Lorentz inner product, □. Generi- 
cally, we will denote these manifolds by 

M 3 = {(w, x, y, z) G IR 4 | w 2 +x 2 +y 2 +az 2 = k~ x , with z > for « < 0}, 
where a is the signum function, 



so M 3 is endowed with the operation 0, meaning • for k > and □ for 
k < 0. 

If R is the radius of the sphere S 3 , then the relationship between 
k > and R is k" 1 = R 2 . As we already mentioned, to have an 
analogue interpretation in the case of negative curvature, M 3 can be 
viewed as a 3-dimensional pseudosphere of imaginary radius iR, such 
that the relationship between k < and iR is k^ 1 = (iR) 2 . 

Let us further define some concepts that will be useful later. 

Definition 1. A great sphere of § 3 is a 2- dimensional sphere of the 
same radius as § 3 . 

Definition [T] implies that the curvature of a great sphere is the same 
as the curvature of S 3 . Great spheres of S 3 are obtained by intersecting 
§ 3 with hyperplanes of IR 4 that pass through the origin of the coordinate 
system. Examples of great spheres are: 



{(w, x, y, z)\w 2 + x 2 + y 2 + z 2 = k 1 }. 



(3) 




for k > 
for k < 0. 



Given the 4-dimensional vectors 



a = (a w , a x , a y , a z ) and b = (b w , b x , b y , b 
we define their inner product as 

(4) a b := a w b w + a x b x + a y b y + cra z b z , 



(5) 



{(w,x,y,z)\x 2 + y 2 + z 2 = k \ w = 0}, 



Relative equilibria in the 3-dimensional curved n-body problem 19 

(6) s l,x = {(w,x,y,z)\w 2 + y 2 + z 2 = kT 1 , x = 0}, 

( 7 ) S% y = {(w,x,y,z)\ w 2 + x 2 + z 2 = k" 1 , y = 0}, 

(8) S 2 KjZ = {(w,x,y,z)\w 2 + x 2 + y 2 = k-\ z = 0}. 

Definition 2. A great circle of a great sphere of is a circle (1- 
dimensional sphere) of the same radius as S 3 K . 

Definition [2] implies that, since the curvature of S>^ is k, the curvature 
of a great circle is n = l/R, where R is the radius of §^ and of the 
great circle. Examples of great circles are: 



(9) 


S 1 

kJ K,WX 


= {(w 




ly 2 


+ z 2 = 


K~ 


1 


W 


= X 


= 0}, 


(10) 


s 1 

K,yz 




x,y,z) 


u; 2 


+ x 2 = 


K 


-1 

? 


y 


= z 


= 0}, 


(11) 


S 1 

ft, MI J/ 


= {(w 






+ z 2 = 


K~ 


-1 

j 


w 


= y 


= 0}, 


(12) 


s 1 

^K,XZ 


= {K 




w 2 


+ x 2 = 


K 


-1 

? 


x 


= z 


= 0}, 


(13) 


s 1 

kJ K,WZ 


= {(«; 




\y 2 


+ z 2 = 


K~ 


-1 


w 


= z 


= 0}, 


(14) 


s 1 


= {(«, 




\y 2 


+ z 2 = 


K~ 


-1 

5 


X 


= y 


= 0}. 



Notice that S* ^ is a great circle for both the great spheres S 2 and 
S 2 x , whereas is a great circle for both the great spheres S 2 and 
S 2 2 . Similar remarks can be made about any of the above pairs of 
great circles. 

Definition 3. Two great circles, G\ and C 2 , of two different great 
spheres of S>^ are called complementary if there is a coordinate system 
wxyz such that 

(15) Ci = S KiWX and C 2 = S^ 2 
or 

(16) d = Sl !Wy and C 2 = S^ xz . 



The representations (15) and (16) of C\ and C 2 give, obviously, all 
the possibilities we have, because, for instance, the representation 



is the same as (15) after we perform a circular permutation of the 
coordinates w, x, y, z. For simplicity, and without loss of generality, we 
will always use representation (15). 
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The pair C\ and C2 of complementary circles of S> 3 , for instance, 
forms, in topological terms, a Hopf link in a Hopf fibration, which is 
the map 

% : S 3 — > S> 2 , h(w, x, y, z) = {w 2 + x 2 — y 2 — z 2 , 2{wz + xy), 2(xz — wy)) 

that takes circles of § 3 to points of § 2 , [12], [59]. In particular, H, 
takes S\ wx to (1,0,0) and S}^ to (—1,0,0). Using the stereographic 
projection, it can be shown that the circles C\ and C2 are linked (like 
any adjacent rings in a chain), hence the name of the pair, [59J. Of 
course, these properties can be expressed in terms of any curvature 
k > 0. Hopf fibrations have important physical applications in fields 
such as rigid body mechanics, |60j, quantum information theory, |64j . 
and magnetic monopoles, [65] . As we will see later, they are also useful 
in celestial mechanics via the curved n-body problem. 

We will show in the next section that the distance between two points 
lying on complementary great circles is independent of their position. 
This remarkable geometric property turns out to be even more surpris- 
ing from the dynamical point of view. Indeed, given the fact that the 
distance between 2 complementary great circles is constant, the magni- 
tude of the gravitational interaction (but not the direction of the force) 
between a body lying on a great circle and a body lying on the com- 
plementary great circle is the same, no matter where the bodies are on 
their respective circles. This simple observation will help us construct 
some interesting, nonintuitive classes of solutions of the curved n-body 
problem. 

In analogy with great spheres of S 3 , we define great hyperboloids of 
H 3 as follows. 

Definition 4. A great hyperboloid o/HI 3 is a 2- dimensional hyperboloid 
of the same curvature as H 3 . 

Great hyperboloids of H 3 are obtained by intersecting H 3 with hy- 
perplanes of M 4 that pass through the origin of the coordinate system. 
Examples of great hyperboloids are: 

(17) U% w = {(w, x, y, z) I x 2 + y 2 - z 2 = k~\ w = 0}, 

(18) H 2 jX = {(w, x, y, z)\w 2 + y 2 -z 2 = k~\ x = 0}, 

(19) U 2 Kty = {(w,x,y,z)\w 2 + x 2 - z 2 = k~\ y = 0}. 

It is interesting to remark that great spheres and great hyperboloids 
are totally geodesic surfaces in S 3 and H 3 , respectively. Recall that if 
through a given point of a Riemannian manifold (such as § 3 and H 3 ) 
we take the various geodesies of that manifold tangent at this point to 
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the same plane element, we obtain a geodesic surface. A surface that 
is geodesic at each of its points is called totally geodesic. 

5.4. The metric. An important preparatory issue lies with introduc- 
ing the metric used on the manifolds S> 3 and H 3 , which, given the 
definition of the inner products, we naturally take as 

fiT 1 / 2 cos -1 (/ca • b), k > 

(20) 4(a,b) := { |a-b|, k = 

-/^T^cosh^OaHb), k < 0, 

where the vertical bars denote the standard Euclidean norm. 

When k — > 0, with either n > or n < 0, then R — >• oo, where i? 
represents the radius of the sphere § 3 or the real factor in the expression 
iR of the imaginary radius of the pseudosphere H 3 . As R — » oo, both 
§ 3 and H 3 become M 3 , and the vectors a and b become parallel, so the 
distance between them is given by the Euclidean distance, as indicated 



in (20). Therefore, in a way, d is a continuous function of k when the 
manifolds § 3 and H 3 are pushed to infinity. 

To get more insight into the fact that the metric in § 3 and H 3 be- 
comes the Euclidean metric in M 3 when k — > 0, let us use the stereo- 
graphic projection, which takes the points of coordinates (w, x, y, z) G 
M 3 , where 

M 3 = {(w, x, y, z) G M 4 | w 2 +x 2 +y 2 +az 2 = with z > for k < 0}, 

to the points of coordinates (W, X, Y) of the 3-dimensional hyperplane 
z = through the bijective transformation 

(21) W=-^-, X=^-, Y RV 



R — az R — az R — az 
which has the inverse 

2R 2 W 2R 2 X 

W = r— ; — zr, X 



y 



R 2 + aW 2 + oX 2 + oY 2 ' R 2 + aW 2 + oX 2 + oY 2 ' 
2R 2 Y R(W 2 + X 2 + Y 2 - aR 2 ) 



R? + aW 2 + aX 2 + aY 2 ' ' R 2 + aW 2 + oX 2 + oY 2 
From the geometric point of view, the correspondence between a point 
of M 3 and a point of the hyperplane z = is made via a straight line 
through the point (0, 0, 0, aR), called the north pole, for both k > 
and n < 0. 

For k > 0, the projection is the Euclidean space IR 3 , whereas for 
k < it is the 3-dimensional solid Poincare ball of radius k^ 1 ^ 2 . The 



22 



Florin Diacu 



metric in coordinates (W, X, Y) is given by 

AR\dW 2 + dX 2 + dY 2 ) 



ds< 



(R 2 + aW 2 + aX 2 + aY 2 ) 2 ' 

which can be obtained by substituting the inverse of the stereographic 
projection into the metric 

ds 2 = dw 2 + dx 2 + dy 2 + adz 2 . 

The stereographic projection is conformal, but it's neither isometric 
nor area preserving. Therefore we cannot expect to recover the exact 
Euclidean metric when k — > 0, i.e. when R — > oo, but hope, never- 
theless, to obtain some expression that resembles it. Indeed, we can 
divide the numerator and denominator of the right hand side of the 
above metric by i? 4 and write it after simplification as 

2 _ 4{dW 2 + dX 2 + dY 2 ) 

S ~ (1 + aW 2 /R 2 + aX 2 /R 2 + aY 2 /R 2 ) 2 ' 

When R — > oo, we have 

ds 2 = A(dW 2 + dX 2 + dY 2 ) 1 
which is the Euclidean metric of R 3 up to a constant factor. 



Remark 1. From (20) we can conclude that if, for k > 0, C\ and C2 
are two complementary great circles, as described in Definition [3j and 
a G C\, b £ C2, then the distance between a and b is 

4(a,b) = fiT 1/2 7r/2. 

This fact shows that the distance between two complementary circles 
is constant. 

Since, to derive the equations of motion, we will apply a variational 
principle, we need to extend the distance from the 3-dimensional man- 
ifolds of constant curvature S 3 . and H 3 . to the 4-dimensional ambient 
space. We therefore redefine the distance between a and b as 

Ka b k > 

(22) r/ K (a.b):: { I, k = 

KaHb K < 0. 




v^aEWrebHb ' 



Notice that on § 3 we have y/nsL ■ a = y/nh h = 1 and on H 3 we have 
v^a □ a = \J Kh □ b = 1, which means that the new distance reduces 
to the previously defined distance on the corresponding 3-dimensional 
manifolds of constant curvature. 
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5.5. Unified trigonometry. Following the work of Carinena, Ranada, 
and Santander, [9] , we will further define the trigonometric ^-functions, 
which unify circular and hyperbolic trigonometry. The reason for this 
step is to obtain the equations of motion of the curved n-body problem 
in both constant positive and constant negative curvature spaces. We 
define the K-sine, sn K , as 

{k~ 1/2 smn 1/2 x if k > 
x if k = 
(-k)- 1 / 2 sinh(- Ky/ 2 x if k<0, 

the K-cosine, csn K , as 

{cos k 1 I 2 x if k > 
1 if k = 
cosh(-K) 1/2 a; if k < 0, 

as well as the K-tangent, tn K , and /t-cotangent, ctn K , as 

, ( v sn K (x) csn K (x) 
tn K (x) := -— and ctn K (x) := -— , 



csn K (x) sn K {x) 

respectively. The entire trigonometry can be rewritten in this unified 
context, but the only identity we will further need is the fundamental 
formula 

(23) Ksn 2 K (x) + csn 2 K (x) = l. 

Notice that all the above trigonometric K-functions are continuous 
with respect to k. In the above formulation of the unified trigonometric 
K-functions, we assigned no particular meaning to the real parameter 
k. In what follows, however, k will represent the constant curvature of 
a 3-dimensional manifold. Therefore, with this notation, the distance 



(20) on the manifold can be written as 

4(a,b) = (^-^[(wj^aeb] 
for any a, b <G and k^O. 



6. Equations of motion 

The main purpose of this section is to derive the equations of motion 
of the curved n-body problem on the 3-dimensional manifolds M^. To 
achieve this goal, we will introduce the curved potential function, which 
also represents the potential energy of the particle system, present and 
apply Euler's formula for homogeneous functions to the curved po- 
tential function, and introduce the variational method of constrained 
Lagrangian dynamics. After we derive the equations of motion of the 
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curved n-body problem, we will show that they can be put in Hamil- 
tonian form. 

6.1. The potential. Since the classical Newtonian equations of the 
n-body problem are expressed in terms of a potential function, our 
next goal is to define a potential function that extends to spaces of 
constant curvature and reduces to the classical potential function in 
the Euclidean case, i.e. when k = 0. 

Consider the point particles (also called point masses or bodies) of 
masses mi,m2 ! ...,m n > in IR 4 , for k > 0, and in IR 3,1 , for k < 
0, whose positions are given by the vectors = (wi,Xi,yi, z^), i = 
1,2, ... ,n. Let q = (qi, q 2 , . . . , q n ) be the configuration of the system 
and p = (pi,p 2 , . . . ,p n ), with p^ = m^, i = 1,2,..., n, be the 
momentum of the system. We define the gradient operator with respect 
to the vector q, as 

V qi := (d Wi ,d x .,dy.,ad Zi ). 

From now on we will rescale the units such that the gravitational 
constant G is 1. We thus define the potential of the curved n-body 
problem, which we will call the curved potential, as the function — U K , 
where 

j n n 

1=1 j = l,j^i 

stands for the curved force function, and q = (qi, . . . , q n ) is the config- 
uration of the system. Notice that, for k = 0, we have ctn (c/ (qj, q^)) = 
|Oi ~~ which means that the curved potential becomes the classi- 

cal Newtonian potential in the Euclidean case. Moreover, U K — > Uq as 
k — y 0, whether through positive or negative values of n. Nevertheless, 
we cannot claim that U K is continuous with respect to k at in the 
usual way, but rather in a degenerate sense. Indeed, when k — > 0, the 
manifold M 3 . is pushed to infinity, so the above continuity with respect 
to k must be understood in this restricted way. In fact, as we will 
further see, the curved force function, U K , is homogeneous of degree 0, 
whereas the Newtonian potential, U , defined in the Euclidean space, 
is a homogeneous function of degree —1. 

Now that we defined a potential that satisfies the basic limit condi- 
tion we required of any extension of the n-body problem beyond the 
Euclidean space, we emphasize that it also satisfies the basic properties 
the classical Newtonian potential fulfills in the case of the Kepler prob- 
lem, as mentioned in the Introduction: it obeys Bertrand's property, 
according to which every bounded orbit is closed, and is a solution 
of the Laplace-Beltrami equation, [49J, the natural generalization of 
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Laplace's equation to Riemannian and pseudo-Riemannian manifolds. 
These properties ensure that the cotangent potential provides us with 
a natural extension of Newton's gravitational law to spaces of constant 
curvature. 

Let us now focus on the 7^ 0. A straightforward computation, 

which uses the fundamental formula (23), shows that 

(24) W= E E ^^M, I#0 , 



2 

i=i j=i,j& 1 1 (j — (j 



KqiOqj 



V^qiOqi y/ tqj 0qj 



an expression that is equivalent to 
(25) 



In fact, we could simplify U K even more by recalling that kc^ © q« = 
1, % = 1,2, ...,n. But since we still need to compute VU K , which 
means differentiating U K , we we will not make that simplification here. 



6.2. Euler's formula for homogeneous functions. In 1755, Leon- 
hard Euler proved a beautiful formula related to homogeneous func- 
tions, [36]. We will further present it and show how it applies to the 
curved potential. 

Definition 5. A function F: M. m — > R is called homogeneous of degree 
a G 1R if for all rj ^ and q G M m ; we have 

F( V q) = V a F(q). 

Euler's formula shows that, for any homogeneous function of degree 
a G K, we have 

q . VF(q) = aF(q) 

for all q G R m . 

Notice that U K (r]q) = U K (q) = rfU K (q) for any 77 7^ 0, which means 
that the curved potential is a homogeneous function of degree zero. 
With our notations, we have m = 3n, therefore Euler's formula can be 
written as 

q© VF(q) = aF(q). 
Since a = for U K with k ^ 0, we conclude that 



(26) 



q©W K (q) = 0. 
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We can also write the curved force function as U K (q) = \ Y17=i ^(O.*)) 
where 

n T(ljT(lj\<JK) ' - — , 

y V V^qiOqiA/^qjOqj / 

are also homogeneous functions of degree 0. Applying Euler's formula 
for functions F : 1R 3 — > R, we obtain that q* V qi ?7*(q) = 0. Then 
using the identity V qi C/ re (q) = V qi Z7*(qi), we can conclude that 

(27) q^OVq^q) = 0, i = l,...,n. 

6.3. Constrained Lagrangian dynamics. To obtain the equations 
of motion of the curved ra-body problem, we will use the classical vari- 
ational theory of constrained Lagrangian dynamics, [39] . According to 
this theory, let 

L = T — V 

be the Lagrangian of a system of n particles constrained to move on a 
manifold, where T is the kinetic energy and V is the potential energy 
of the system. If the positions and the velocities of the particles are 
given by the vectors q«,qi, i = 1,2, ...,n, and the constraints are 
characterized by the equations fa = 0, i = l,2,...,n, respectively, 
then the motion is described by the Euler-Lagrange equations with 
constraints, 

where A*, i = 1,2, ... ,n, are the Lagrange multipliers. To obtain 
this formula we have to assume that the distance is defined in the 
entire ambient space. Using this classical result, we can now derive the 
equations of motion of the curved n-body problem. 

6.4. Derivation of the equations of motion. In our case, the po- 



tential energy is V = —U K , given by the curved force function (24), 
and we define the kinetic energy of the system of particles as 

1 - 

T K (q,q) := - rriife qi)(«qi q;). 

i=i 

The reason for introducing the factors nc[i q^ = 1, % = 1,2, ... ,n, 
into the definition of the kinetic energy will become clear in Subsection 



6.5 Then the Lagrangian of the curved n-body system has the form 
L K (q, q) = T K (q, q) + U K (q). 
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So, according to the theory of constrained Lagrangian dynamics dis- 
cussed above, which requires the use of a distance defined in the am- 



bient space — a condition we satisfied when we produced formula (22) 
the equations of motion are 



(29) 



d_ ( dL^ 
dt \ dqi 



9Lk X(t) df * 



o, 



where /* = qi — k' 1 is the function that characterizes the con- 
straints /* = 0, i = 1,2, ... ,n. Each constraint keeps the body of 
mass rrii on the surface of constant curvature k, and X l K is the Lagrange 
multiplier corresponding to the same body. Since qjOq^ = k -1 implies 
that q« q« = 0, it follows that 



d_ ( dLn 
dt \ dqi 



rriiq^Kqi qi) + 2mi(Kqi qi) = m*qj 



This relationship, together with 
dL K 



dqi 



miK(qi qi)qi + V^U K (q) 



implies that equations (29) are equivalent to 



(30) miqi-miK(qi0qi)qi-V qi £/ K (q)-2A l K (i)qi = 0, i = l,...,n. 
To determine notice that = /* = 2qi qi + 2(qj qi), so 

(31) qi q % = -qi qi- 



Let us also remark that ©-multiplying equations (30) by qi and using 



Euler's formula (27), we obtain that 



m<(qi © qi) - rrn(qi qi) - qi V qi ?7 K (q) = 2A l K qi © qi = 2k \* k , 



which, via (31), implies that A* = — nm^qi qi). Substituting these 



values of the Lagrange multipliers into equations (30), the equations of 
motion and their constraints become 

(32) miqi = V qi f/ K (q) - miK(qi © qi)qi, q i Oq i = K~ 1 , K ^ 0, 

i — I, . . . , n. 
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The qj-gradient of the curved force function, obtained from (24), has 
the form 



qi©q 7 ' 
■=q;Qq, 



1i 



(33) V q ,£«q) 



E 



^Kqi 0qi -y/ftq.,- q,, 



(7 



(7 



V^qiOq; y/ nqjOqj 



3/2 ' 



which is equivalent to 
(34) 

v j j ( , — V* m i m j\ K \ 3/2 ( K< ij qj)[( K q» Q qpg? - (^q» qj)q» 
v qi «w - 2^ [(j(/tq . Q q . )(req . ^ _ ^ qj)2]3/2 

Using the fact that /tq^ q^ = 1, we can write this gradient as 

(35) v q! f/ K (q) = > j "--y ^ v "*: 2 r 3 /? / '" J . * ^ °- 



[cr - <r («qj qj) ] 



E : 

Sometimes we can use the simpler form ([35j) of the gradient, but 
whenever we need to exploit the homogeneity of the gradient or have 



to differentiate it, we must use its original form (34). Thus equations 



(32) and (34) describe the n-body problem on surfaces of constant 



curvature for k ^ 0. Though more complicated than the equations of 



motion Newton derived for the Euclidean space, system (32) is simple 
enough to allow an analytic approach. 

6.5. Hamiltonian formulation. It is always desirable to place any 
new problem into a more general theory. In our case, like in the classical 
n-body problem, the theory of Hamiltonian systems turns out to be the 
suitable framework. 

The Hamiltonian function describing the motion of the n-body prob- 
lem in spaces of constant curvature is 

# K (q,p)=T K (q,p)-£/ K (q). 

Then the Hamiltonian form of the equations of motion is given by the 
system with constraints 

q, = V Pi H K (q, p) = m^p, 
(36) { pi = -V qi # K (q, p) = V qi C/ K (q) - m^/^P* © P*)q*> 



q* 0q; 



q» Pi = o, k^o, 



1,2 n. 



The configuration space is the manifold (M^) n , where, recall, 
represents E> 3 K or H^. Then (q, p) G T* 



Kj n , where T" 



is the 
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cotangent bundle, which represents the phase space. The constraints 
K<\i q« = 1, qj pj = 0, i = 1,2,. . . ,n, which keep the bodies 
on the manifold and show that the position vectors and the momenta 
of each body are orthogonal to each other, reduce the 8n-dimensional 



system (36) by 2n dimensions. So, before taking into consideration 
the first integrals of motion, which we will compute in Section [7| the 
phase space has dimension 6n, as it should, given the fact that we are 
studying the motion of n bodies in a 3-dimensional space. 

6.6. Invariance of great spheres and great hyperboloids. In the 

Euclidean case, planes are invariant sets for the equations of motion. 
In other words, if the initial positions and momenta are contained in 
a plane, the motion takes place in that plane for all time for which 
the solution is defined. We can now prove the natural analogue of this 
result for the curved n-body problem. More precisely, we will show 
that, in and H^, the motion can take place on 2-dimensional great 
spheres and 2-dimensional great hyperboloids, respectively, if we assign 
suitable initial positions and momenta. 

Proposition 1. Let n > 2 and consider the point particles of masses 
mi,m 2 , ... ,m n > in Mj*. Assume that M 2 is any 2-dimensional 
submanifold of having the same curvature, i. e. a great sphere for 
k > or a great hyperboloid for k < 0. Then, for any nonsingular 
initial conditions (q(0),p(0)) G (M^) n x (T(M 2 ) n , where x denotes 
the cartesian product of two sets and T(M. 2 K ) is the tangent space of 
M 2 ., the motion takes place in M 2 K . 

Proof. Without loss of generality, it is enough to prove the result for 
M l,wi where 

M l,w ■= {O; z> V, z)\x 2 + y 2 + az 2 = k -1 , w = 0} 

is the great 2-dimensional sphere S^, for k > 0, and the great 2- 
dimensional hyperboloid H^, for k < 0, both of which we defined in 
Subsection 6.6 as ^ and (17), respectively. Indeed, we can obviusly 



restrict to this case since any great sphere or great hyperboloid can be 
reduced to it by a suitable change of coordinates. 

Let us denote the coordinates and the momenta of the bodies m i? i 
1,2, ... ,n, by 

qi = {w i ,x i ,y il Zi) and p 4 = {r h s i5 u h Vi), £ = l,2,...,n, 
which when restricted to M^ to and T(M^ U ,), respectively, become 
qi = (0,Xi,yi,Zi) and p; = (0, s h u h i = 1, 2, . . . , n. 
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Relative to the first component, w, the equations of motion (36) have 
the form 

Wi = m~ l ri, 

( qi0qi = « _1 , qiOpi = 0, k^O, £ = 1,2, ...,n. 

For our purposes, we can view this first-order system of differential 
equations as linear in the variables Wi, j-j, i = 1, 2, . . . , n. But on w , 
the initial conditions are Wi(0) = ^(0) = 0, £ = 1,2, . . . ,n, therefore 
Wi(t) = Ti{t) = 0, £ — 1, 2, . . . , n, for all t for which the corresponding 
solutions are defined. Consequently, for initial conditions (q(0), p(0)) e 
x T(M^), the motion is confined to w , a remark that completes 
the proof. □ 

7. First integrals 

In this section we will determine the first integrals of the equations of 
motion. These integrals lie at the foundation of the reduction method, 
which played an important role in the theory of differential equations 
ever since mathematicians discovered the existence of functions that 
remain constant along solutions. The classical ra-body problem in IR 3 
has 10 first integrals that are algebraic with respect to q and p, known 
already to Lagrange in the mid 18th century, [S2]- In 1887, Heinrich 
Bruns proved that there are no other first integrals, algebraic with 
respect to q and p, [7j. Let us now find the first integrals of the curved 
n-body problem. 

7.1. The integral of energy. The Hamiltonian function provides the 
integral of energy, 

(37) H K (q,p) = h, 



where h is the energy constant. Indeed, ©-multiplying equations (32) 
by q,, we obtain 

n n n , 

^miqiOqi = ^[V qi £/ K (q)]Oqi-^mi/t(qiOqj)qiOqi = —U K (q(t)). 

i=l i=l i=l 



Then equation (37) follows by integrating the first and last term in the 
above equation. 

Unlike in the classical Euclidean case, there are no integrals of the 
centre of mass and the linear momentum. This is not surprising, giv- 
ing the fact that n-body problems obtained by discretizing Einstein's 
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field equations lack these integrals too, [31], [38], [50] , [51]. Their ab- 
sence, however, complicates the study of our problem since many of 
the standard methods don't apply anymore. 

One could, of course, define some artificial centre of mass for the 
particle system, but this move would be to no avail. Indeed, forces do 
not cancel each other at such a point, as it happens in the Euclidean 
case, so no advantage can be gained from introducing this concept. 

7.2. The integrals of the total angular momentum. As we will 



show below, equations (36) have six angular momentum integrals. To 
prove their existence, we need to introduce the concept of bivector, 
which generalizes the idea of vector. If a scalar has dimension 0, and a 
vector has dimension 1, then a bivector has dimension 2. A bivector is 
constructed with the help of the wedge product a A b, defined below, 
of two vectors a and b. Its magnitude can be intuitively understood 
as the oriented area of the parallelogram with edges a and b. The 
wedge product lies in a vector space different from that of the vectors 
it is generated from. The space of bivectors together with the wedge 
product is called a Grassmann algebra. 
To make these concepts precise, let 

e, = (1, 0, 0, 0), e x = (0, 1, 0, 0), e y = (0, 0, 1, 0), e z = (0, 0, 0, 1) 
denote the elements of the canonical basis of IR 4 , and consider the 



vectors u = (u w ,u x ,u y ,u z ) and v = (v w ,v x ,v y ,v z ). We define the 
wedge product (also called outer product or exterior product) of u and 
v of IR 4 as 

u A v := (u w v x - u x v w )e w A e x + [u w v y - u y v w )e w A e y + 

(38) iu w v z - u z v w )e w Ae z + (u x v y - u y v x )e x A e y + 

(u x v z - u z v x )e x Ae z + (u y v z - u z v y )e y A e z , 

where e w A e x ,e w A e y ,e w A e z , e x A e y ,e x A e z ,e y A e z represent the 
bivectors that form a canonical basis of the exterior Grassmann algebra 
over IR 4 (for more details, see, e.g., |32j). In IR 3 , the exterior product is 
equivalent with the cross product. 

Let us define YH=i m i c li A qi to be the total angular momentum of 
the particles of masses m 1 ,m 2 ,...,m n > in IR 4 . We will further 
show that the total angular momentum is conserved for the equations 
of motion, i.e. 



(39) ^2 m i^i A 4i 



=i 



where c Cwx^w A Q x ~\~ c wy Q w A Q y ~\~ c wz g w A g z -\- c xy G x A s y ~\~ c xz s x A g z -\- 
c yz e y A e z , with the coefficients c wx ,c wy ,c wz ,c xy ,c xz ,c yz G IR. Indeed, 
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relations (39) follow by integrating the identity formed by the first and 



last term of the equations 



n n 



(40) j2 m ^ A ^ = Yl E 



i=i 

n 

E 



mivrij | K| 3 / 2 qj A 
[(j - a(Kq, qj) 2 ] 3 / 2 



E 



mimj\K,\ 3 / 2 (Kqi qj) 
[a - a(nqi qj) 2 ] 3 / 2 



m^q; qi 



qi A qi = 0, 



obtained after A-multiplying the equations of motion (32) by q^. The 



last of the above identities follows from the skew-symmetry of the A 
operation and the fact that q^ A q^ = 0, i — 1, . . . , n. 



On components, the 6 integrals in (39) can be written as 



(41) ^m^WjXj - WiXi) = c wx , E 

i=l i=l 



m^Wiyi - Wiyi) = c wy , 



(42) E 
i=i 

n 

(43) E 



i=l i=l 
n n 

Tfli(XiZi XiZi) C X zi ^ TTlj {%)jZj ijiZij Cyz^ 

i=l i=l 



The physical interpretation of these six integrals is related to the ge- 
ometry of IR 4 . The coordinate axes Ow, Ox, Oy, and Oz determine 
six orthogonal planes, wx,wy,wz,xy,xz, and yz. We call them basis 
planes, since they correspond to the exterior products e w Ae x , e w A e y , 
e 
e 



A e 2 , A e„, e x A e 2 , and e„ A e 2 , respectively, of the basis vectors 



Of 



Then the constants c 



mi! ^-"wyi '-'■wzi '-'xyi *-xzi ^yz 



Ciiz mea- 



sure the rotation of an orbit relative to a point in the plane their indices 
define. This point is the same for all 6 basis planes, namely the origin 
of the coordinate system. 

To clarify this interpretation of rotations in IR 4 , let us point out that, 
in IR 3 , rotation is understood as a motion around a pointwise invariant 
axis orthogonal to a basis plane, which the rotation leaves globally (not 
pointwise) invariant. In IR 4 , there are infinitely many axes orthogonal 
to this plane, and the angular momentum is the same for them all, since 



each equation of (43) depends only on the 2 coordinates of the plane 



and the corresponding velocities. It is, therefore, more convenient to 
think of these rotations in IR 4 as taking place around a point in a plane, 
in spite of the fact that the rotation moves points outside the plane too. 

Whatever sense of rotation a scalar constant of the angular momen- 
tum determines, the opposite sign indicates the opposite sense. A zero 
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scalar constant means that there is no rotation relative to the origin in 
that particular plane. 

Using the constraints Kqi q 8 = 1, % = 1,2, ... ,n, we can write 



system (36) as 
(44) 



m j \K\"'"[q j - (nqi q^q,] 



3/2 



. , [a-a^Oq,) 2 ] 3 / 2 



[niii qi)qi, % = 1, 2, . . . , n, 



which is the form of the equations of motion we will mostly use in 
this paper. The sums on the right hand side of the above equations 
represent the ith gradient of the potential after the constraints are 
taken into account. 

If we regard these equations as a first order system with constraints, 

{q, = m~V', 
Pi - 2^=1 [ff - ff{Kqi0q ^ 2 ]3/2 — («q* © Wh 
q, q, = 1, q, pi = 0, i = 1, 2, . . . , n, 

the dimension of the phase space, after taking into account the integrals 
of motion described above, is 6n — 7. 

8. Singularities 

Before we begin the study of relative equilibria, it is important 
to know whether there exist impossible configurations of the bod- 
ies. The answer is positive, and these configurations occur when sys- 



tem (44) encounters singularities, i.e. if at least one denominator in 



the sum on the right hand sides of the equations occurring in sys- 



tem (44) vanishes. In other words, a configuration is singular when 



(/■cqi qj) =1, i, j = 1, 2, . . . , n, i ^ j, which is the same as saying 
that nqi © q^ = 1 or nqi Q qj = —1, i,j = l,2,...,n, i ^ j. The fol- 
lowing result shows that the former case corresponds to collisions, i.e. to 
configurations for which at least two bodies have identical coordinates, 
whereas the latter case occurs in § 3 , but not in H 3 , and corresponds 
to antipodal configurations, i.e. when at least two bodies have coordi- 
nates of opposite signs. These are impossible initial configurations all 
the solutions we will introduce in this paper must avoid. 

Proposition 2. (Collision and antipodal configurations) Con- 
sider the 3-dimensional curved n-body problem, n > 2, with masses 
m 1 ,m 2 , . . . , m n > 0. Then, in S> 3 ; if there are i,j G {1, 2, ... , n}, i ^ j , 
such that nqi • qj = 1, the bodies and rrij form a collision configu- 
ration, and if Kqi q^ = — 1, they form an antipodal configuration. In 
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if there are i,j G {1, 2, ... , n}, i ^ j , such that /tq, □ qj = 1, the 
bodies rrii and rrij form a collision configuration, whereas configurations 
with nc[i □ qj = — 1 don't exist. 

Proof. Let us first prove the implication related to collision configura- 
tions for k > 0. Assume that there exist i,j G {1,2, ... ,n}, i ^ j, 
such that nc[i q, = 1, relationship that can be written as 

(46) K(wiWj + XiXj + UiUj + ZiZj) = 1. 

But since the bodies are on the coordinates satisfy the conditions 

I 2 , 2 , 2 i 2\ / 2 i 2 i 2 i 2\ 1 

K[Wi +x i +y i +z i ) = K{ Wj + Xj + y j + z-) = 1. 
Consequently we can write that 

(wiWj + x^j + yiyj + ZiZjf = (w^ + x\ + y\ + crzf)(w] + x) + y) + zj), 

which is the equality case of Cauchy's inequality. Therefore there is a 
constant r / such that Wj = tw%, Xj = TXi, yj = ryi, and Zj = TZi. 
Substituting these values in equation (46), we obtain that 

(47) Kr(wiWj + XiXj + ytyj + z { z 3 ) = I. 



Comparing (|46|) and (|47|), it follows that r = 1, so Wi — Wj, Xi — x 3 , 



yi = yj, and Zi = Zj, therefore the bodies m; and rrij form a collision 
configuration. 

The proof of the implication related to antipodal configurations for 
k > is very similar. Instead of relation (46), we have 

K,(wiWj + XiXj + yiyj + z t Zj) = -1. 

Then, following the above steps, we obtain that r = —1, so Wi = —Wj, 
Xi = —Xj, yi = —yj, and Zi = —Zj, therefore the bodies m 8 and re- 
form an antipodal configuration. 

Let us now prove the implication related to collision configurations 
in the case k < 0. Assume that there exist i,j G {1,2, .. . ,n}, i ^ j, 
such that K<\i q,- = 1, relationship that can be written as 

K(WiWj + XiXj + ViVj - ZiZj) = 1, 

which is equivalent to 

(48) WiWj + XiXj + yiyj — fiT 1 = ZiZj. 

But since the bodies are on H^, the coordinates satisfy the conditions 

/ 2 i 2 , 2 2\ / 2 i 2 , 2 2\ i 

k[w\ + x i +y i - z i ) = K[Wj + Xj + y j - Zj) = 1, 
which are equivalent to 

(49) wj + x\ + yf - kT x = zi and w) + x] + y) - k~ 1 = z). 
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From ( 48 ) and ( 49 ) , we can conclude that 

(wiWj + XiXj + yiyj - fiT 1 ) 2 = (w 2 + x 2 + y 2 - kT x ){w) + x 2 + y 2 - fiT 1 ), 

which is equivalent to 

(wiXj - WjXif + (Wiyj - wjyi) 2 + (x^j - Xjvi) 2 

-K -1 [(u>i - wj) 2 + (xi - xj) 2 + (yi - yj) 2 } = 0. 
Since — k,^ 1 > 0, it follows from the above relation that Wi = Wj, = 
Xj, and yi = yj. Then relation (49) implies that z 2 = z 2 . But since for 
k < all z coordinates are positive, we can conclude that Zq^ — Z j j SO 
the bodies rrii and rrij form a collision configuration. 

For k < 0, we now finally prove the non-existence of configurations 
with nqi E\qj = —1, i, j £ {1,2,..., n}, i ^ j. Let us assume that they 
exist. Then 



(50) WiWj + x^j + yiVj = ZiZ 



j 



K 



(51) ^ + Xj + y { = z i + k and Wj + Xj + yj = Zj + k 
According to Cauchy's inequality, we have 

(wiWj + xiXj + yiyjf < (w 2 + x 2 + y 2 )(w 2 + x 2 + y 2 ), 



which, using (50) and (51), becomes 



(Z.;Z 



K' 1 ) 2 < (Z 2 + K- X ){Z 2 + K 



iii j _^ ^j i iii J\^j 

Since — k^ 1 > 0, this inequality takes the form 

+ zj) 2 < 0, 

which is impossible because z iy Zj > 0, a contradiction that completes 
the proof. □ 

It is easy to construct solutions ending in collisions. Place, for in- 
stance, 3 bodies of equal masses at the vertices of an Euclidean equilat- 
eral triangle, not lying on the same geodesic if k > 0, and release them 
with zero initial velocities. The bodies will end up in a triple collision. 
(If, for k > 0, the bodies lie initially on a geodesic and have zero ini- 
tial velocities, they won't move in time, a situation that corresponds 
to a fixed-point solution for the equations of motion, as we will show 
in Section 12.) The question of whether there exist solutions ending 
in antipodal or hybrid (collision-antipodal) singularities is harder, and 
it was treated in [2S] and [2H]- But since we are not touching on this 
subject when dealing with relative equilibria, we will not discuss it fur- 
ther. All we need to worry about in this paper is to avoid placing the 
bodies at singular initial configurations, i.e. at collisions for k ^ or 
at antipodal positions for k > 0. 
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Part 3. ISOMETRIES AND RELATIVE EQUILIBRIA 

9. Isometric rotations in S> 3 and H 3 

This section introduces the isometric rotations in § 3 and H 3 , since 
they play an essential role in defining the relative equilibria of the 
curved n-body problem. There are many ways to express these rota- 
tions, but their matrix representation will suit our goals best, as it did 



in Subsection 5.1 for the 2-dimensional Weierstrass model of hyperbolic 
geometry. 

For k > 0, the isometric transformations of S 3 . are given by the el- 
ements of the Lie group SO (4) of M 4 that keep § 3 invariant. They 
consist of all orthogonal transformations of the Lie group 0(4) repre- 
sented by matrices of determinant 1, and have the form PAP -1 , with 
P e 50(4) and 

/ cosfl -sin9 \ 
sin 9 cos# 
cos — sin 
\ sin0 cos0 / 

where 9, <f> 6 R. We will call these rotations ^-positive elliptic-elliptic 
if 9 7^ and 7^ 0, and K-positive elliptic if 9 7^ and = (or 9 = 
and 7^ 0, which is a possibility we will never discuss since it perfectly 
resembles the previous one). When 9 = = 0, A is the identity matrix, 
so no rotation takes place. The above description is a generalization to 
§ 3 of Euler's Fixed Axis Theorem for As we will next explain, the 
reference to a fixed axis is, from the geometric point of view, far from 
suggestive in M 4 . 



(52) A 



The form of the matrix A given by (52) shows that the K-positive 
elliptic-elliptic transformations have two circular rotations, one relative 
to the origin of the coordinate system in the plane wx and the other 
relative to the same point in the plane yz. In this case, the only fixed 
point in M 4 is the origin of the coordinate system. The K-positive 
elliptic transformations have a single rotation around the origin of the 
coordinate system that leaves infinitely many axes (in fact, an entire 
plan) of 1R 4 pointwise fixed. 

For k < 0, the isometric transformations of W\ are given by the 
elements of the Lorentz group Lor (3, 1), a Lie group in the Minkowski 
space M 3 ' 1 . Lor(3, 1) is formed by all orthogonal transformations of 
determinant 1 that keep H 3 invariant. The elements of this group are 
the ^-negative elliptic, K-negative hyperbolic, and K-negative elliptic- 
hyperbolic transformations, on one hand, and the K-negative parabolic 
transformations, on the other hand; they can be represented as matrices 
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of the form PBP- 1 and PCP~ 



(53) 



B 



respectively, with P e Lor(3, 1), 
\ 



(54) 



C 



/ cos 9 — sin 9 

sin 9 cos 9 

coshs sinhs 

\ sinh s cosh s J 



( 1 \ 

o i -i i 

i l-e 2 /2 e 2 /2 

v o i -e/2 i+e/2 j 

The K-negative elliptic, k- 



where 9,s,£ are some fixed values in 
negative hyperbolic, and re-negative elliptic-hyperbolic transformations 
correspond to 9 ^ and s = 0, to 9 = and s ^ 0, and to 9 ^ and 
s 7^ 0, respectively. The above description is a generalization to HEof 



the Fixed Axis Theorem for Ef^, which we presented in Subsection 5.1 



Again, as in the case of the group SO(A) } the reference to a fixed axis 
has no real geometric meaning in R 3 ' 1 . 

Indeed, from the geometric point of view, the K-negative elliptic 
transformations of 1R 3,1 are similar to their counterpart, K-positive ellip- 
tic transformations, in IR 4 , namely they have a single circular rotation 
around the origin of the coordinate system that leaves infinitely many 
axes of IR 3 ' 1 pointwise invariant. The K-negative hyperbolic transfor- 
mations correspond to a single hyperbolic rotation around the origin 
of the coordinate system that also leaves infinitely many axes of IR 3,1 
pointwise invariant. The k- negative elliptic- hyperbolic transformations 
have two rotations, a circular one about the origin of the coordinate 
system, relative to the wx plane, and a hyperbolic one about the origin 
of the coordinate system, relative to the yz plane. The only point they 
leave fixed is the origin of the coordinate system. Finally, parabolic 
transformations correspond to parabolic rotations about the origin of 
the coordinate system that leave only the w axis pointwise fixed. 



10. Some geometric properties of the isometric rotations 

In this section we aim to understand how the previously defined 
isometries act in S 3 and H 3 . In fact we will be interested only in 
the transformations represented by the matrices A and B, defined in 



(52) and (53), respectively. As we will see in Subsection 13.6, the re- 



negative parabolic rotations represented by the matrix C, defined in 
(54), generate no relative equilibria in H 3 , so we don't need to worry 



about their geometric properties for the purposes of this paper. 
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Since our earlier work focused on the curved n-body problem in the 
2-dimensional manifolds S 2 . and EI 2 ., we would like to see whether the 
above rotations preserve 2-dimensional spheres in §^ and and 2- 
dimensional hyperboloids in Hj^. We begin with the spheres. 

10.1. Invariance of 2-dimensional spheres in and ET^. Let us 
start with the K-positive elliptic-elliptic rotations in and consider 
first great spheres, which can be obtained, for instance, by the intersec- 
tion of S>^ with the hyperplane z — 0. We thus obtain the 2-dimensional 
great sphere 

(55) S% z = {(w, x, y, 0)| w 2 + x 2 + y 2 = kT 1 }, 

already defined in pi). Let (w, x, y, 0) to be a point on S 2 2 . Then the 



/ 


W L 


\ 


( 


cos 9 


— sin 9 








\ 


I 


w 


\ 




Xl 






sin 9 


cos 9 












X 






yi 












COS (f> 


— sine/ 






y 




\ 




) 


\ 








sin <j) 


cos</> 


) 


\ 





/ 



K-positive elliptic-elliptic transformation (52) takes (w,x, y,0) to the 
point (w 1 ,x 1 ,y 1 ,z 1 ) given by 



(56) 



/ w cos 9 — x sin 9 \ 
w sin 9 + x cos 9 
y cos 
V 1/ sin J 
Since, in general, y is not zero, it follows that z\ = ysm<ft = only 
if = 0, a case that corresponds to K-positive elliptic transforma- 
tions. In case of a /t-positive elliptic-elliptic transformation, the point 
(wi,X\, yi,Zi) does not lie on S 2 z because this point is not in the hyper- 
plane z — 0. Without loss of generality, we can always find a coordinate 
system in which the considered sphere is S 2 2 . We can therefore draw 
the following conclusion. 

Remark 2. For every great sphere of S^, there is a suitable system 
of coordinates such that ^-positive elliptic rotations leave the great 
sphere invariant. However, there is no system of coordinates for which 
we can find a K-positive elliptic-elliptic rotation that leaves a great 
sphere invariant. 

Let us now see what happens with non-great spheres of Such 
spheres can be obtained, for instance, by intersecting with a hy- 
perplane z = zo, where \zq\ < kT v I 2 and z ^ 0. These 2-dimensional 
sphere are of the form 



(57) 



S 2 



{(w, x, y, z) \ w 2 + x 2 



y 



Zq, z 
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where kq = — z 2 ) -1 / 2 is its curvature. 

Let (w, x, y, Zq) be a point on a non-great sphere S 2 0j ., o , given by 
some z as above. Then the ^-positive elliptic-elliptic transformation 
(52) takes the point (w, x, y, z ) to the point (w 2 , x 2 , y 2 , z 2 ) given by 



(5* 



/ 


w 2 


\ 




f 


cos 6 


— sin# 










x 2 








sin 9 


cos# 










2/2 














COS0 


— sin 


\ 


z 2 


) 




\ 








sin 


COSq 



\ f w 

X 



J 



\ 



y 

\z j 



( w cos 9 — x sin 9 \ 
w sin 9 + x cos 9 
y cos 4> — z sin 
y y sin + z cos y 

Since, in general, y is not zero, it follows that z 2 = y sin (fi+z cos = z 
only if = 0, that corresponds to K-positive elliptic transfor- 

mations. In the case of a K-positive elliptic-elliptic transformation, the 
point {w 2 , x 2 , y 2 , z 2 ) does not lie on S 2 KQ ZQ because this point is not in 
the hyperplane z = zq. Without loss of generality, we can always re- 
duce the question we posed above to the sphere S^ 0j2[) . We can therefore 
draw the following conclusion, which resembles Remark [2} 



Remark 3. For every non-great sphere of S^, there is a suitable system 
of coordinates such that K-positive elliptic rotations leave that non- 
great sphere invariant. However, there is no system of coordinates for 
which there exists a ^-positive elliptic-elliptic rotation that leaves a 
non-great sphere invariant. 

Since in we have z > 0, 2-dimensional spheres cannot be centred 
around the origin of the coordinate system. We therefore look for 2- 
dimensional spheres centred on the z axis, with z > |k| -1 / 2 (because 
z = |k| -1//2 is the smallest possible z coordinate in BP, attained only by 



the point (0, 0, 0, \k\ 
with a plane z = 



-1/2) 



Such spheres can be obtained by intersecting 



Zq, where z > \k\ 



-1/2 



and they are given by 



(59) 



i2,h 
'k ,z 



{(w, x,y, z)\ w 2 + x 2 + y 2 



Zn + K 



where h indicates that the spheres lie in a 3-dimensional hyperbolic 



space, and kq 



-1X-1/2 



> is the curvature of the sphere. 



Let (w, x, y, zq) be a point on the sphere S 2 .'/ 1 . Then the k- negative 



elliptic transformation B, given by (53) with s = 0, takes the point 
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(w, x, y, z ) to the point (w 3 , x 3 , y 3 , z 3 ) 



(60) 



/ U>3 \ 

X 3 

\Z 3 ) 



( cos 9 
sin# 


V 



— sin 9 
cos 9 







given by 
\ 




1 




/ w \ 

X 

y 

\ zo ) 



\ 

which also lies on the sphere S 



/ w cos 9 — x sin 9 \ 
w sin 9 + x cos 9 

y 

z 



2,h 

K, ,ZQ- 



) 

Indeed, since z 3 = z and w 3 + 
x l + V% = z o + ^ means that (w 3 , x 3 ,y 3 ,z 3 ) also lies on the sphere 

C2,h 

Since for any 2-dimensional sphere of we can find a coordinate 
system and suitable values for kq and zq such that the sphere has the 
form S^' o ft 2o , we can draw the following conclusion. 

Remark 4. For every 2-dimensional sphere of H^, there is a system 
of coordinates such that K-negative elliptic rotations leave the sphere 
invariant. 

Let us further see what happens with K-negative hyperbolic transfor- 
mations in H^. Let (w,x,y,z ) be a point on the sphere S^\, o . Then 
the K-negative hyperbolic transformation B, given by (53) with 9 = 0, 
takes the point (w, x, y, z ) to (1V4, x±, y^ z±) given by 



(61) 



/ lt>4 \ 

X,l 

1/4 
\Z 4 J 



( 1 




V° 
/ 



00 \ 

1 

cosh s sinh s 

sinh s cosh s J 



( w \ 

X 

y 

\z ) 



XV 

X 



\ 



y cosh s + zq sinh s 
V y sinh s + z cosh s J 

which does not lie on S^ ZQ . Indeed, since z 3 = y sinh s + z cosh s = z 
only for s = 0, a case we exclude because the above transformation is 
the identity, the point (w 4 , x 4 , y 4 , 24) does not lie on a sphere of radius 
\J Zq + Therefore we can draw the following conclusion. 

Remark 5. Given a 2-dimensional sphere of curvature kq = (zq + 



K 



-1/2 



with 2 > \k\ 



-1/2 



in H 3 K , 



there is no coordinate system 



for which some K-negative hyperbolic transformation would leave the 
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sphere invariant. Consequently the same is true about k- negative 
elliptic-hyperbolic transformations. 

We will further see how the problem of invariance relates to 2- 
dimensional hyperboloids in H^. Let us start with K-negative elliptic 
transformations in H^. 

10.2. Invariance of 2-dimensional hyperboloids in H^. Let us 

first check whether K-negative elliptic rotations preserve the great 2- 
dimensional hyperboloids of Hj^. For this consider the 2-dimensional 
hyperboloid 



(62) 



H 



K.y 



{(w, x, 0, z) | w 2 + x 2 — z 2 



K 



already defined in (19), and obtained by intersecting M. 3 K with the hy- 



perplane y — 0. Let (w, x, 0, z) be a point on H^ y . Then a K-negative 
elliptic rotation takes the point (w,x,0, z) to the point (w^, x$, y$, z 5 ) 
given by 

/ cos6 -sin6 \ 
sin^ cos^ 
10 
1 



(63) 



x 5 
\z 5 ) 



\ 

( w cos e — x sin 6 \ 
w sin e + x cos e 


z 



( w\ 



w 

X 



z 



\ z J 



which, obviously, also belongs to Since for any 2-dimensional 

hyperboloid of curvature k we can find a coordinate system such that 
the hyperboloid can be represented as H 2 Ky , we can draw the following 
conclusion. 

Remark 6. Given a 2-dimensional hyperboloid of curvature k in H^, 
there is a coordinate system for which the hyperboloid is invariant to 
K-negative elliptic rotations. 

Let us further check what happens in the case of K-negative hyper- 
bolic rotations. Consider the 2-dimensional hyperboloid of curvature k 
given by 



(64) 



H 



{(0,x,y, z) 



2 , 2 

x + y 



K 



}■ 



already defined in (17), and obtained by intersecting with the hyper- 
plane w = 0. Let (0, x, y, z) be a point on H 2 W - Then a /t-negative hy- 
perbolic rotation takes the point (0,x,y,z) to the point (u> 6 , x 6 , y 6 , z & ) 
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w 6 


\ 


( 1 











\ 


I ° ^ 











1 










X 














cosh s 


sinh s 




y 


\ 




) 


v° 





sinhs 


coshs 


/ 


\ z I 








/ 







\ 







given by 



(65) 



x 

y cosh s + z sinh s 
y y sinh s + z sinh s J 

which, obviously, also belongs to H^ w . Since for any 2-dimensional 
hyperboloid of curvature k we can find a coordinate system such that 
the hyperboloid can be represented as , we can draw the following 
conclusion. 

Remark 7. Given a 2-dimensional hyperboloid of curvature k in H^, 
there is a coordinate system for which the hyperboloid is invariant to 
re-negative hyperbolic rotations. 

Remark 8. The coordinate system in Remark [6] is different from the 
coordinate system in Remark [7], so re- negative elliptic- hyperbolic trans- 
formation don't leave 2-dimensional hyperboloids of curvature k invari- 
ant in H^. 

The next step is to see whether re-negative elliptic rotations preserve 



the 2-dimensional hyperboloids of curvature kq = —(y^ ■ 
of H^. For this consider the 2-dimensional hyperboloid 



K 



-1/2 



(66) 



H 



«0,S/0 



{(w,x,y,z) 



2 i 2 
W + X 



vl, y 



yo}, 

obtained by intersecting HI^ with the hyperplane y = y , with y ^ 0. 
Let (w, x, y , z) be a point on H^ y . Then a re-negative elliptic rotation 



takes the point {w } x, yo, z) to the point (w 7 , x 7 , y 7j zj) 



(67) 



/ w 7 \ 

x 7 

yi 

\z 7 / 



( cos6 -sin6 \ 

sin 9 cos 9 

10 

1 



w 

X 

yo 

z 



given by 
\ 



V z I 



( w cos 9 — x sin 9 \ 
w sin 9 + x cos 9 

yo 

z 



which, obviously, also belongs to H^ 0)J/Q . Since for any 2-dimensional 
hyperboloid of curvature re we can find a coordinate system such that 
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the hyperboloid can be represented as H 2 , we can draw the following 
conclusion. 

Remark 9. Given a 2-dimensional hyperboloid of curvature kq = 
~{Vo ~ ft -1 ) -1 / 2 7^ k in H^, there is a coordinate system for which 
the hyperboloid is invariant to K-negative elliptic rotations. 

Consider further the 2-dimensional hyperboloid of curvature k,q = 
— (wq — /t -1 ) -1 / 2 ^ k given by 

( 68 ) U l ,w = i( w , x , V, z ) \ x 2 + y 2 -z 2 = k- 1 -wl,w = wo}, 

obtained by intersecting with the hyperplane w = w ^ 0. Let the 
point (wo, x, y, z) lie on H 2 o ^ wo . Then a K-negative hyperbolic rotation 
takes the point (wq, x, y, z) to the point (w 8 , x$, y&, z$) given by 



(69) 



x$ 

ys 

\ z 8 J 



( 1 









1 
cosh s sinh s 
sinh s cosh s J 

Wo \ 

X 

y cosh s + z sinh s 
y y sinh s + z sinh s J 



\ 



( Wo 

X 

y 

\ z J 



which, obviously, also belongs to 



Since for any 2-dimensional 



hyperboloid of curvature k = —(wq — k^ 1 ^ 2 )^ 1 we can find a coordinate 
system such that the hyperboloid can be represented as H 2 m , we can 
draw the following conclusion. 

Remark 10. Given a 2-dimensional hyperboloid of curvature k,q = 
' m BE, there is a coordinate system for which the 



- (Wn 



K 



hyperboloid is invariant to ^-negative hyperbolic rotations. 

Remark 11. Since the coordinate system in Remark [9] is different from 

K-negative elliptic-hyperbolic 



10 



the coordinate system in Remark 
transformation don't leave 2-dimensional hyperboloids of curvature 
Kg 7^ k invariant in ET^. 

11. Relative equilibria 

The goal of this section is to introduce the concepts we will explore 
in the rest of the paper, namely the relative equilibrium solutions (also 
called relative equilibrium orbits or, simply, relative equilibria) of the 
curved n-body problem. For relative equilibria, the particle system 
behaves like a rigid body, i.e. all the mutual distances between the 
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point masses remain constant during the motion. In other words, the 
bodies move under the action of an element belonging to a rotation 
group, so, in the light of Section [9j we can define 6 types of relative 
equilibria in M^: 2 in §^ and 4 in H^. In each case, we will bring the 
expressions involved in these natural definitions to simpler forms. We 
will later see that 1 of the 4 types of relative equilibria we define in 
does not translate into solutions of the equations of motion. In what 
follows, the upper T will denote the transpose of a vector. 



11.1. Definition of /t-positive elliptic relative equilibria. The 

first kind of relative equilibria we will introduce in this paper are in- 
spired by the K-positive elliptic rotations of 

Definition 6. (^-positive elliptic relative equilibria) Let q° = 

(Qij Q25 • • • 5 Qn) ^ e a nonsingular initial position of the point particles 



of masses mi, m 2 
K > 0, where q° = 
the form q 



,0 



> 0,71 > 
,0 ^0 



2, on the manifold S K , i.e. for 
Then a solution of 



, n. 



(A[<£) T ,A[<£] T , 



(70) 



A(t) 



( cos at 
sin at 


V 



i = 1,2, 
*4jq°] T ) of system (44), with 

0\ 




1 



— sin at 
cos at 








1 





/ 



where a 7^ denotes the frequency, is called a (simply rotating) k- 
positive elliptic relative equilibrium. 

Remark 12. In A, the elements involving trigonometric functions 
could well be in the lower right corner instead of the upper left corner 
of the matrix, but the behaviour of the bodies would be similar, so we 
will always use the above form of the matrix A. 



Ifr< =: VK ) 2 + (x?) 



such that u> 



Tj COS di, 



^ 2 , we can find constants a, 
x? = r,; sinaj, i — 1, 2, 



1,2, 



n. Then 



A(t)[q, 



/ w® cos at — x° sin at \ 
sin at + x° cos at 



\ 



( Ti cos di cos at 



Ti sin dj sin 



at \ 



Ti cos a» sin at + n sin ck cos at 



3? 



/ 

^ Ti cos(at + aj) \ 
Ti sin (at + d. 



\ 



3° 



/ 



1,2, 
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11.2. Definition of K-positive elliptic-elliptic relative equilibria. 

The second kind of relative equilibria we introduce are inspired by the 
^-positive elliptic-elliptic rotations of S^. 

Definition 7. (/{-positive elliptic-elliptic relative equilibria) Let 



q° = (q'j 1 , q°, . . . , q°) be a nonsingular initial position of the bodies of 
masses mi, m 2 , . . . , m n > 0, n > 2, on the manifold S^, i.e. for n > 0, 
, , x° , y®, z°), i = 1, 2, . . . , n. T/ien a solution of the form 
<" i r « r -° lT £[q°] T ) of system (0), 



w/iere q° = (//•'/../• //'/. :; 



q=(B[q?] T ,B[q§] 



(71) 



B(t) 



( cos at 
sin at 





— sin at 
cos at 










cos (3t 
sin fit 



\ 






- sin (3t 
cos (3t J 

where a, (3 ^ are the frequencies, is called a (doubly rotating) k- 
positive elliptic- elliptic relative equilibrium. 



- 

el, % = 1,2 



0N2 



Pi 



a/ (yf) 2 + (<2°) 2 , we can find constants 



Pi cos6j, and z\ 



B(t)H 



OlT 



n, such that u>° = TjCosaj,a;° = rjSinaj,|/° 
Pi sin bi, i — 1, 2, . . . , n. Then 

/ cos at — x° sin at \ 
sin at + x° cos at 
y® cos /3t — 2° sin /3t 
\ sin /3t + 2° cos (3t J 

/', so 1 (/ ; sin at \ / /•, (•()>(/>/ • (•/, ) \ 



1,2 



/ Ti cos Oi cos at 

r j cos a>{ sin at + r j sin cos at 
Pi cos 6j cos /3t — pj sin 6, sin (3t 
y Pi cos bi sin /3t + pj sin 6, cos /3t y 

,71. 



rj sin(at + Oj) 
Pi cos(/3t + 6j) 
^ pi sin(/3t + 6j) y 



11.3. Definition of K-negative elliptic relative equilibria. The 

third kind of relative equilibria we introduce here are inspired by the 
k- negative elliptic rotations of H^. 

Definition 8. (^-negative elliptic relative equilibria) Let q° = 



,0 • 



, q^J be a nonsingular initial position of the point particles 
of masses mi, m,2, 

(C[ q ?p,e[ q ^,.... 



. ,m n > 0,n > 2, in M^, i.e. for k < 0, where 
i = 1,2, ... ,n. Then a solution of the form q = 

01T\ 



(72) 



C(t) 



C[q"] J ) of system (|44j), with 

( cos at — sin at \ 

sin at cos at 

10 

1 



V 
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where a ^ is the frequency, is called a (simply rotating) n-negative 
elliptic relative equilibrium. 



If r, = 
that w? 



(73) 



w 



0\2 



+ (X 



0\2 



we can find a,- e 



1,2, 



n, such 



Tj COS di, X 







C(*)[q? 



OlT 



n sin Oj, z = 1, 2, . . . , n, so 

/ w® cos at — a;° sin at \ 
sin at + cos at 



V 



/ Ti cos a,- cos at 



Ta sin a,- sin 



at \ 



V 

1,2, 



Tj, cos a, sin at + r j sin aj cos at 

71. 



/ 

^ rj cos(at + Oj) 
rj sin (at + a,) 

v 



/ 



11.4. Definition of /{-negative hyperbolic relative equilibria. 

The fourth kind of relative equilibria we introduce here are inspired 
by the /c-negative hyperbolic rotations of H^. 

Definition 9. (/{-negative hyperbolic relative equilibria) Letq = 

q°) be a nonsingular initial position of the bodies of masses 



mi, 7712, • 



771„ > 0, 71 > 2, 



positions are q^ 



o 



system (44) o/ t/ie /orm q = (D[q ( j ) ] T , P[q2] T , . 

/ 1 
10 



i.e. for k < 0, where the initial 
i = 1,2, ... ,n. Then a solution of 

with 



(74) 



V(t) 



\ 



IV) 



cosh/3t sinh/3t 
V sinh/3t cosh/3t I 

where (3 ^ denotes the frequency, is called a (simply rotating) k- 
negative hyperbolic relative equilibrium. 

If rji := a/ (z®) 2 — (Vi) 2 , we can find constants h L G K, i = 1, 2, . . . , n, 
such that y° = r\i sinh 6j and = r\i cosh fej, i = 1, 2, . . . , n. Then 



/ 



OlT 



7C 



X? 



\ 



/ 



70 



?/° cosh bt + z? sinh 6t 
V 7/f sinh 6t + z 4 ° cosh fet y 

\ / 



rji sinh 6j cosh bt + 77j cosh 6j sinh bt 
y 77j sinh bi sinh 6t + 77j cosh bi cosh 6t y 



,-,.u 

?7i sinh(6t + bi 
y 77, cosh(6t + bi 
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,77. 



11.5. Definition of K-negative elliptic-hyperbolic relative equi- 
libria. The fifth kind of relative equilibria we introduce here are in- 
spired the k- negative elliptic-hyperbolic rotations of H^. 



Definition 10. (/{-negative elliptic-hyperbolic relative equilib- 
ria) Let q° = (q°, q°, . . . , q°) be a nonsingular initial position of the 
point particles of masses mi, 777.2, ■ ■ ■ , Ta n > 0,n > 2, in i.e. for 
k < 0, where q° = (tu°, yf, z°), i = 1,2, ... , n. Then a solution of 

(£[q?] T ,£[ q2 T, 



system (44) of the form q 



0lT\ 



with 



(75) 



£(t) 



/ cos at 
sin at 


V 



— sin at 
cos at 











cosh fit sinh fit 

sinh fit cosh fit J 



where a,fi 7^ denote the frequencies, is called a (doubly rotating) 
K-negative elliptic-hyperbolic relative equilibrium. 



w 



0\2 



+ 



,0\2 



if n 

aj, 6j G R, i = 1 
7/j sinh bi, and z° = r]i cosh 6 



(.;•',')-, r/j := \/Jzfy~^ljjf)^ , we can find constants 
such that = rjCosaj,a;° = rjSinaj,|/° = 



n. 



1,2, 



£(t)[q: 



/ 



n. 



Then 
sin at \ 



cos at — x£ 
sin at + x° cos at 

/? 



7/" cosh fit + sinh /3t 
y y° sinh /3t + z° cosh /3t y 



/ 



Ti cos a, 



cos at — r j sin aj sin at \ 



r j cos ai sin at + 7*j sin aj cos at 
r/j sinh /3, cosh /3t + 77, cosh $ sinh fit 
\ rji sinh /3, sinh fit + r\i cosh fii cosh /3t / 



1,2, 



,77. 



/ Ti cos (at + ai) \ 
Ti sin (at + ai) 
rji sinh(/3t + bi) 
\ r]i cosh( fit + h) J 



11.6. Definition of /{-negative parabolic relative equilibria. The 

sixth class of relative equilibria we introduce here are inspired by the 
K-negative parabolic rotations of H^. 



Definition 11. (/t-negative parabolic relative equilibria) Con- 
sider a nonsingular initial position q° = (q°, q°, • • • , q°) 0/ t/ie pomt 
particles of masses m^mj, . . . , m n > 0,n > 2, on t/ie manifold M.^, i.e. 
for k < 0, where q° 



i, 

7° 



hi Z i)i 



1,2, 



,77. 



JTien a solution 
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(76) 



of the form q = (J r [q ( l] T , J 7 ^] 7 , . . . , J 7 ^] 7 ") of system (44), with 

( 1 \ 

1 -t t 
t 1 — 1 2 /2 t 2 /2 

\ o t -t 2 /2 i + 1 2 /2 y 

is called a (simply rotating) n-negative parabolic relative equilibrium. 

For simplicity, we denote := w®,Pi := x°,7i := 5, := z®, i = 
1,2, ... ,n. Then parabolic relative equilibria take the form 



/ 



( 



\ 



X 



n. 



x *t + y»(l-t 2 /2) + z»t 2 /2 
\x°t-y°t 2 /2 + z°(l + t 2 /2) ; 

oh \ 
Pi + (Si - n)t 
7i + Pit + & - 7i )t 2 /2 
\6i + Pit + (Si - n)t 2 /2, ) 

11.7. Formal expressions of the relative equilibria. To summa- 
rize the previous findings, we can represent the above 6 types of relative 
equilibria of the 3-dimensional curved n-body problem in the form 



1,2, 



q = (qi, q2, ■ ■ ■ , q n ), q» = (w i: x h y h % = 1, 2, 



n, 



(77) 



[k > 0, elliptic] 



w-i(t) = Vi cos(at + Oi) 
Xi(t) = Ti sin(ai + tij) 
Ui(t) = Hi (constant) 
Zi(t) = Zi (constant), 

with w 2 + x 2 = r 2 , r 2 + y 2 + z\ = z = 1, 2, . . . , n; 

' ' Wi(t) = Ti cos(at + ai) 
x.i(t) = Ti sin(at + dj) 
Vi(t) = PiCos(pt + hi) 
K Zi(t) = pism(pt + bi), 

with w 2 + x 2 = r 2 , yf + z 2 = pj, so rf + p\ = k~ x , i = l,2,..., n; 

Wi(t) = Ti cos(at + Oj) 



(7 



[« > 0, elliptic— elliptic] 



(79) 



[k < 0, elliptic] : < 



Xi(t) = Tj sin(at + a,) 
= Vi (constant) 
K Zi(t) = Zi (constant), 
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with W~ + x\ 



rt + Vi 



1,2,..., n; 



with 



[k < 0, hyperbolic] : 



(81) [ac < 0, elliptic— hyperbolic] 



Wi(t) = Wi (constant) 
Xi(t) = Xi (constant) 
yi(t) = riismh((3t + b { ) 
^(t) = rji cosh(/3t + bi) , 

Ki . 1 — 1 . 2, . . . . . Th'. 

u>i(t) = r, cos(at + a$) 
Xi(t) = n sin(at + a,) 
Vi(t) = r/isinh(/3t + bi) 
ji(t) = ?7iCosh(/3t + bi), 



with 



wf + X? 



Vi-h= ~Vi .son- Wi 



1,2, 



, n: 



(82) [k < 0,parafeo/ic] : 



' Wi(t) = ai (constant) 

Xi{t) = (3i + (5i - 7i)t 

Vi(t) = ji + Pit + (5i-ji)t 2 /2 
[ Zi (t) =6i + f3 i t + (6 i - ll )t 2 /2, 



with aj + ft + t? - <S f 2 



1,2, 



,n. 



12. Fixed Points 



In this section we introduce the concept of fixed-point solution of 
the equations of motion, show that fixed points exist in 8^, provide a 
couple of examples, and finally prove that they cannot show up in 
and in hemispheres of 8^. 

12.1. Existence of fixed points in 8^. Although the goal of this 
paper is to study relative equilibria of the curved n-body problem in 
3-dimensional space, some of these orbits can be generated from fixed- 
point configurations by imposing on the initial positions of the bodies 
suitable nonzero initial velocities. It is therefore necessary to discuss 
fixed points as well. We start with their definition. 



Definition 12. A solution of system (36) is called a fixed point if it is 
a zero of the vector field, i.e. Pi(t) = V qi i7 K (q(£)) = for all f 6 1, 
i = 1, 2, ... ,n. 



In [21], [2S], and [2Z], we showed that fixed points exist in SJ 2 ., but 
they don't exist in H 2 .. Examples of fixed points are the equilateral 
triangle of equal masses lying on any great circle of S 2 . in the curved 3- 
body problem and the regular tetrahedron of equal masses inscribed in 
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§ 2 in the 4-body case. There are also examples of fixed points of non- 
equal masses. We showed in [22] that, for any acute triangle inscribed 
in a great circle of the sphere § 2 , there exist masses mi, 777.2,777.3 > 
that can be placed at the vertices of the triangle such that they form a 
fixed point, and therefore can generate relative equilibria in the curved 
3-body problem. All these examples can be transferred to S 3 .. 

12.2. Two examples of fixed points specific to S 3 .. We can con- 
struct fixed points of § 3 for which none of its great spheres contains 
them. A first simple example occurs in the 6-body problem if we take 
6 bodies of equal positive masses, place 3 of them, with zero initial 
velocities, at the vertices of an equilateral triangle inscribed in a great 
circle of a great sphere, and place the other 3 bodies, with zero initial 
velocities as well, at the vertices of an equilateral triangle inscribed 
in a complementary great circle (see Definition |3| of another great 
sphere. Some straightforward computations show that 6 bodies of 
masses 777.1 — m 2 — m 3 — m i — m 5 — m e =: m > 0, with zero 
initial velocities and initial conditions given, for instance, by 

0, zi = 0, 

0, z 2 = 0, 



Wi 


= «r 1/2 , 


Xi 


= 0, 


yi = 


w 2 










2 ' 


2 ' 


w 3 




X3 


V3K- 1 / 2 


ys = 


2 ' 


2 ' 


W4 


= 0, 




= 0, 


2/4 = 


w 5 


= 0, 


x 5 


= 0, 


2/5 = 


w 6 


= 0, 


x 6 


= 0, 


Ve = 



0, z 3 = 0, 

fiT 1/2 , z 4 = 0, 
re" 1 / 2 v/3re- 1/2 

z% 



2 ' 2 ' 

re- 1 / 2 v^re" 1 / 2 



form a fixed point. 

The second example is inspired from the theory of regular polytopes, 
[Hj, [TS]. The simplest regular polytope in IR 4 is the pentatope (also 
called 5-cell, 4-simplex, pentachrone, pentahedroid, or hyperpiramid) . 
The pentatope has Schlafli symbol {3, 3, 3}, which translates into: 3 
regular polyhedra that have 3 regular polygons of 3 edges at every 
vertex (i.e. 3 regular tetrahedra) are attached to each of the pentatope's 
edges. (From the left to the right, the numbers in the Schlafli symbol 
are in the order we described them.) 

A different way to understand the pentatope is to think of it as the 
generalization to IR 4 of the equilateral triangle of IR 2 or of the regular 
tetrahedron of IR 3 . Then the pentatope can be constructed by adding 
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to the regular tetrahedron a fifth vertex in M 4 that connects the other 
four vertices with edges of the same length as those of the tetrahedron. 
Consequently the pentatope can be inscribed in a sphere § 3 , in which it 
has no antipodal vertices, so there is no danger of encoutering singular 
configurations for the fixed point we want to construct. Specifically, 
the coordinates of the 5 vertices of a pentatope inscribed in the sphere 
S 3 . can be taken, for example, as 

w 1 = /sT 1/2 , xi = 0, yi= 0, zx = 0, 



W2 = — , x 2 = , 2/2 =0, z 2 = 0, 

«* = -— -3 = -^-, ^3 = ^^, ,3 = 0, 

V5K- 1 / 2 ^- 1/2 V^K- 1 / 2 

— X4 -~^7T' VA ~~^w ZA -^vr> 

VSk- 1 / 2 VSk- 1 / 2 VZk- 1 / 2 

Straightforward computations show that the distance from the origin 
of the coordinate system to each of the 5 vertices is k -1 / 2 and that, for 
equal masses mi = m,2 = = = =: m > 0, this configuration 



K 


-1/2 




4 


K 


-1/2 




4 


K 


-1/2 




4 


K 


-1/2 


4 



produces a fixed-point solution of system (44) for k > 0. Like in the 
previous example, this fixed point lying at the vertices of the pentatope 
is specific to S 3 . in the sense that there is no 2-dimensional sphere that 
contains it. 

It is natural to ask whether other convex regular polytopes of 1R 4 
can form fixed points in S 3 if we place equal masses at their vertices. 
Apart from the pentatope, there are five other such geometrical ob- 
jects: the tesseract (also called 8-cell, hypercube, or 4-cube, with 16 
vertices), the orthoplex (also called 16-cell or hyperoctahedron, with 8 
vertices), the octaplex (also called 24-cell or polyoctahedron, with 24 
vertices), the dodecaplex (also called 120-cell, hyperdodecahedron, or 
polydodecahedron, with 600 vertices), and the tetraplex (also called 
600-cell, hypericosahedron, or polytetrahedron, with 120 vertices). All 
these polytopes, however, are centrally symmetric, so they have an- 
tipodal vertices. Therefore, if we place bodies of equal masses at their 
vertices, we encounter singularities. Consequently the only convex reg- 
ular polytope of 1R 4 that can form a fixed point if we place equal masses 
at its vertices is the pentatope. 

12.3. Nonexistence of fixed points in H 3 and hemispheres of 
S 3 .. We will show further that there are no fixed points in H 3 . or in 
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any hemisphere of S 3 K . In the latter case, for fixed points not to exist 
it is necessary that at least one body is not on the boundary of the 
hemisphere. 

Proposition 3. (No fixed points in H^) For k < 0, there are no 
masses mi, m 2 , . . . , m n > 0, n > 2, that can form a fixed point. 

Proof. Consider n bodies of masses m 1; m 2 , . . . , m n > 0, n > 2, lying on 
ET^, in a nonsingular configuration (i.e. without collisions) and with zero 
initial velocities. Then one or more bodies, say, mi,m 2 , . . . ,m^ with 
k < n, have the largest z coordinate. Consequently each of the bodies 
mk+i, • • • , m n will attract each of the bodies mi, mj, . . . , m k along a ge- 
odesic hyperbola towards lowering the z coordinate of the latter. For 
any 2 bodies with the same largest z coordinate, the segment of hyper- 
bola connecting them has points with lower z coordinates. Therefore 
these 2 bodies attract each other towards lowering their z coordinates 
as well. So each of the bodies mi, m 2 , . . . , m k will move towards lower- 
ing their z coordinate, therefore the initial configuration of the bodies 
is not fixed. □ 

Proposition 4. (No fixed points in hemispheres of S^) For k > 0, 

there are no masses mi, m 2 , ■ ■ ■ , m n > 0, n > 2, that can form a fixed 
point in any closed hemisphere o/§^ (i.e. a hemisphere that contains 
its boundary), as long as at least one body doesn't lie on the boundary. 

Proof. The idea of the proof is similar to the idea of the proof we gave 
for Proposition |3j Let us assume, without loss of generality, that the 
bodies are in the hemisphere z < and they form a nonsingular initial 
configuration (i.e. without collisions or antipodal positions), with at 
least one of the bodies not on the boundary z — 0, and with zero initial 
velocities. Then one or more bodies, say, mi,m 2 , ■ ■ ■ , m^, with k < n (a 
strict inequality is essential to the proof), have the largest z coordinate, 
which can be at most 0. Consequently the bodies mk+i, ■ ■ ■ ,m n have 
lower z coordinates. Each of the bodies m^+i, ■ ■ ■ , m n attract each of 
the bodies mi, m 2 , . . . , m k along a geodesic arc of a great circle towards 
lowering the z coordinate of the latter. 

The attraction between any 2 bodies among mi, m 2 , . . . , m^ is either 
towards lowering each other's z coordinate, when z < or along the 
geodesic z = 0, when they are on that geodesic. In both cases, however, 
composing all the forces that act on each of the bodies mi, m 2 , . . . , m k 
will make them move towards a lower z coordinate, which means that 
the initial configuration is not fixed. This remark completes the proof. 

□ 
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Part 4. CRITERIA AND QUALITATIVE BEHAVIOUR 

13. Existence criteria for the relative equilibria 

In this section we establish criteria for the existence of K-positive 
elliptic and elliptic-elliptic as well as K-negative elliptic, hyperbolic, and 
elliptic-hyperbolic relative equilibria. These criteria will be employed 
in later sections to obtain concrete examples of such orbits. We close 
this section by showing that K-negative parabolic relative equilibria do 
not exist in the curved n-body problem. 

13.1. Criteria for /t-positive elliptic relative equilibria. We pro- 
vide now a criterion for the existence of simply rotating K-positive ellip- 
tic orbits and then prove a corollary that shows under what conditions 
such solutions can be generated from fixed-point configurations. 

Criterion 1. (K-positive elliptic relative equilibria) Consider the 
point particles of masses m\,m<i, . . . ,m n > 0, n > 2, in i.e. for 
k > 0. Then system (44) admits a solution of the form 

q = (qi, q2, • • • , q n ), q* = (w h x h y h % = 1, 2, . . . , n, 

Wi(t) = Ti cos(at + ai),Xi(t) = sin(at + Oi),yi(t) = y h z^t) = z h 

with wf + xf = rf, rf + yf + zf = k~ x , and Zi constant, i — 1,2, . . . ,n, 
i.e. a (simply rotating) K-positive elliptic relative equilibrium, if and 
only if there are constants r iy a^y^ Zi, i — 1, 2, . . . , n, and such 
that the following 4n conditions are satisfied: 

/oo\ <^ nijK^irj cos aj - KVijVi cos ai) 2 2 

(83) ^ -± (1 - j^.W 2 = ^ - l ) ar * cos <*> 

(84) ^ (1 - k*vI)W = ^ ~ )a n Sm ah 



3=1 



1 - « 2 i/? f ) 3 / 2 

3=1 ^ U ' 



lti\ mjK Zl2 (Zj KUjjZj) _ 2 2 



1 K'l-'fj ) ''' 



i = l,2,...,n, where = r^j cos(<2j — a,j) + y^j + ZjZj, i,j 
1,2..., n, i^j. 
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Proof. Consider a candidate of a solution q as above for system (44). 
Some straightforward computations show that 

Vij ■= <ii <ij = TiTj cos(aj - a,j) + ViVj + ZiZj, i, j = 1, 2 . . . , n, i ^ j, 

qi © qi = a 2 r 2 , i = 1, 2, . . . , n, 

Wi = —a 2 ri cos(at + Oj), Xj = — a 2 ?^ sin(at + Oj), 

Sl< = Zi = 0, i = 1,2, . . . ,n. 
Substituting the suggested solution and the above expressions into sys- 



tem (44), for the w coordinates we obtain conditions involving cos(crf + 
<2j), whereas for the x coordinates we obtain conditions involving sin(at+ 
<2j). In the former case, using the fact that cos(at + aj) = cos at cosa^ — 
sin at sin a^, we can split each equation in two, one involving cos at and 
the other sin at as factors. The same thing happens in the latter case 
if we use the formula sin(at + a«) = sin at cos + cos at sin a.j. Each 



of these equations are satisfied if and only if conditions (83) and (84) 



take place. Conditions (85) and (86) follow directly from the equa- 



tions involving the coordinates y and z. This remark completes the 
proof. □ 

Criterion 2. (/{-positive elliptic relative equilibria generated 
from fixed-point configurations) Consider the point particles of 
masses mi,m 2 , . . . ,m n > 0, n > 2, in S>^ ; i.e. for k > 0. Then, for any 



a^O, system (|44j) admits a solution of the form (77): 
q = (qi, Q2, • • • , Qn), = (w h x h y h z t ), i = 1, 2, 

Wi(t) = Ti cos(at + ai), Xi(t) = Ti sin(at + a t ), y^t) = y h z^t) = z h 

with wf + xf = rf, rf + yf + zf = k~ 1 , and y^ Zi constant, i — 1, 2, . . . , n, 
generated from a fixed point, i.e. a (simply rotating) K-positive elliptic 
relative equilibrium generated from the same initial positions that would 
form a fixed point for zero initial velocities, if and only if there are con- 
stants Ti, ai, yi, Zi, i — 1, 2, . . . , n, such that the following 4n conditions 
are satisfied: 

rrijK 3 / 2 {tj cos a, - KU^n cos a*) 

(87) ^ (1 - )3/2 - °' 

^ m i ^ 3 / 2 (r i sin aj - KUjf\ sin a,) _ 

(88) ^ (1 - k*4) 3 / 2 " ' 



Relative equilibria in the 3-dimensional curved n-body problem 



55 



19) 



E 

3=1 



rrijK 



kV 2 ) 3 / 2 



(90) 



ErrijK 
— 

3 = 1 



3/2 (Zj - KUijZi) 



^2)3/2 



1,2, 



,n, where Vi. 



y - nrjcoslai - aj) + yii/j + Z{Zj, i,j = 
, 2, . . . , n, i 7^ j, and one of the following two properties takes place: 

(i) n = fiT 1/2 for all i G {1,2 ... , n}, 

(ii) there is a proper subset Ic {1,2, 



i G X and r,- 



-1/2 



sttc/i t/iat Ti = /or a// 



/or a// j G {l,2,...,n}\X. 



Proof. We are seeking a simply rotating elliptic relative equilibrium, 
as in Criterion [TJ that is valid for any ct^O. But the solution is also 
generated from a fixed-point configuration, so the left hand sides of 



equations (83), (84), (85), and ( 



conditions (87), (88), (89), and (90). However, the right hand sides of 



equations (83), (84), (85), and (86) necessarily vanish, thus leading to 



86) must also vanish, so we have the 



An conditions: 



l)a 2 rj cos aj 
l)o; 2 r,; sin 
0, 



0, i 



0, % 



1,2, 



1,2, 
1,2,. 



, n, 



na 2 r 2 Zi = 



, i = 1,2,. . . ,n. 



Since cr ^ and there is no 7 G R such that the quantities sin 7 and 
cos 7 vanish simultaneously, the above 4n conditions are satisfied in 
each of the following cases: 

for all 



(a) ri = (consequently w t 
% G {1,2,. ..,n}, 

(b) rj = k -1 / 2 (consequently u> 2 + x 2 
% G {1,2. ..,n}, 

(c) there is a proper subset 1c {1,2, ... ,n} such that r. 



and y 2 + 2 2 
= k _1 and ?/j = 



0) for all 







con- 



sequently yf + zf = k~ 1 ) for all i G X and r,- = re -1 / 2 (consequently 

Vj = z i = 0) f° r a ^ i e {!> 2, . . . , n} \ X. 

In case (a), we recover the fixed point, so there is no rotation of 
any kind, therefore this case does not lead to any simply rotating re- 
positive elliptic relative equilibrium. As we will see in Theorem [2j case 
(b), which corresponds to (i) in the above statement, and case (c), 
which corresponds to (ii), lead to relative equilibria of this kind. This 
remark completes the proof. □ 
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13.2. Criteria for /t-positive elliptic-elliptic relative equilibria. 

We can now provide a criterion for the existence of doubly rotating 
K-positive elliptic-elliptic relative equilibria and a criterion about how 
such orbits can be obtained from fixed-point configurations. 

Criterion 3. (^-positive elliptic-elliptic relative equilibria) Con- 
sider the bodies of masses mi,m 2 , . . . ,m n > 0,n > 2, in S 3 ., i.e. for 
k>0. Then system (44) admits a solution of the form (78): 

q = (qi, q2, • • • , q*), q; = (w u x h y h z t ), % = 1, 2, . . . , n, 

Wi(t) = Ti cos(at + ai), Xi(t) = Ti sm(at + a^), 

Vi{t) = Pi cos(/3t + bi), Zi(t) = pi sin(/3t + h), 

with w 2 + x\ = r 2 ,y 2 + z 2 p 2 ,r 2 + p 2 = kT 1 , i = 1,2,..., n, i.e. a 
(doubly rotating) K-positive elliptic- elliptic relative equilibrium, if and 
only if there are constants r^, pi,ai,bi, i = 1,2, ... ,n, and a,/? ^ 0, 
such that the following An conditions are satisfied 

(91 i 

Em,K 3 ' 2 (r,- cosa ? - — fix^r; cos a^) . 2 2 ^,22 2\ 
n _ W 2 = ( K£ * r i +K P Pi~ a ) r i COSa i, 

j=i 

(92) 

n 

EllljK,' IT,- bill Uj — IXUJiil i bill Uj I ,99 ^9 9 9n 

^ ^_ j,,, 2W 2 %1 = (™ 2 r 2 + Kfp\ - a 2 ) ri sma, 

3=1 
(93) 

n 

Eniji% ' ijy,cusw,- — i%uJiAUi<jjisun . 9 9 ^99 ^9 N 
^ _ X,2 ^/2 = ( Ka T i + "Pfi ~ 

3=1 

(94) 

n 

^ ^ ^-^—-f— ^ = + ./3 2 p 2 - f) P , sin 6, 

i = 1,2, ... ,n, where = r^j cos(a,j — a,) + piPj cos(6j — bj), i,j 
1,2, ...,n, i^j, and rf + p 2 = K" 1 , % — l,2,...,n. 



mjK 3 / 2 (rj cos a,- — Ku^ri 


COS Ctj) 


(1 - k 2 u%) 3 I 2 




mjK 3 ^ 2 (rj smaj — noj^ri 


sin aj) 


(1 - K 2 U% )3/2 




irijK 3 ^ 2 (pj cos 6j — KUijPi 


cos 6j) 


(1 - k 2 u; 2 ) 3 / 2 




mjK 3 / 2 (pj sinbj — Ktu i jp l 


sin 6j) 



Proof. Consider a candidate q as above for a solution of system (44). 
Some straightforward computations show that 

Uij := qiOq, = nrj cos(a i -a i )+PiPj cos(b—bj), i, j = 1, 2 . . . , n, i ^ j, 
qi ^ = a 2 r 2 + (3 2 p 2 , i = l,2,...,n, 
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W; 



-a 2 ri cos(at + di 



X ; 



-a 2 ri sin(at + ai 



-/3 2 piCos(/3t + bi),Zi = -/3 2 pism(l3t + b { ), i 



1,2, 



n. 



Substituting q and the above expressions into system (44), for the w 
coordinates we obtain conditions involving cos(at + aj), whereas for the 
x coordinates we obtain conditions involving sin(at + aj). In the former 
case, using the fact that cos(at+cij) = cos at cos a^ — sin at sina^, we can 
split each equation in two, one involving cos at and the other sin at as 
factors. The same thing happens in the latter case if we use the formula 
sin(at + di) = sin at cos aj + cos at sin a^. Each of these equations are 



satisfied if and only if conditions (91) and (92) take place. 



For the y coordinate, the substitution of the above solution leads to 
conditions involving cos(/3 + bi), whereas for z coordinate it leads to 
conditions involving sin(/3 + bi). Then we proceed as we did for the w 



and x coordinates and obtain conditions (93) and (94). This remark 
completes the proof. □ 

Criterion 4. (/t-positive elliptic-elliptic relative equilibria gen- 
erated from fixed-point configurations) Consider the point parti- 
cles of masses mi,m 2 ,. . . , m n > 0,n > 2, in §> 3 , i.e. for n > 0. Then, 



for any a, (3 ^ ; system (44) admits a solution of the form (|78j): 

q = (qi, q2, • • • , q»), q* = (m, x h y i} z { ), 



1,2, 



Wj(t) = Tj cos(at + a,i) 
Vi(t) = PiCos((3t + bi) 



with wf + x • 



2^2 



rt,Vt + ztP 



Xi(t) = Ti sin(at + aj), 
Zi(t) = pisin(f3t + bi), 

= G {1,2, ... ,n}, gener- 



ated from a fixed-point configuration, i.e. a (doubly rotating) K-positive 
elliptic- elliptic relative equilibrium generated from the same initial po- 
sitions that would form a fixed point for zero initial velocities, if and 



only if there are constants 7*j, p^, a*, bi 
relationships below are satisfied: 

n 

05) E 



1TljK 3 / 2 (rj COS Qj 



1,2, 



KUJijTi cos a« 



3=1 



1 - K 2 c4)3/2 



such that the 4n 



(96) 



£ 

3=1 
3& 



3/2, 



Tj sm a,j — KUijri sm a,{ 



K 2 w 2)3/2 



(97) 



E 

3=1 

3& 



mjK 



3 ^ 2 (pj cos bj — KUijPi cos bi 



^4)3/2 
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Y> rrijK^ipj sin - ftxu^pj sin bj) _ 

1 j 2^ (1-^2)3/2 

z = 1, 2, . . . , n, where cUy = r^j cos(aj — ay) + pip./ cos(6j — z,j = 
1,2, ... ,n, i 7^ j ; and rf + p 2 = R 2 = k -1 , i = 1,2, ... ,n, and, 
additionally, one of the following properties takes place: 

(i) there is a proper subset J C {1, 2, . . . , n} such that = /or a// 
i E J and pj = /or a/Z j G {1, 2, . . . , n} \ J , 

(ii) the frequencies a, (3 ^ satisfy the condition \a\ = \/3\. 

Proof. A fixed-point configuration requires that the left hand sides of 



equations (91), (92), (93), and (94) vanish, so we obtain the conditions 



(95), (96), (97), and (98). A relative equilibrium can be generated from 



a fixed-point configuration if and only if the right hand sides of (91) 



(|92j), (|93j), and (|94j) vanish as well, i.e. 

(kojV 2 + k(3 2 p 2 — a 2 )ri cosaj = 0, 

/ 2 2 , n2 2 2\ n 

[Ken r i J rKpp i —a )Ti smai = U, 
(na 2 r 2 + n(3 2 p 2 — [3 2 )pi cos 6« = 0, 

(Ka 2 r 2 + n(3 2 p 2 — (3 2 )pi sin&j = 0, 

where r 2 + p 2 = k~ x , i — 1, 2, . . . , n. Since there is no 7 E R such that 
sin 7 and cos 7 vanish simultaneously, the above expressions are zero in 
one of the following circumstances: 

(a) r« = 0, and consequently pi = k~ x I 2 , for all i £ {1,2, ... , n}, 

(b) pi = 0, and consequently r\ = k~ x I 2 for all « G {1,2,..., n}, 

(c) there is a proper subset J C {1,2, ... ,n} such that r$ = 
(consequently p, = ft -1 / 2 ) for all i E J and p 3 - = (consequently 
Tj = k- 1 ' 2 ) for all j E {1,2, . . . ,n}\ J, 

(d) Ka 2 r 2 + Kf3 2 p 2 -a 2 = na 2 r 2 + Kf3 2 p 2 - f3 2 = 0, i E {1, 2, ... , n}. 
Cases (a) and (b) correspond to simply rotating /t-positive relative 

equilibria, thus recovering condition (i) in Criterion [2] Case (c) cor- 
responds to (i) in the above statement. Since, from Definition [6j it 
follows that the frequencies a and (3 are nonzero, the identities in case 
(d) can obviously take place only if a 2 = (3 2 , i.e. \a\ = \(3\ 7^ 0, so 
(d) corresponds to condition (ii) in the above statement. This remark 
completes the proof. □ 

13.3. Criterion for K-negative elliptic relative equilibria. We 

further consider the motion in and start with proving a criterion for 
the existence of simply rotating K-negative elliptic relative equilibria. 
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Criterion 5. (^-negative elliptic relative equilibria) Consider the 
point particles of masses m\,m2, ■ ■ ■ ,m n > 0, n > 2, in EH 3 , i.e. for 
k < 0. Then system (44) admits solutions of the form (l79|).' 



q = (qi, q2, • • • , q«), q» = {w u %i, tu, Zi),i = i,2,...,n, 

Wiit) = n cos(at + en), Xi(t) = n sin(at + a*), yi(t) = y h Zi(t) = z h 

with wf + xf = rf, rf + yf — zf = k~ x , and yi, Zi constant, i = 1, 2, . . . , n, 
i.e. (simply rotating) K-negative elliptic relative equilibria, if and only 
if there are constants rj, Oj, yi, Zi, % = 1, 2, . . . , n, and a ^ 0, such that 
the following An conditions are satisfied: 

rnJ/d 3 / 2 ^,- cosgu — «ei,Ti cosaj) . 9 . 9 

m,-|/c| 3 / 2 (r 7 - sin a,- — near, sina«) , 9 9 
(100) £ J ' ' r K 2 e 2 l 1)3 /2 = - l)«^sma,, 



J=l V ^ 



(ioi) E 



Ka 2 rjyi, 



(102) E 



-2 C 2 



i = l,2,...,n, where = r^cos^ - a_,-) + y^- - z^-, i,j 
1,2, ...,n, i ^ j. 



Proof. Consider a candidate q as above for a solution of system (44). 
Some straightforward computations show that 

tij ■= q* qj = nrj cos(ai - dj) + yiy-j - ZiZj, i,j = 1, 2 . . . , n, i ^ j, 

q { © ^ = a 2 r 2 , j = 1,2, . . . ,n, 

u)i = —a 2 r-i cos(at + a*), Xj = — aVj sin(at + a*), 

j/i = 2j = 0, i = 1,2, . . . ,n. 

Substituting q and the above expressions into the equations of motion 
(44), for the w coordinates we obtain conditions involving cos(at + ai), 
whereas for the x coordinates we obtain conditions involving sin(at + 
aj). In the former case, using the fact that cos(at + aj) = cos at cos dj — 
sin at sin a,, we can split each equation in two, one involving cos at 
and the other sin at as factors. The same thing happens in the latter 
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case if we use the formula sin (at + a*) = sin at cos a» + cos at sin 



Each of these equations are satisfied if and only if conditions (99) and 



(100) take place. Conditions (101) and (102) follow directly from the 



equations involving the coordinates y and z. This remark completes 
the proof. □ 

13.4. Criterion for K-negative hyperbolic relative equilibria. 

We continue our study of the hyperbolic space with proving a cri- 
terion that shows under what conditions simply rotating K-negative 
hyperbolic orbits exist. 

Criterion 6. (/{-negative hyperbolic relative equilibria) Con- 
sider the point particles of masses mi, 777.2, • • • , m n > 0, n > 2, in H^ 7 
i.e. for k < 0. Then the equations of motion (44) admit solutions of 



the form (80): 

q = (qi,q2, 



,qnj,qi 



[Wi, Xi } yi, Zij 



1,2, 



Wi(t) = Wi (constant), Xi(t) = Xi (constant), 

yi{t) = rji sinh(/3t + k), Zi(t) = T]i cosh(/3t + h), 

with yf — zf = —rjf, wf + xf — r/f = k" 1 , i = 1, 2, . . . , n, i.e. (simply 
rotating) K-negative hyperbolic relative equilibria, if and only if there are 
constants r]i,Wi,Xi, % — 1,2, . . . ,n, and (3^0, such that the following 
An conditions are satisfied: 

003) t ^te^ ^ 



(104) ± m ^^_^ = 



£ = 2 + 2 s . nh 

(106) £ '">l 3/2 ('^°^'' J _-p^.^h''.) = (K1;? + mshK 

i = 1,2, ... ,n, where //^ = WiWj + XiXj — r/ir/j cosh(6j — bj), i,j = 
1,2, ... ,n, iy^j. 
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Proof. Consider a candidate q as above for a solution of system (44). 
Some straightforward computations show that 

fiij := q 4 Hqj = WiWj + XiXj - T]iT]j cosh(^ - bj) , i, j = 1, 2, . . . , n, i j, 

qU Hq, = /3 2 Vi, % = l,2,...,n, 

Wi = = 0, 

jji = (3 2 r]i sinh(/3t + bi), z\ = fi 2 ^ cosh(/3i + bi), i = 1, 2, . . . , n. 
Substituting q and the above expressions into the equations of motion 



(44), we immediately obtain for the w and x coordinates the equations 



(103) and (104), respectively For the y and z coordinates we obtain 
conditions involving sinh(/3i + bi) and cosh(/3t + bi), respectively. In 
the former case, using the fact that sinh(/3t + bi) = sinh fit cosh 6, + 
cosh /3t sinh 6j, we can split each equation in two, one involving sinh fit 
and the other cosh fit as factors. The same thing happens in the later 
case if we use the formula cosh(/3t+6j) = cosh fit cosh 6j + sinh fit sinh 6j. 



Each of these conditions are satisfied if and only of conditions (105) 



and (106) take place. This remark completes the proof. □ 



13.5. Criterion for /t-negative elliptic-hyperbolic relative equi- 
libria. We end our study of existence criteria for relative equilibria in 
hyperbolic space with a result that shows under what conditions simply 
rotating k- negative elliptic- hyperbolic orbits exist. 

Criterion 7. (/{-negative elliptic-hyperbolic relative equilibria) 

Consider the point particles of masses m\,m,2, ■ ■ ■ , m n > 0, n > 2, in 



i.e. for k < 0. Then the equations of motion (44) admit solutions 



of the form (81): 



q = (qi, q2, ■ ■ ■ , q n ), q* = (w l , x { , y u z t ), i = l, 2, . . . , n, 

Wi(t) = Ti cos(at + cii), Xi(t) = Ti sin (at + a*), 

yi(t) = rji sinh(/3t + bi), z^t) = rji cosh(/3t + bi), 

i.e. (doubly rotating) K-negative elliptic-hyperbolic relative equilibria, if 
and only if there are constants r,- t , r/j, a^, bi, i — 1,2 ... , n, and a, fi ^ 0, 
such that the following An conditions are satisfied: 
(107) 

£ mM^r ^-^cosa,) = + _ ^ 



(i2 



Florin Diacu 



(108) 

•n 

E 

2=1 

(109) 

n 

E 



3/2, 



r j sin Oj 



K7^rj sin a. 



l) 3 / 2 



/ 2 2, o2 2 2\ 

[na r { + Kfi rji — a ) 



a m sin a,- 



rrij\K 



3/2 



(?7j sinh bj — K'jijTji sinh fe. 



(110) 

n 

E 

3=1 



!)3/2 



k« ^ +/3 )?7j sinhfoj 



3/2 



(?7j cosh 6j — Kjijrji cosh 6j 



( K 2 T 2 _ 1)3 /2 



z = 1, 2, . . . , n, where 7^ 
1,2, ...,n, i ^ j. 



r^rj cos^aj 



na r i +n/3 rj i +/i )r/j cosh bi 



rjiTjj cosh(bi - bj), i,j 



Proof. Consider a candidate q as above for a solution of system (44). 
Some straightforward computations show that 

7y := qjHqj = r^j cos(a i -a i )-^ j cosh(6 i -6 j ), z, j = 1, 2, . . . , n, i ^ j, 



q { B i\i = a z rf + /3\ 2 
— a 2 ri cos(at + a*), Xi 



-- 1,2,. ..,71, 

-aVj sin (at + 



j/i = (3 2 r]i sinh(/3t + ^ = fi 2 7]i cosh(/3i + i = 1, 2, 



n. 



Substituting these expression and those that define q into the equa- 
tions of motion (44), we obtain for the w and x coordinates conditions 
involving cos(at + aj) and sin(at + aj), respectively. In the former case, 
using the fact that cos(at + a{) = cos at cos dj — sin at sin a,, we can 
split each equation in two, one involving cos at and the other sin at as 
factors. The same thing happens in the latter case if we use the formula 
sin(at + af) = sin at cos aj + cos at sin dj. Each of these equations are 
satisfied if and only if conditions (107) and (108) take place. 

For the y and z coordinates we obtain conditions involving sinh(/3t + 
bi) and cosh(/3t+6j), respectively. In the former case, using the fact that 
sinh(/3t+6j) = sinh fit cosh 6j+cosh fit sinh bi, we can split each equation 
in two, one involving sinh fit and the other cosh fit as factors. The same 
thing happens in the later case if we use the formula cosh(/3t + bj) = 
cosh fit cosh bi + sinh fit sinh 6j. Each of these conditions are satisfied 
if and only of conditions (109) and (110) take place. This remark 
completes the proof. □ 
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13.6. Nonexistence of /t-negative parabolic orbits. The same as 
in the curved n-body problem restricted to H 2 , parabolic relative equi- 
libria do not exist in H^. The idea of the proof is similar to the one 
we used in [26]: it exploits the basic fact that a relative equilibrium 
of parabolic type would violate the conservation law of the angular 
momentum. Here are a formal statement and a proof of this result. 

Proposition 5. (Nonexistence of /t-negative parabolic relative 
equilibria) Consider the point particles of masses mi,m2, . . . ,m n > 
0,n > 2, in i.e. for k < 0. Then system (44) does not admit solu- 
tions of the form (82), which means that K-negative parabolic relative 



equilibria do not exist in the 3-dimensional curved n-body problem. 



Proof. Checking a solution of the form (82) into the last integral of 
(|43l), we obtain that 



C34 = ^rriiijjiZi-yiZi) = ^ 



rrii 



i=i 



i=l 



7< + fat + {5 { - 7^) 



t 1 



\Pi+{5i-ii)A 



E 



8=1 



Si + Pit + (5i - 7i 



[Pi + (Si - Ti )t] 



miPiili - $i) ~ ^1 m i{li - Sift, i e {1, 2, . . . , n}. 



i=l 



i=l 



Since C34 is constant, it follows that 7^ = 5 i? j £ {1,2,..., n}. But from 



(82) we obtain that af + Pf = k 1 < 0, a contradiction which proves 



that parabolic relative equilibria don't exist. This remark completes 
the proof. □ 



14. Qualitative behaviour of the relative equilibria in 

In this section we will describe some qualitative dynamical prop- 
erties for the K-positive elliptic and K-positive elliptic-elliptic relative 
equilibria, under the assumption that they exist. (Examples of such 
solutions will be given in Sections 17 and 18 for various values of n 
and of the masses mi, m.2, . . . , m n > 0.) For this purpose we start with 
some geometric-topologic considerations about S^. 

14.1. Some geometric topology in Consider the circle of radius 
r, in the wx plane of IR 4 and the circle of radius p, in the yz plane of 
IR 4 , with r 2 + p 2 = hT x . Then T 2 p is the cartesian product of these 
two circles, i.e. a 2-dimensional surface of genus 1, called a Clifford 
torus. Since these two circles are submanifolds embedded in IR 2 , T 2 is 
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embedded in IR 4 . But T 2 p also lies on the sphere S 3 , which has radius 
R = k~ 1 / 2 . Indeed, we can represent this torus as 

(111) T 2 p = {(w 7 x,y,z) | r 2 + p 2 = K'\0<6, ( f ) <27r}, 

where w = rcosO, x = rsin^, y = pcos0, and z = psin0, so the 
distance from the origin of the coordinate system to any point of the 
Clifford torus is 

(r 2 cos 2 # + r 2 sin 2 # + p 2 cos 2 + p 2 sin 2 0) 1/2 = (r 2 + p 2 ) 1/2 = k~ 1/2 = R. 

When r (and, consequently, p) takes all the values between and R, 
the family of Clifford tori such defined foliates § 3 (see Figure [lj. Each 
Clifford torus splits § 3 into two solid tori and forms the boundary 
between them. The two solid tori are congruent when r = p = R/y/2. 
For the sphere § 3 , this is the standard Heegaard splittingnof genus 1. 



Figure 1. A 3-dimensional projection of a 4- 
dimensional foliation of the sphere §^ into Clifford tori. 

Unlike regular tori embedded in IR 3 , Clifford tori have zero Gaussian 
curvature at every point. Their flatness is due to the existence of an 
additional dimension in IR 4 . Indeed, cylinders, obtained by pasting two 
opposite sides of a square, are flat surfaces both in IR 3 and IR 4 . But 
to form a torus by pasting the other two sides of the square, cylinders 
must be stretched in IR 3 . In IR 4 , the extra dimension allows pasting 
without stretching. 



6 A Heegaard splitting, named after the Danish mathematician Poul Heegaard 
(1871-1943), is a decomposition of a compact, connected, oriented 3-dimensional 
manifold along its boundary into two manifolds having the same genus g, with 
.9 = 0,1,2,... 
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14.2. Qualitative properties of the relative equilibria in S^. The 

above considerations allow us to state and prove the following result, 
under the assumption that K-positive elliptic and K-positive elliptic- 
elliptic relative equilibria exist in S^. 

Theorem 1. (Qualitative behaviour of the relative equilibria in 

E>D Assume that, in the 3-dimensional curved n-body problem, n> 3, of 
curvature k > 0, with bodies of masses mi,m2, . . . , m n > 0, K-positive 
elliptic and K-positive elliptic- elliptic relative equilibria exist. Then the 
corresponding solution q may have one of the following dynamical be- 
haviours: 



(i) If q is given by (77), the orbit is a (simply rotating) K-positive 
elliptic relative equilibrium, with the body of mass mi moving on a 
(not necessarily geodesic) circle Ci,i = 1,2, ... ,n, of a 2-dimensional 
sphere in in the hyperplanes wxy and wxz, the circles Ci are parallel 
with the plane wx; another possibility is that some bodies rotate on a 
great circle of a great sphere, while the other bodies stay fixed on a 
complementary great circle of another great sphere. 



(ii) If q is given by (78), the orbit is a (doubly rotating) K-positive 
elliptic- elliptic relative equilibrium, with some bodies rotating on a great 
circle of a great sphere and the other bodies rotating on a complemen- 
tary great circle of another great sphere; another possibility is that each 
body rrii is moving on the Clifford torus T^ ., i = 1,2, ... ,n. 

Proof, (i) The bodies move on circles Ci, i = 1, 2, . . . , n, because, by 



(77), the analytic expression of the orbit is given by 

q = (qi,q2,- • • ,Qn), q» = (w^x^y^Zi), i = l,2,...,n, 

Wi(t) = n cos(at + ai), Xi{t) = n sin(at + ai),yi{t) = y { , Zi{t) = z u 

with wf + xl = rf, rf + yf + zf = re -1 , and y iy Zi constant, i = 1, 2, . . . , n. 
This proves the first part of (i), except for the statements about par- 
allelism. 

In particular, if some bodies lie on the circle 

SU = {(°.°.2/»^ \y 2 + z 2 = K- 1 }, 

with yi(t) = yi = constant and Zi(t) — Z{ — constant, then the elliptic 
rotation, which changes the coordinates w and x, does not act on them, 
therefore the bodies don't move. This remark proves the second part 
of statement (i). 

To prove the parallelism statement from the first part of (i), let us 
first remark that, as the concept of two parallel lines makes sense only 
if the lines are contained in the same plane, the concept of two parallel 
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planes has meaning only if the planes are contained in the same 3- 
dimensional space. This explains our formulation of the statement. 
Towards proving it, notice first that 

n n 

C wx = ^2 m i( w iXi - ™i x i) = miT i 
i=l i=l 

and 

n 

i=i 

These constants are independent of the bodies' position, a fact that 
confirms that they result from first integrals. To determine the values 
of the constants c wyi c wz , c xy , and c xz , we first compute that 

n n 

i=l i=l 
n n 

Cwz = ^2 m i( w iZi ~ ™i z i) = a ^2 m i r i Z i sin ( at + a i), 
i=l i=l 
n n 

Cx V = ^2 m i( X iVi ~ X iVi) = a ^2 m » r *^» cos ( at + a i), 
i=l i=l 
n n 

Cxz = ^2 m i( X i Z i ~ X i Z i) = a ^2 m i r i Zj - COS ( at + a i) ■ 
i=l i=l 

Since they are constant, the first integrals must take the same value 
for the arguments t — and t = ir/a. But at t — 0, we obtain 

n n 

C wy = a^2 m i r iVi Sm a h C wy = a ^2 mir ' iZi Sil1 
i=l i=l 

n n 

c xy = a^2 m i r iVi cos a h c xz = a^2 m i r 'i z i cos a h 
1=1 1=1 

whereas at £ = 7t/q;, we obtain 

n n 

C wy = ~Oi^2 m » r »fi Sil1 a *> C ™f = m » r ^» Sil1 a *' 

i=l i=l 

n n 

c xy = -a^2 m i r iVi c °s a?, c xz = -a^2 m i r %Zi cos a,. 
i=i i=i 
Consequently, c„, y = c wz = c xy = c xz = 0. Since, as we already showed, 
c yz = 0, it follows that the only nonzero constant of the total angular 
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momentum is c wz . This means that the particle system has nonzero 
total rotation with respect to the origin only in the wx plane. 

To prove that the circles Ci,i = 1,2, ... ,n, are parallel with the plane 
wx in the hyperplanes wxy and wxz, assume that one circle, say C\, 
does not satisfy this property. Then some orthogonal projection of C\ 
(within either of the hyperplanes wxy and wxz) in at least one of the 
other base planes, say xy, is an ellipse, not a segment — as it would 
be otherwise. Then the angular momentum of the body of mass mi 
relative to the plane xy is nonzero. Should other circles have an elliptic 
projection in the plane xy, the angular momentum of the corresponding 
bodies would be nonzero as well. Moreover, all angular momenta would 
have the same sign because all bodies move in the same direction on 
the original circles. Consequently c xy 7^ 0, in contradiction with our 
previous findings. Therefore the circles Ci,i = 1,2, ... ,n, must be 
parallel, as stated. 

(ii) When a K-positive elliptic-elliptic (double) rotation acts on a 
system, if some bodies are on a great circle of a great sphere of S^, while 
other are on a complementary great circle of another great sphere, then 
the former bodies move only because of one rotation, while the latter 
bodies move only because of the other rotation. The special geometric 
properties of complementary circles leads to this kind of qualitative 
behaviour. 

To prove the other kind of qualitative behaviour, namely that the 
body of mass rrii of the doubly rotating K-positive elliptic-elliptic rel- 
ative equilibrium moves on the Clifford torus , i = 1,2, ... ,n, it 
is enough to compare the form of the orbit given in (78) with the 
characterization (111) of a Clifford torus. This remark completes the 
proof. □ 



14.3. Qualitative properties of the relative equilibria generated 
from fixed-point configurations in We will further outline the 
dynamical consequences of Criterion [2] and Criterion |4j under the as- 
sumption that K-positive elliptic and K-positive elliptic-elliptic relative 
equilibria, both generated from fixed-point configurations, exist in S^. 
This theorem deals with a subclass of the orbits whose qualitative be- 
haviour we have just described. 

Theorem 2. (Qualitative behaviour of the relative equilibria 
generated from fixed-point configurations in Consider in 
i.e. for k, > 0, the point particles of masses mi, m 2 , . . . , m n > 0, n > 2. 
Then a relative equilibrium q generated from a fixed point configuration 
may have one of the following characteristics: 
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(i) q is a (simply rotating) K-positive elliptic orbit for which all bodies 
rotate on the same great circle of a great sphere ofE>^.; 

(ii) q is a (simply rotating) K-positive elliptic orbit for which some 
bodies rotate on a great circle of a great sphere, while the other bodies 
are fixed on a complementary great circle of a different great sphere; 

(Hi) q is a (doubly rotating) K-positive elliptic- elliptic orbit for which 
some bodies rotate with frequency a 7^ on a great circle of a great 
sphere, while the other bodies rotate with frequency (3 ^ on a com- 
plementary great circle of a different sphere; the frequencies may be 
different in size, i.e. \a\ ^ \(3\; 

(iv) q is a (doubly rotating) K-positive elliptic- elliptic orbit with fre- 
quencies a, (5 7^ equal in size, i.e. \a\ = \(3\. 

Proof, (i) From conclusion (i) of Criterion [2j a (simply rotating) k- 
positive elliptic relative equilibrium of the form 

q = (qi, q2, q«), q* = (w i} x { , y h % = 1, 2, . . . , n, 

Wi = Ti cos(at + ai),Xi(t) = r, t sm[pd + a.^y^t) = y h z^t) = z u 

with wf+xf = rf , rf + yf + z 2 = k" 1 and yi, Z{ constant, i — 1, 2, . . . , n, 
generated from a fixed-point configuration, must satisfy one of two 
additional conditions (besides the initial 4n equations), the first of 
which translates into 

Ti = K~ 1/2 , Z = 1,2, ... ,71. 

This property implies that — Z{ — 0, i — 1, 2, . . . , n, so all bodies 
rotate along the same great circle of radius /t~ 1//2 , namely S 1 Kyz thus 
proving the statement in this case. 

(ii) From conclusion (ii) of Criterion |2j there is a proper subset X C 
{l,2,...,n} such that r« = for all i E X and rj = kT 1 ! 2 for all 
j E {l,2,...,n}\X. 

The bodies for which r, = must have Wi = Xi = and yf + z 2 = k^ 1 , 
so they are fixed on the great circle S\ wx , since yi and Zi are constant, 
j el, and no rotation acts on the coordinates w and x. 

As in the proof of (i) above, it follows that the bodies with rj = 
kT 1 ! 2 , j E {1, 2, . . . , n} \ X rotate on the circle S* , which is comple- 
mentary to wx so statement (ii) is also proved. 

(iii) From conclusion (i) of Criterion |4j a (doubly rotating) K-positive 
elliptic-elliptic relative equilibrium of the form 

q = (qi, Q2, • • • , q n ), q* = (m, x i} y u i = 1, 2, . . . , n, 



Wi(t) = Ti cos(at + cij), Xi(t) = Ti sin(at + a*), 

Ui{t) = Pi cos(/3t + bi), Zi{t) = pi sin(/3t + h), 
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with wf + x\ = rf, yf + zf = tf, rf + pf = re \ i = 1,2, ... ,n, 
generated from a fixed-point configuration, must satisfy one of two 
additional conditions (besides the initial 4n equations), the first of 
which says that there is a proper subset J C {1,2, ... ,n} such that 
n = for all z e J and pj = for all j G {1, 2, . . . , nj \ J. But 
this means that the bodies m.j with i & J have Wj = = and 
y 2 + zf = pf, so one rotation acts along the great circle wx , while the 
bodies with m%,i E. {1,2, ... ,n}\J' satisfy the conditions wf + xf = rf 
and yi = Zi = 0, so the other rotation acts on them along the great circle 
S^yz, which is complementary to S^ jU , x . Moreover, since the bodies are 
distributed on two complementary circles, there are no constraints on 
the frequencies a,/3 ^ 0, so they can be independent of each other, a 
remark that proves the statement. 

(iv) From statement (d) in the proof of Criterion |1J a (doubly rotat- 
ing) ^-positive elliptic-elliptic relative equilibrium may exist also when 
the bodies are not necessarily on complementary circles but the fre- 
quencies satisfy the condition \a\ = \/3\, that concludes the last 
statement of this result. □ 

15. Qualitative behaviour of the relative equilibria in M 3 k 

In this section we will describe some qualitative dynamical properties 
of the K-negative elliptic, hyperbolic, and elliptic-hyperbolic relative 
equilibria, under the assumption that they exist. (Examples of such 
solutions, which prove their existence, will be given in Sections 19, 
20, and 21.) For this purpose we start with some geometric-topologic 
considerations about ET^. 

15.1. Some geometric topology in H^. Usually, compact higher- 
dimensional manifolds have a richer geometry than non-compact man- 
ifolds of the same dimension. This is also true about if compared to 
H^. Nevertheless, we will be able to characterize the relative equilibria 
of in geometric-topologic terms. 

The surface we are introducing in this section, which will play for 
our dynamical analysis in the same role as the Clifford torus in S>^, 
is homoeomorphic to a cylinder. Consider a circle of radius r in the 
wx plane of IR 4 and the upper branch of the hyperbola r 2 — rf = kT 1 
in the yz plane of M 4 . Then we will call the surface C 2 ^ obtained 
by taking the cartesian product between the circle and the hyperbola 
a hyperbolic cylinder since it surrounds equidistantly a branch of a 
geodesic hyperbola in Hf^. Indeed, we can represent this cylinder as 

(112) C 2 ={(w,x,y,z) \r 2 -r, 2 = K-\ < 9 < 2tt, f e R}, 
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where w = rcosO, x = rsin^, y = T^sinli^, z = r\ cosh £. But the hy- 
perbolic cylinder Q 2 rr) also lies in HI 3 because the coordinates w, x, y, z, 
endowed with the Lorentz metric, satisfy the equation 

2,2,2 2 2 2 -1 

W +x +y —z =T —7] = K . 

As in the case of S> 3 , which is foliated by a family of Clifford tori, HI 3 
can be foliated by a family of hyperbolic cylinders. The foliation is, of 
course, not unique. But unlike the Clifford tori of 1R 4 , the hyperbolic 
cylinders of 1R 3,1 are not flat surfaces. In general, they have constant 
positive curvature, which varies with the size of the cylinder, becoming 
zero only when the cylinder degenerates into a geodesic hyperbola. 

15.2. Qualitative properties of the relative equilibria in HI 3 . 
The above considerations allow us to state and prove the following 
result, under the assumption that ^-negative elliptic, hyperbolic, and 
elliptic- hyperbolic relative equilibria exist. Notice that, on one hand, 
due to the absence of complementary circles, and, on the other hand, 
the absence of fixed points in HI 3 ., the dynamical behaviour of relative 
equilibria is less complicated than in S 3 . 

Theorem 3. (Qualitative behaviour of the relative equilibria 

in HI 3 ) In the 3- dimensional curved n-body problem, n > 3, of curva- 
ture k < 0, with bodies of masses m 1 ,m 2 , ■ ■ ■ ,m n > 0, every relative 
equilibrium q has one of the following potential behaviours: 



(i) if q is given by (79), the orbit is a (simply rotating) K-negative 
elliptic relative equilibrium, with the body of mass mi moving on a 
circle C l , i = 1, 2, . . . , n, of a 2-dimensional hyperboloid in HI 3 ; in the 
hyperplanes wxy and wxz, the planes of the circles C l are parallel with 
the plane wx; 



(ii) if q is given by (80), the orbit is a (simply rotating) K-negative 
hyperbolic relative equilibrium, with the body of mass rrii moving on 
some (not necessarily geodesic) hyperbola Hi of a 2-dimensional hy- 
perboloid in H 3 , i = 1,2, ... ,n; in the hyperplanes wyz and xyz, the 
planes of the hyperbolas C l are parallel with the plane yz; 



(Hi) if q is given by (81), the orbit is a (doubly rotating) K-negative 



elliptic-hyperbolic relative equilibrium, with the body of mass mi moving 
on the hyperbolic cylinder Cf , i — 1, 2, . . . ,n. 

Proof, (i) The bodies move on circles, Ci, i — 1,2, . . . ,n, because, by 



(79), the analytic expression of the orbit is given by 

q = (qi, Q2, • • • , q»), q* = Oi, x h y h i = 1, 2, . . . , n, 
Wi(t) = r { cos(at + a*), x { {t) = r; sin(at + a*), y^t) = y h z { {t) = z { 
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with wf+xj = rf , rf+yf — zf = kT 1 , and jji, z% constant, i = 1, 2, . . . , n. 
The parallelism of the planes of the circles in the hyperplanes wxy and 
wxz follows exactly as in the proof of part (i) of Theorem [TJ using the 
integrals of the total angular momenta. 

(ii) The bodies move on hyperbolas, Hi, i = 1, 2, . . . , n, because, by 
(80), the analytic expression of the orbit is given by 

q = (qi, Q2, • • , q«), q* = {w t , x i: y i} Zi), i = 1, 2, . . . , n, 

Wi(t) = Wi (constant), %i(t) = Xi (constant), 

Vi(t) = rji sinh(/3t + bi), Zi{t) = rji cosh(/3t + bi), 

with yf - zf = -rjf, wf + x?-rft = K~ 1 , i = 1, 2, . . . , n. 

Let us now prove the parallelism statement for the planes containing 
the hyperbolas Hi. For this purpose, notice that 

n 

C wx = ^ m i{ W i X i - ™i X i) = 
i=l 

and 

n n 

C V Z = /^ijy^i ~ Vi Z i) = -Py^^iVi- 
i=l i=l 

These constants are independent of the bodies' position, a fact that 
confirms that they result from first integrals. To determine the values 
of the constants c wy , c wz , c xy , and c xz , we first compute that 



Cwy 


i=i 


Wiiji 


- WiVi) 


= f3 rriiWirii cosh(/3t + bi), 
i=i 


Cwz 


n 
i=l 


[WiZi 


- WiZi) 


n 

= P/, miWirji sinh(/3t + bi), 
i=i 


Cxy 


n 
i=l 


{xiVi 


- XiVi) 


n 

= py^ j rriiXirii cosh(/3t + bi), 

i=l 


Cxz 


n 

= 


\XiZi 


X%Zi) 


n 

= (3 2j rriiXiTji sinh(/3t + bi). 
i=i 



We next show that c wy = 0. For this, notice first that, using the 
formula cosh(/3t + bi) = cosh bi cosh {3t + sinh bi sinh {3t, we can write 

(113) c wy = f3[A(t) + B(t)}, 

where 

n 

A(t) = mjWjrii cosh bj cosh {3t, 



i=l 
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and 

n 

B{t) = niiWirji sinh 6, sinh j3t. 
»=i 

But the function cosh is even, whereas sinh is odd. Therefore A is even 
and .B is odd. Since c wy is constant, we also have 

(114) c wy = 0[A(-t) + B(-t)} = 0[A{t) - B(t)}. 



From (113) and (114) and the fact that /3 ^ 0, we can conclude that 

c wy = f3A(t) and B(t) = 0, 

so jj.B{t) = 0. However, j^B{t) = (3A(t), which proves that c wy = 0. 

The fact that c xy = can be proved exactly the same way. The 
only difference when proving that c wz = and c xz = is the use of 
the corresponding hyperbolic formula, sinh(/3t + bj) = sinh /3t cosh&j + 
cosh (3t sinh b\. In conclusion, 

Cwx Cwy Cwz c xy c xz and c yz ~/~ 0? 

which means that the hyperbolic rotation takes place relative to the 
origin of the coordinate system solely with respect to the plane yz. 

Using a similar reasoning as in the proof of (i) for Theorem [TJ it can 
be shown that the above conclusion proves the parallelism of the planes 
that contain the hyperbolas % in the 3-dimensional hyperplanes wyz 
and xyz. 

(iii) To prove that (doubly rotating) K-negative elliptic-hyperbolic 
solutions move on hyperbolic cylinders, it is enough to compare the 
form of the orbit given in (80) with the characterization (112) of a 
hyperbolic cylinder. □ 
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Part 5. EXAMPLES 

Since Theorems [TJ |2j and [3] provide us with the qualitative behaviour 
of all the five classes of the relative equilibria that we expect to find 
in §^ and H^, we know what kind of rigid-body- type orbits to look 
for in the curved n-body problem for various values of n > 3. Ideal, 
of course, would be to find them all, but this problem appears to be 
very difficult, and it might never be completely solved. As a first step 
towards this (perhaps unreachable) goal, we will show that each type 
of orbit described in the above criteria and theorems exists for some 
values of n and mi, m 2 , . . . , m n > 0. 

To appreciate the difficulty of the above mentioned problem, we re- 
mark that its Euclidean analogue is that of finding all central configu- 
rations for the Newtonian potential. The notoreity of this problem has 
been recognized for at least seven decades, [82]. In fact, we don't even 
know whether, for some given masses m 1; m 2 , . . . , m n > 0, with n > 5, 
the number of classes of central configurations (after we factorize the 
central configurations by size and rotation) is finite or notQ and, if it is 
infinite, whether the set of classes of central configurations is discrete 
or contains a continuum. The finiteness of the number of classes of cen- 
tral configurations is Problem 6 on Steven Smale's list of mathematics 
problems for the 21st century, [80]. Its analogue in our case would 
be that of deciding whether, for given masses, mi,m2, • • • ,m n > 0, 
the number of classes of relative equilibria of the 3-dimensional curved 
n-body problem is finite or not. 

16. Examples of ^-positive elliptic relative equilibria 

In this section, we will provide specific examples of K-positive ellip- 
tic relative equilibria, i.e. orbits on the sphere §^ that have a single 
rotation. The first example is that of a 3-body problem in which the 3 
equal masses are at the vertices of an equilateral triangle that rotates 
along a not necessarily great circle of a great or non- great sphere. The 
second example is that of a 3-body problem in which the 3 non-equal 
masses move at the vertices of an acute scalene triangle along a great 
circle of a great sphere. The third example is that of a relative equilib- 
rium generated from a fixed-point configuration in a 6-body problem 
of equal masses for which 3 equal masses move along a great circle of 
a great sphere at the vertices of an equilateral triangle, while the other 
3 masses, which are equal to the others, are fixed on a complementary 

7 Alain Albouy and Vadim Kaloshin have recently proved the finiteness of central 
configurations in the planar 5-body problem, except for a negligible set of masses, 
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great circle of another great sphere at the vertices of an equilateral 
triangle. The fourth, and last example of the section, generalizes the 
third example to the case of acute scalene, not necessarily congruent, 
triangles and non-equal masses. 

16.1. A class of ^-positive elliptic relative equilibrium of equal 
masses moving on a 2-dimensional sphere. In the light of Re- 
mark [2J we expect to find solutions in § 3 that move on 2-dimensional 
spheres. A simple example is that of the Lagrangian solutions (i.e. 
equilateral triangles) of equal masses in the curved 3-body problem. 
Their existence, and the fact that they occur only when the masses are 
equal, was first proved in [26] and [27] . So, in this case, we have n = 3 
and mi = m 2 = m 3 =: m. The solution of the corresponding system 
(44) we are seeking is of the form 

(115) q = (qi,q 2 ,q 3 ), 



q 

= 

yi{t) = 

w 2 {t) = 

y 2 (f) = 

w 3 (t) = 

ys(t) = 

with r 2 + y 2 + z 2 = 
Criterion [IJ we have 



Wi , Xi , Xji , Z% J , I 



1,2,3, 



r cos ut, 


Xlit) 


= r sin cut, 


y (constant), 


z x {t) 


= z (constant), 


r cos(wt + 27r/3), 


x 2 (t) 


= r sm(ut + 27r/3), 


y (constant), 


z 2 (*) 


= z (constant), 


r cos(wt + 47r/3), 


x 3 (t) 


= r sm(ut + 47r/3), 


y (constant), 


23® 


= z (constant), 



K 



Consequently, for the equations occurring in 



2/1 = 2/2 : 



: r 2 = r 3 =: r, ai = 0, a 2 = 27r/3, a 3 = 47r/3, 
Vz ='■ V (constant), z\ = z 2 = Z3 =: z (constant). 



Substituting these values into the equations (83), (|84]), (|85 
obtain either identities or the equation 

o 8m 



861, we 



a 



y/3r 3 (4-3Kr 2 ) 3 / 2 ' 
Consequently, given m > 0, re (0, n~ l l 2 ), and y, z with r 2 + y 2 + z 2 



K 



" 1 , we can always find two frequencies, 

2 / 2m 
r 



v/3r(4 - 3/tr 2 ) 3 /2 



and a 2 = — - 



2m 



\/3r(4-3/tr 2 ) 3 / 2 ' 



such that system (44) has a solution of the form (115). The positive 
frequency corresponds to one sense of rotation, whereas the negative 
frequency corresponds to the opposite sense. 
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Notice that if r = k 1//2 , i.e. when the bodies move along a great 



circle of a great sphere, equations (83), (84), (85), (86) are identities 



for any a 6 t, so any frequency leads to a solution. This phenomenon 
happens because, under those circumstances, the motion is generated 
from a fixed-point configuration, a case in which we can apply Criterion 
[2j whose statement is independent of the frequency. 

The bodies move on the great circle S^ y2 of a great sphere only if 
y = z = 0. Otherwise they move on non-great circles of great or non- 
great spheres. So we can also interpret this example as existing in the 
light of Remark [3j which says that there are K-positive elliptic rotations 
that leave non-great spheres invariant. 

The constants of the angular momentum are 

Cwx 3ttik oj ^ and c W y c wz c X y Cxz ^yz 0? 

which means that the rotation takes place around the origin of the 
coordinate system only relative to the plane wx. 



16.2. Classes of K-positive elliptic relative equilibria of non- 
equal masses moving on a 2-dimensional sphere. It is natural 
to ask whether solutions such as the one in the previous example also 
exist for non-equal masses. The answer is positive, and it was first 
answered in |22j, where we proved that, given a 2-dimensional sphere 
and any triangle inscribed in a great circle of it (for instance inside the 
equator z — 0), there are masses, m 1; m 2 , m 3 > 0, such that the bodies 
form a relative equilibrium that rotates around the z-axis. 

We will consider a similar solution here in which moves on the 
great circle S\ yz of the great sphere S 2 or S 2 Z . The expected analytic 
expression of the solution depends on the shape of the triangle, i.e. it 
has the form 



(116) 





q = 


(qi,q2,qs), q* = 


(tfj, Xj, 




z = l,2,3, 


w 1 (t) 


= fiT 


^ 2 cos(u;t + ai), 


Xl(t) 


= K 


_1/2 sin(a;t + ai) 


yi(t) 


= 0, 




z x {t) 


= 0, 




w 2 (t) 


= K~ 


" 1//2 cos(o;t + a 2 ), 


x 2 (t) 


= K 


_1/2 sm(ujt + a 2 ) 


V2 (t) 


= o, 






= o, 




w 3 (t) 


= K~ 


^ 2 cos(a;t + a 3 ), 


x 3 (t) 


= K 


~ 1/2 sm(ut + a 3 ) 


y 3 (t) 


= o, 




zs(*) 


= o, 





where the constants 01,02,03, with < ai < a 2 < a% < 2ir, determine 
the triangle's shape. The other constants involved in the description 
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of this orbit are 
(117) 



n = r 2 = r 3 



We can use now Criterion |2| to prove that (116) is a K-positive elliptic 
relative equilibrium for any frequency w^O. Indeed, we know from [22] 
that, for any shape of the triangle, there exist masses that yield a fixed 
point on the great circle S* ya ., so the corresponding equations (87), 



(88), (89), (90) are satisfied. Since conditions (117) are also satisfied 



0. 



the proof that (116) is a solution of (44) is complete. 
The constants of the angular momentum integrals are 

c wx = {mi + m 2 + m^K^uj ^ and c wy = c wz = c xy = c xz = c yz 

which means that the bodies rotate in IR 4 around the origin of the 
coordinate system only relative to the plane xy. 

16.3. A class of /{-positive elliptic relative equilibria not con- 
tained on any 2-dimensional sphere. The following example of a 
(simply rotating) K-positive elliptic relative equilibrium in the curved 
6-body problem corresponds to the second type of orbit described in 
part (i) of Theorem [TJ and it is interesting from two points of view. 
First, it is an orbit that exists only in but cannot exist on any 
2-dimensional sphere. Second, 3 bodies of equal masses move on a 
great circle of a great sphere at the vertices of an equilateral triangle, 
while the other 3 bodies of masses equal to the first stay fixed on a 
complementary great circle of another great sphere at the vertices of 
an equilateral triangle. 

So consider the equal masses 

mi = 7712 = "^3 = "^4 = ^5 = ^6 ='■ 771 > 0. 

Then a solution as described above has the form 

q = (qi,q2,q3,q4,q5,qe), q* = (w^x^y^Zi), \ 



1,2, 



6. 



Wi 

ih 

V2 
W 3 
V3 



k 2 cos at, 
0, 

cos(at + 2tt/3), 

0, 

fiT 5 cos(at + 47r/3), 
0, 

0, 

_ i 

K 2 , 



Xi 


= K 


Zl 


= o, 


X 2 


= K 


Z 2 


= o, 


X3 


= K 


Z3 


= o, 


X4 


= 0, 


Zi 


= 0, 



k 2 sin at, 



sin(at + 27r/3), 



sin(at + 47r/3), 



Relative equilibria in the 3-dimensional curved n-body problem 



77 



/, ^ — w, 

_ 1 

K 2 

2/5 = n , z 5 



w 5 = 0, x 5 = 0, 

2 ' " J 2 ' 

w G = 0, x 6 = 0, 

\/3/c _ 2 

Ve = jr-, ^6 



2 7 2 
A straightforward computation shows that this attempted orbit, which 
is generated from a fixed-point configuration, satisfies Criterion[2j there- 



fore it is indeed a solution of system (44 ) forn = 6 and for any frequency 
a ^ 0. 

The constants of the angular momentum are 

Cwx Sinn ol ^ and c W y c wz Cxy Cxz ^yz 0, 

which implies that the rotation takes place around the origin of the 
coordinate system only relative to the wx plane. 

16.4. Classes of K-positive elliptic relative equilibria with non- 
equal masses not contained on any 2-dimensional sphere. To 

prove the existence of K-positive elliptic relative equilibria with non- 
equal masses not contained in any 2-dimensional sphere, we have only 



to combine the ideas of Subsections |16.2| and |16.3| More precisely, 
we consider a 6-body problem in which 3 bodies of non-equal masses, 
m 1 ,m 2 ,m 3 > 0, rotate on a great circle (lying, say, in the plane wx) 
of a great sphere at the vertices of an acute scalene triangle, while 
the other 3 bodies of non-equal masses, m 4 ,m 5 ,m 6 > 0, are fixed 
on a complementary great circle of another great sphere (lying, as a 
consequence of our choice of the previous circle, in the plane yz) at 
the vertices of another acute scalene triangle, which is not necessarily 
congruent with the first. 

Notice that, as shown in [22], we must first choose the shapes of the 
triangles and then determine the masses that correspond to them, not 
the other way around. The reason for proceeding in this order is that 
not any 3 positive masses can rotate along a great circle of a great 



sphere. Like the solution offered in Subsection 16.3, this orbit exists 
only in S^, but cannot exist on any 2-dimensional sphere. Its analytic 
expression depends on the shapes of the two triangles, i.e. 

(118) q = (qi,q 2 ,q3,q4,q5,q6), qi = (wi,Xi,yi,Zi), i = 1,2, 6, 

Wx = k^ 5 cos(at + <2i), x\ = K~i sin(a;t + ai), 

Vl = 0, zi = 0, 
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w 2 


= K 


2 cos(at + a 2 ), 


x 2 


= K 


"2 sin(at + a 2 ), 


V2 


= o, 




Z 2 


= 0, 




W 3 


= k~ 


2 cos (at + a 3 ), 


x 3 


= K~ 


2 sin(at + a 3 ), 


2/3 


= o, 




z 3 


= 0, 




Wi 


= 0, 




XA 


= 0, 




Va 


= K~ 


5 COS 64, 


Za 


= K~ 


5 sin 64, 


w 5 


= 0, 




x 5 


= 0, 




2/5 


= K~ 


5 COS 65, 


Z5 


= K~ 


2 sin6 5 , 


we 


= o, 




x 6 


= 0, 




V6 


= K~ 


5 COS 6 6) 


z% 


= Kf 


5 sin6 6 . 



where the constants ai,a 2 , and a 3 , with < a x < a 2 < a 3 < 2n, and 
64, 65, and 6 6 , with < 64 < b 5 < b 6 < 2n, determine the shape of the 
first and second triangle, respectively. 

For t — 0, we obtain the configuration given by the coordinates 



Wi 


= K~ 


2 cosai, 


Xi 


= K~ 


2 sin ai, 


yi 


= 0, 




Z\ 


= 0, 




w 2 


= K~ 


1 

'2 cosa 2 , 


x 2 


= K~ 


1 

'2 sin a 2 , 


1/2 


= 0, 




z 2 


= 0, 




w 3 


= K~ 


1 

2 cosa 3 , 


x 3 


= K~ 


1 

2 sin a 3 , 


2/3 


= 0, 




Z3 


= 0, 




WA 


= 0, 




XA 


= 0, 




1/4 


= K~ 


5 COS 64, 


za 


= K~ 


5 sin 64, 


w 5 


= 0, 




x 5 


= 0, 




2/5 


= K~ 


^ COS 65, 


z 5 


= K~ 


2 sin 65, 


Wq 


= 0, 




Xq 


= 0, 




2/6 


= K~ 


2 COS &6, 


za 


= K~ 


2 sin 6 6 . 



We prove next that this is a fixed-point configuration. For this purpose, 
let us first compute that 





V\ 2 


= v 2X 


= COs((2i — 


a 2 ), 






V\ 3 


= v 3X 


= cos(ai — 


«3), 






v 23 


= v 32 


= cos(a 2 — 


03), 




via 


= Va\ ■- 


= v\h -- 


= ^51 = ^16 


= ^61 


= 0, 


v 2 a 


= Va 2 ■- 




= ^52 = ^26 


= ^62 


= 0, 


v-za 


= Vaz - 


= v 3h - 


- ^53 — ^36 


= ^63 


= 0, 




VA5 


= ^54 


= cos(6 4 — 


05), 






VA6 


= ^64 


= COS (64 — 


b 6 ), 






"56 


= ^65 


= C0S(6 5 — 


be). 
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Since y\ = y% = y 3 = z\ = z 2 = Z3 = 0, it follows that, at t = 0, the 
equations involving the force for the coordinates W\, x\, w^, x 2 , W3, X3 
involve only the constants ra\, m 2 , 777.3, &\, ^2, 03- In other words, for 
these coordinates, the forces acting on the masses m\,m 2 , and m 3 do 
not involve the masses m 4 ,m 5 , and m 6 . But the bodies m\,m 2 , and 
m 3 lie on the great circle S* ^, which can be seen as lying on the 
great sphere S^. This means that, by applying the result of [22], 
the bodies mi, rri2, and 777.3 form an independent fixed-point configura- 
tion. Similarly, we can show that the masses morris, and m& form an 
independent fixed-point configuration. Therefore all 6 bodies form a 
fixed-point configuration. 

Consequently, we can use Criterion [2] to check whether q given by 



(118) is a K-positive elliptic relative equilibrium generated from a fixed- 
point configuration. We can approach this problem in two ways. One 
is computational, and it consists of using the fact that the positions at 
t = form a fixed-point configuration to determine the relationships 
between the constants mi, m 2 , 777.3, cti, a 2 , and a 3 , on one hand, and the 
constants 777.4, Ta^i m 6, b\, 65, and be, on the other hand. It turns out 



that they reduce to conditions (87), (88), (89), and (90). Then we only 
need to remark that 

r i — r 2 — r 3 — k" 1 ^ 2 and r 4 = r 5 = r 6 = 0, 

which means that condition (ii) of Criterion [2] is satisfied. The other 
approach is to invoke again the result of [22] and a reasoning similar 
to the one we used to show that the position at t = is a fixed-point 



configuration. Both help us conclude that q given by ( 118 ) is a solution 



of system (44) for any 

The constants of the angular momentum are 

c W x = (rri! + m 2 + m 3 )/t _1 a ^ and c wy = c wz = c xy = c xz = c yz = 0, 

which means, as expected, that the bodies rotate in M 4 around the 
origin of the coordinate system only relative to the plane wx. 

17. Examples of k-positive elliptic-elliptic relative 

equilibria 

In this section we construct examples of /t-positive elliptic-elliptic 
relative equilibria, i.e. orbits with two elliptic rotations on the sphere 
S^. The first example is that of a 3-body problem in which 3 bodies 
of equal masses are at the vertices of an equilateral triangle, which 
has two rotations of the same frequency. The second example is that 
of a 4-body problem in which 4 equal masses are at the vertices of a 
regular tetrahedron, which has two rotations of the same frequency. 
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The third example is that of 5 equal masses lying at the vertices of a 
pentatope with double rotation. This is the only case of the a regular 
polytope that allows relative equilibria, because the 5 other existing 
regular polytopes of M 4 have antipodal vertices, so they introduce sin- 
gularities. As in the previous example, this motion cannot take place 
on any 2-dimensional sphere. The fourth example is that of a 6-body 
problem, with 3 bodies of equal masses rotating at the vertices of an 
equilateral triangle along a great circle of a great sphere, while the other 
3 bodies, of the same masses as the others, rotate at the vertices of an 
equilateral triangle along a complementary great circle of another great 
sphere. The frequencies of the two rotations are distinct, in general. 
The fifth example generalizes the fourth example in the sense that the 
triangles are scalene, acute, not necessarily congruent, and the masses 
as well as the frequencies of the rotations are distinct, in general. 

17.1. Equilateral triangle with equal masses moving with equal 
frequencies. The example we will now construct is that of a (doubly 
rotating) ^-positive elliptic-elliptic equilateral triangle of equal masses 
in the curved 3-body problem in S>^ for which the rotations have the 
same frequency. Such solutions cannot be found on 2-dimensional 
spheres. So we consider k > and the masses mi = m 2 = m 3 =: 



m > 0. Then the solution we check for system (44) with n = 3 has the 
form: 

(119) q = (qi, q 2 , q 3 ), q* = (w it x h y h Zi), i = 1, 2, 3, 



w x = 


- r cos at, 


Xx 


= r sin at, 


yi = 


- p cos at, 


Z\ 


= p sin at, 


w 2 = 


- r cos (at + 27r/3) , 


x 2 


= r sin (at + 27r/3) , 


V2 = 


- pcos(at + 27r/3), 


Z2 


= p sin(at + 27r/3), 


W 3 = 


: r cos (at + 47r/3) , 


^3 


= r sin (at + 47r/3) , 


2/3 = 


- pcos(at + 47r/3), 


23 


= p sin(at + 47r/3), 


with r 2 + p 2 
the 3 bodies 


= k . For t = 0, the above 
the coordinates 


attempted solution gives for 


Wi = r, 


xi = 0, 


yi = 


= p, z x = 0, 


w 2 = - 


r r-\/3 
2' X2 ~ 2 ' 


2/2 = 


p pv^ 
"2' Z2 ~ 2 ' 



r 



'\/3 p pV3 



w 3 = -r, x 3 = — , y 3 = --, 
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which is a fixed-point configuration, since the bodies have equal masses 
and are at the vertices of an equilateral triangle inscribed in a great 
circle of a great sphere. Consequently, we can use Criterion [4] to check 



whether a solution of the form (119) satisfies system (44) for any a^O. 



A straightforward computation shows that the first 4n conditions are 
satisfied. Moreover, since the two rotations have the same frequency, 
it follows that condition (ii) of Criterion El is verified, therefore (119) is 



indeed a solution of system (44) for any q^O. 



The angular momentum constants are 



3mar , c wy 



0, c v 



-xy 



Smarp, c xz = 0, 



-yz 



3marp, 
- 3map 2 . 



which show that rotations around the origin of the coordinate system 
take place relative to 4 planes: wx,wz,xy, and yz. Consequently the 
bodies don't move on the same circle, but on the same Clifford torus, 
namely T 2 , that agrees with the qualitative result described in 

part (ii) of Theorem [T] 

Notice that for r = k^ 1 and p = 0, the orbit becomes a (simply 
rotating) K-positive elliptic relative equilibrium that rotates along a 
great circle of a great sphere in S^, i.e. an orbit such as the one we 
described in Subsection 116.11 



17.2. Regular tetrahedron with equal masses moving with e- 
qual-size frequencies. We will further construct a (doubly rotating) 
K-positive elliptic-elliptic solution of the 4-body problem in in which 
4 equal masses are at the vertices of a regular tetrahedron that has 
rotations of equal frequencies. So let us fix k > and m\ = mi = = 
777,4 —'■ Tn > and consider the initial position of the 4 bodies to be 



given as in the first example of Subsection 12.2, i.e. by the coordinates 





= 0, 


d 1 


= o, 




o, 




4 = 






w° 2 


= 0, 


• r 2 


= 0, 


vl = 


K 

3 


1/2 

? 


4 = 


1 

K 

3 


-1/2 


w% 


= 0, 


<r° 
3 


3 


vl = 


K 

3 


-1/2 

: 




1 

K 

3 


-1/2 


wl 


= 0, 


X 4 


3 


vl = 


K 

3 


-1/2 

j 


z A — 


1 

K 

3 


-1/2 



which is a fixed-point configuration. Indeed, the masses are equal and 
the bodies are at the vertices of a regular tetrahedron inscribed in a 
great sphere of S^. 
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For this choice of initial positions, we can compute that 
T\ = r 2 = 0, pi = p 2 = fiT 1 2 



■ 



V6 _ 1/2 _ 1/2 

which means that we expect that 2 bodies (corresponding to m,\ and 
m 2 ) move on the Clifford torus with r = and p = k^ 1 / 2 (i.e. the only 
Clifford torus, in the class of a given foliation of which is also a 
great circle of S\, see Figure fl]), while we expect the other 2 bodies to 
move on the Clifford torus with r = ^k" 1 / 2 and p = ^/t -1 / 2 . 

These considerations allow us to obtain the constants that determine 
the angles. Indeed, a± and a 2 can take any values, 

a 3 = 3n/2, a 4 = ir/2, b\ = ir/2, 

and b 2 ,b 3 , 6 4 are such that 

sm6 2 = --, cos& 2 = — ^— , 
cos 6 3 = — — , sm& 3 = — — 

COSO4 = — — , sm04 = — — 
o o 



which means that 63 = 64. 

We can now compute the form of the candidate for a solution gener- 
ated from the above fixed-point configuration. Using the above values 
of Ti, pi, cii and bi, i — 1,2, 3, 4, we obtain from the equations 

w° = n cos a», x°i = n sin a i; y° = p { cos 6 i; z° = sin & i( z = 1, 2, 3, 4, 

that the candidate for a solution is given by 

(120) q = (qi, q 2 , q 3 , q 4 ), qi = (u>i, %h Vh z i), i = 1, 2 , 3, 4, 

t«i = 0, x\ — 0, 

yi = fiT 1//2 cos(at + vr/2), zi = k~ 1 ^ 2 sin(at + vr/2), 

w 2 = 0, x 2 = 0, 

?/2 = k~ 1//2 cos(at + 62), z 2 = k -1 / 2 sin(at + b 2 ), 

\/G v/6 
w 3 = ^-fiT 1/2 cos(/3t + 3vr/2), x 3 = ^-fiT 1/2 sin(/3t + 3tt/2), 

2/3 = ^«r 1/2 cos(/3t + 63), ^3 = ^ «T 1/2 sin(/3t + 63), 
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W4 



2/4 



^"-^COsOft + TTytt), 



£4 



-AC 



a/3 

*~ 1/2 cos(/3t + £> 4 ), z 4 = ^-K- 1/2 sin(/3t + 6 4 ) 



VE 
T 

v/3 



1/2 



sin(/% + 7r/2), 



3 v 3 

If we invoke Criterion [4], do a straightforward computation, and use 
the fact that the frequencies of the two rotations have the same size, 



i.e. are equal in absolute value, we can conclude that q, given by (120) 



satisfies system (44), so it is indeed a solution of the curved 4-body 
problem for k > 0. 

Straightforward computations lead us to the following values of the 
angular momentum constants: 

4 8 



for /3 



ix ^TTLQK , C W y 

a and 
- i -1 

Cwx ^TYlOtK, , C W y 



^xy 



-xy 



Cxz 



0, c. 



ii'- 



0, Cy Z 



-man 



-man, 



for (3 = —a, a fact which shows that rotations around the origin of the 
coordinate system takes place only relative to the planes vox and yz. 

17.3. Regular pentatope with equal masses and equal-size fre- 
quencies. We will next construct a (doubly rotating) K-positive elliptic- 
elliptic solution of the 5-body problem in S^, in which 5 equal masses 
are at the vertices of a regular pentatope that has two rotations of 
equal-size frequencies. So let us fix k > and m 1 = m 2 = m 3 = 
m 4 = m 5 =: m > and consider the initial position of the 5 bodies 



to be given as in the second example of Subsection 12.2, i.e. by the 
coordinates 

0, 



w° = 




~0 _ 

x l — 


0, 


vl = 


0. 


A 


w° 2 = 




T - 
x 2 — 


ViEk- 1 / 2 


y°2 = 


0. 


z\ 


4 ' 


4 


w° = 


K~ 1 ' 2 


4 = 




y° 3 = 




z 3 


4 ' 


4V3 ' 












vl = 


VEk- 1 / 2 


z o 


4 ' 


4^ ' 


2V6 ' 






x 5 — 




yl = 


V5K- 1 / 2 


z o 

z 5 


4 ' 


4v/3 ' 


2VQ ' 



0. 

2V2 ' 

V5K- 1 / 2 

which is a fixed-point configuration because the masses are equal and 
the bodies are at the vertices of a regular pentatope inscribed in 
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For this choice of initial positions, we can compute that 
n = r 2 = k~ 1/2 , pi = p 2 = 0, 

r 3 = r 4 = r 5 = — — , p 3 = p 4 = p 5 = — -^=— , 

which means that we expect that 2 bodies (corresponding to m\ and 
m 2 ) move on the Clifford torus with r = kT 1 ! 2 and p = (i.e. the only 
Clifford torus, in a class of a given foliation of which is also a great 
circle of §^), while we expect the other 3 bodies to move on the Clifford 
torus with r = -^/t -1 / 2 and p = ^k~ 1 / 2 . 

These considerations allow us to obtain the constants that determine 
the angles. We obtain that 

di = 0, 

a 2 is such that 

f . Vl5 

cosa 2 = — -, sma 2 = — 

and a 3 , a 4 , a 5 are such that 



cos a 3 = — — , sin a 3 

V6 . VlO 

cosa 4 = , sma 4 = 

4 4 



\/6 . VTO 
cosa 5 = — — , sma 5 = — , 

which means that We further obtain that, since pi = 

p 2 = 0, the constants b\ and b 2 can be anything, in particular 0. Further 
computations lead us to the conclusion that 

b 1 =b 2 = b 3 = 0, h = 2tt/3, b 5 = An/3. 

We can now compute the form of the candidate for a solution gener- 
ated from the above fixed-point configuration. Using the above values 
of rj, pi, di and bi, i — 1,2, 3, 4, 5, we obtain from the equations 

u>° = Tj cos cij, a;° = sin a^, y° = p^ cos&j, z° = p^ sin 6j, i = 1,2, 3, 4, 5 

that the candidate for a solution is given by 

(121) q = (qi, q 2 , q 3 , q 4 ), q* = {w h x i: y h z t ), i = 1, 2, 3, 4, 5 

Wi = kT 1 ! 2 cos at, rri = /sT 1 / 2 sin at, 

Vi = 0, = 0, 
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w 2 = 


K- 1 ' 2 


V2 = 


0, 




1 




K 






2/3 — 


7!" 




1 


W4 = 


AC 

\/6 




v/5 . 


V4 = 


Vq 




1 


W 5 = 


Vq 






2/5 = 


\/6 



cos(at + 02), X2 = k 1 ^ 2 sin(at + a 2 ), 

Z2 = 0, 

~ 1 ^ 2 cos(at + a 3 ), x 3 = -=k -1 ' 2 sin(at + a 3 ), 

V6 

~ 1//2 cos (3t, z 3 = ^-=kT x I 2 sin /3t, 

V6 

_1 / 2 cos(at + 04), X4 = -zK -1 ' 2 sin(ai + 04), 

V6 

- 1/2 cos(/3t + 2tt/3), z 4 = ^fiT 1/2 sin(/3t + 2tt/3), 

V6 

~ 1//2 cos(at + 05), £5 = -zft -1 '' 2 sin(at + a 5 ), 

V6 

" 1/2 cos(/3t + 471-/3), z 5 = ^k- 1 / 2 sm(pt + 4tt/3), 

V6 

If we invoke Criterion [4j do a straightforward computation, and use the 
fact that the frequencies of the two rotations have the same size, i.e. 
I a I = |/3 1, we can conclude that q, given by (121), satisfies system (95), 
(96), (97), (98) and condition (ii), so it is indeed a (doubly rotating) 
/t-positive elliptic-elliptic relative equilibrium generated from a fixed 
point, i.e. a solution of the curved 5-body problem for k > and any 
value of a and with |a| = 7^ 0. 

A straightforward computation shows that the constants of the an- 
gular momentum are 

5-i _ _ _ _ n A -1 

Cwx ^TTlhZ O;, C W y Cwz Cxy Cxz U, Cyz ^TTLK (X 

for (3 = a and 

5 -1 _ n 5 -1 

Cwx —TTlfc Oi y C W y C wz C X y C xz U, Cy Z —TTlfc Oi 

for = —a, which means that the bodies rotate around the origin of 
the coordinate system only relative to the planes wx and yz. 

17.4. Pair of equilateral triangles with equal masses moving 
with distinct-size frequencies. We now construct an example in 
the 6-body problem in which 3 bodies of equal masses move along a 
great circle at the vertices of an equilateral triangle, while the other 3 
bodies of masses equal to those of the previous bodies move along a 
complementary circle of another great sphere, also at the vertices of an 
equilateral triangle. So consider k > 0, the equal masses mi = m 2 = 
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m 3 = m 4 = m 5 = m 6 =: m > 0, and the frequencies a, (3, which, in 
general, we can take as distinct, a ^ (3. Then a candidate for a solution 
as described above has the form 



q = 


= (qi 


,q2,q3,q4,q5,qe), 






yuZi), i = 1,2,.. 




= k~ 


i 

2 cos at, 


Xi = 




i 

2 sin at, 


Vi 


= 0, 




Zl = 


o. 




w 2 


= K~ 


2~cos(at + 27r/3), 


x 2 = 




^sin(at + 27r/3), 


2/2 


= o, 




Z 2 = 


o, 




w 3 


= k~ 


^ cos(at + 47r/3), 


%3 = 


K~ 


^ sin(at + 47r/3), 


2/3 


= 0, 




Z3 = 


o. 




u>4 


= o, 




Xi = 


o. 




2/4 


= k~ 


2 cos 


Za = 


K~ 


2 sin /3t, 




= 0, 




x 5 = 


0. 




2/5 


= K~ 


^cos(/3t + 27r/3), 


Z 5 = 


K~ 


'5sin(/3t + 27r/3), 


we 


= o, 




x 6 = 


o, 




2/6 




2 ~cos(/3t + 47r/3), 


z& = 


K~ 


3 sin(/3t + 4tt/3). 



For t = 0, we obtain the fixed-point configuration specific to S^, similar 
to the one constructed in Subsection 112.21 

Wi = 0, xi = kT 1/2 , yi = 0, zi = 0, 

w 2 = — , x 2 = , y 2 = 0, z 2 = 0, 

w 3 = — , x 3 = , y 3 = 0, z 3 = 0, 

u> 4 = 0, x 4 = 0, 2/4 = k" 1 / 2 , z 4 = 0, 

u> 5 = 0, x 5 = 0, y 5 = — , z 5 



2 ' J 2 ' 
w e = 0, X 6 = 0, |/ 6 = — , z 6 = . 



To prove that q given by (122) is a solution of system (44), we can 
therefore apply Criterion [4} A straightforward computation shows that 
the 4n conditions (95), (96), (97), (98) are satisfied, and then we can 
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observe that condition (i) is also verified because 

n = r 2 = r 3 = k~ 1/2 , Pl = p 2 = p 3 = 0, 

r 4 = r 5 = r 6 = 0, p 4 = Ps = Pe = fiT 1/2 . 



Consequently q given by (122) is a K-positive elliptic-elliptic relative 



equilibrium of the 6-body problem given by system (44) with n = 6 
for any a, (3 ^ 0. If a//3 is rational, that corresponds to a set 

of frequency pairs that has measure zero in IR 2 , the orbits are periodic. 
In general, however, a/ (3 is irrational, so the orbits are quasiperiodic. 
Though quasiperiodic relative equilibria were already discovered for the 
classical equations in M 4 , [2], [13], this is the first example of such an 
orbit in a 3-dimensional space. 

A straightforward computation shows that the constants of the an- 
gular momentum integrals are 

c-uox 3ttik (y ^ 0, Cy Z Sttik f3 ^ 0, c W y c wz c X y c xz 0, 

which means that the rotation takes place around the origin of the 
coordinate system only relative to the planes wx and yz. 

Notice that, in the light of [22], the kind of example constructed 
here in the 6-body problem can be easily generalized to any (n + m)- 
body problem, n,m > 3, of equal masses, in which n bodies rotate 
along a great circle of a great sphere at the vertices of a regular n-gon, 
while the other m bodies rotate along a complementary great circle of 
another great sphere at the vertices of a regular m-gon. The same as 
in the 6-body problem discussed here, the rotation takes place around 
the origin of the coordinate system only relative to 2 out of 6 reference 
planes. 

17.5. Pair of scalene triangles with non-equal masses moving 
with distinct frequencies. We will now extend the example con- 



structed in Subsection 17.4 to non-equal masses. The idea is the same 



as the one we used in Subsection 16.4 based on the results proved in 



according to which, given an acute scalene triangle inscribed in a 
great circle of a great sphere, we can find 3 masses such that this con- 
figuration forms a fixed point. The difference is that we don't keep the 
configuration fixed here by assigning zero initial velocities, but make it 
rotate uniformly, thus leading to a relative equilibrium. In fact, since 
we are in a 6-body problem, 3 bodies of non-equal masses will rotate 
along a great circle of a great sphere at the vertices of an acute scalene 
triangle, while the other 3 bodies will rotate along a complementary 
great circle of another great sphere at the vertices of another acute 
scalene triangle, not necessarily congruent with the other one. 
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So consider the masses mi, m 2 , m 3 , 777,4, ^5, me > 0, which, in general, 
are not equal. Then a candidate for a solution as described above has 
the form 



(123) q = (qi, q 2 , q 3 , q 4 , qs, Qe), Qi = ^i, 2/i, z<), i 



1,2, 



,6, 



U>1 


= ft 


2 cos(at + ai), 


£1 


= K 


"2 sin (at + ai), 


2/1 


= 0, 




Zi 


= 0, 






= k~ 


2 cos(at + a 2 ), 


£2 


= k~ 


"2 sin(at + a 2 ), 


2/2 


= 0, 




^2 


= 0, 




w 3 


= K~ 


2 cos(at + a 3 ), 


£3 


= «f 


"2 sin(at + a 3 ), 


2/3 


= 0, 




Z3 


= 0, 




u>4 


= 0, 




£4 


= 0, 




2/4 


= fiT 


2 cos(/3t + 64), 


^4 


= k~ 


"a sin(/3t + 64), 




= 0, 




£5 


= 0, 




2/5 


= K~ 


2 cos(/3t + 65), 


^5 


= K~ 


■3Bin(^ + 6 6 ), 




= 0, 




X 6 


= 0, 




2/6 


= K~ 


2 cos(/3t + fo 6 ), 




= ft" 


"a sin (/ft + 6 6 ), 



where the constants ai,a 2 , and a 3 , with < ai < a 2 < a 3 < 2tt, and 
61, 62, and 63, with < 64 < 65 < &6 < 27r, determine the shape of the 
first and second triangle, respectively. 

Notice that for t = 0, the position of the bodies is the fixed-point 



configuration described and proved to be as such in Subsection 16.4 



Therefore we can apply Criterion 4 to check whether q given in (123) 



is a K-positive elliptic-elliptic relative equilibrium. Again, as in Subsec- 



tion [16^ we can approach this problem in two ways. One is computa- 
tional, and it consists of using the fact that the positions at t — form 
a fixed-point configuration to determine the relationships between the 
constants mi, 777,2, rn^, 01, o>2, and a 3 , on one hand, and the constants 
m 4 , 777,5, m 6 , 64, 6 5 , and 6 6 , on the other hand. It turns out that they 



reduce to conditions (95), (96), (97), and (98). Then we only need to 
remark that 



ri = r 2 



r 3 



" 1//2 and r 4 = r 5 = r 6 = 0, 



which means that condition (i) of Criterion 4] is satisfied. The other 
approach is to invoke again the result of [22] and a reasoning similar 
to the one we used to show that the position at t = is a fixed-point 
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configuration. Both help us conclude that q given by ( 123 ) is a solution 



of system (44) for any a, f3 ^ 0. 

Again, when a/fl is rational, that corresponds to a negligible 

set of frequency pairs, the solutions are periodic. In the generic case, 
when a/ (3 is irrational, the solutions are quasiperiodic. 

A straightforward computation shows that the constants of the an- 
gular momentum integrals are 



(mi + m 2 + msjK 
, = 0, 



-i 



/3^0, 



^uiy '-'Miz ^xj/ ^x 

which means that the rotation takes place around the origin of the 
coordinate system only relative to the planes wx and yz. 



18. Examples of k-negative elliptic relative equilibria 
The class of examples we construct here is the analogue of the one 



presented in Subsection |16.1| in the case of the sphere, namely La- 
grangian solutions (i.e. equilateral triangles) of equal masses for k < 0. 
In the light of Remark |9j we expect that the bodies move on a 2- 
dimensional hyperboloid, whose curvature is not necessarily the same 
as the one of H^. The existence of these orbits, and the fact that they 
occur only when the masses are equal, was first proved in |26j and |27j . 
So, in this case, we have n = 3 and mi = m2 = m 3 =: m > 0. The 



solution of the corresponding system (44) we are seeking is of the form 



(124) q -- 


= (qi,q2,qs), = 


(tfj, y%i 


, z = l,2,3, 


Wl {t) = 


r cos ut, 


Xl(t) = 


r sin ut, 


yi{t) = 


y (constant), 


z x {t) = 


z (constant), 


w 2 (t) = 


r cos(wt + 27r/3), 


x 2 {t) = 


r sm(ut + 27r/3), 


y*{t) = 


y (constant), 


z 2 (t) = 


z (constant), 


Mt) = 


r cos(u;t + 47r/3), 


x 3 (t) = 


r sm(ut + 47r/3), 


ys(t) = 


y (constant), 


zs(t) = 


z (constant), 


with r 2 + y 2 — 


z 2 = n . Consequently, for the equations occurrin 


Criterion [5j we have 







r 1 = r 2 



y\ = y2 = 2/3 



= r 3 =: r, a x = 0, a 2 
-: y (constant) 



z\ = z 2 



2n/3,a 3 = 4vr/3, 

: z (constant). 



Z3 



Substituting these values into the equations (99), (100), (101), (102) 



we obtain either identities or the same equation as in Subsection 16.1 
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namely 



a 



2 om 



VSr 3 (A- 3/tr 2 ) 3 / 2 ' 

Consequently, given m > 0, r > 0, and y, z with r 2 + y 2 — z 2 = k~ x and 
z > |k|~ 1//2 , we can always find two frequencies, 



2m 2 / 2m 

<xi = —a — i= — and a 2 



V3r(A - 3/tr 2 ) 3 / 2 r V ^(4 - 3/tr 2 ) 3 / 2 : 



such that system (44) has a solution of the form (124). The positive 
frequency corresponds to one sense of rotation, whereas the negative 
frequency corresponds to the opposite sense. 

Notice that the bodies move on the 2-dimensional hyperboloid 

H Lj« = x ' z)\w 2 + x 2 - z 2 = fiT 1 - y%, y = constant}, 

which has curvature kq = — {y 2 — k -1 ) -1 / 2 . When y^ = 0, we have a 
great 2-dimensional hyperboloid, i.e. its curvature is k, so the motion 
is in agreement with Remark [6} 

A straightforward computation shows that the constants of the an- 
gular momentum are 

Cwx 3tjik ol 7^ and c W y c wz c X y Cxz ^yz 0? 

which means that the rotation takes place around the origin of the 
coordinate system only relative to the plane wx. 

19. Examples of k-negative hyperbolic relative equilibria 

In this section we will construct a class of ^-negative hyperbolic 
relative equilibria, for which, in agreement with Remark [7j the bodies 
rotate on a 2-dimensional hyperboloid of the same curvature as H 3 . 
In the 2-dimensional case, the existence of a similar orbit was already 
pointed out in [23], where we have also proved it to be unstable. So let 
us check a solution of the form 

(125) q = (qi, q 2 , q 3 ), q* = (Wi, x h y h z t ), % = 1, 2, 3, 

sinh j3t cosh (3t 

w 1 = 0, x 1 = 0, J/i = -i^ji7T' Zl = ^T/2-' 

w 2 = 0, x 2 = x (constant), y 2 = rj sinh (3t, z 2 = r) cosh j3t, 
W3 = 0, X3 = —x (constant), y$ = 77 sinh (3t, z% = 77 cosh f3t, 

with x 2 — f] 2 = K . Consequently 



Vi = \ K 



1/2 7 V2 = V3='-V (constant), 61 = b 2 = b 3 = 0. 



Relative equilibria in the 3-dimensional curved n-body problem 



91 



We then compute that 

A*ia = /"21 = yUi3 = ^23 = -M~ 1/2 ?7, ^23 = ^32 = ~ ^rj 2 . 

We can now use Criterion [6] to determine whether a candidate of the 
form q given by (125) is a (simply rotating) K-negative hyperbolic rela- 
tive equilibrium. Straightforward computations lead us from equations 
(103), (104), (105), and (106) either to identities or to the equation 

1 — Akt] 2 



A V 3 (\k\t] 2 - I) 3 / 2 ' 



Therefore, given k < and m,x,r] > with x 2 
two non-zero frequencies, 



K 



-1 



there exist 



0i 



1 



1 - 4kt] 2 



2r] y i](\K,\r] 2 — l) 3 / 2 



and /3 2 



1 



Aurj 2 



277 V 1](\k\t] 2 - l) 3 / 2 ' 



such that q given by (125) is a (simply rotating) K-positive hyperbolic 
relative equilibrium. Notice that the motion takes place on the 2- 
dimensional hyperboloid 



H 



{(0,x,y,z) \x +y - z = k }. 



A straightforward computation shows that the constants of the angular 
momentum are 



0, Cy Z 



-m/^fiT 1 + 2r] 2 



which means that the rotation takes place relative to the origin of the 
coordinate system only relative to the plane yz. 



20. Examples of k-negative elliptic-hyperbolic relative 

equilibria 

In this section we will construct a class of (doubly rotating) k- 
negative elliptic-hyperbolic relative equilibria. In the light of Remark 



1 1 , we expect that the motion cannot take place on any 2-dimensional 
hyperboloid of H 3 . In fact, as we know from Theorem |3j relative equi- 
libria of this type may rotate on hyperbolic cylinders, which is also the 
case with the solution we introduce here. 

Consider k < and the masses ni\ = m>2 = =: m > 0. We will 
check a solution of the form 



(126) q = (qi, q 2 , q 3 ), q* = (tUi, x u y h z { ), i = 1, 2, 3, 
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Wi = 0, X\ — 0, 

W2 = r cos at, X2 = r sin at, 
W3 = —r cos at, x 3 = — r sin at, 



Vi 



sinh {3t 
\^T 2 ' 
77 sinh (3t, 
77 sinh f3t, 



z 2 
z 3 



cosh fit 



K 



1/2 



rj cosh 
77 cosh (3t. 



In terms of the form (81) of a hyperbolic solution, (126) is realized 
when 



ri = 0, r 2 



?"3 



'/1 



AC 



-1/2 



0, 03 = 7T, and h\ 



,V2 

b 2 = 



= V3 



=■ V, 
0. 



Substituting these values into the equations (107), (108), (109), (110) 
of Criterion [7] and using the fact that r 2 — r]' 2 = k , we obtain the 
equation 

2 2 m(4| K |77 2 + l) 

P 4r? 3(_ Kr? 2_ 1 )3/2' 

for which there are infinitely many values of a and f3 that satisfy it. 
Therefore, for any k, < 0, masses m x = m 2 = m 3 =: m > 0, and r, rj 



with r 2 



there are infinitely many frequencies a and ft that 



correspond to a K-negative elliptic-hyperbolic relative equilibrium of 
the form (126). 

The bodies m 2 and m.3 move on the same hyperbolic cylinder, namely 
C 2 V , which has constant positive curvature, while mi moves on the 
degenerate hyperbolic cylinder C 2 | re |_ 1/2 , which is a geodesic hyper- 
bola, therefore has zero curvature. The motion is neither periodic nor 
quasiperiodic. 

A straightforward computation shows that the constants of the an- 
gular momentum are 



2mar 2 , c yz 



0. 



2/3?7 , C W y C WZ C X y 

which means that the rotation takes place around the origin of the 
coordinate system only relative to the wx and yz planes. 
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Part 6. CONCLUSIONS 

In this final part of the paper, we aim to emphasize 3 aspects related 
to the curved n-body problem. We will first analyze the issue of sta- 
bility of orbits and then present some open problems that arise from 
the results we obtained here in the direction of relative equilibria and 
the stability of periodic and quasiperiodic solutions. 

21. Stability 

An important problem to explore in celestial mechanics is that of 
the stability of solutions of the n-body problem. Unless the celestial 
orbits proved to exist have some stability properties, they will not be 
found in the universe. Of course, stability is only a necessary, but not 
a sufficient, condition for the existence of certain orbits in nature. If 
the solutions require additional properties, such as equal masses, then 
they are unlikely to show up in the universe either. Therefore some of 
the orbits mathematically proved to be stable can be observed in our 
planetary system, such as the Lagrangian relative equilibria formed by 
the Sun, Jupiter, and the Trojan asteroids, when one mass is very small. 
Other orbits, such as the figure eight solutions, p2], which require equal 
masses, have not been discovered in our solar system, although we have 
strong numerical evidence for their stability. No unstable orbit has ever 
been observed. 

From the mathematical point of view, it is important to distinguish 
between the several levels of stability that appear in the literature. 
Liapunov stability, for instance, is an unlikely property in celestial me- 
chanics. No expert in the field expects to find Liapunov stable orbits 
too often. The reason for this scepticism is that not even the classical 
elliptic orbits of the Euclidean Kepler problem posses this property. 
Indeed, think of two elliptic orbits that are close to each other, each 
of them accommodating a body of the same mass. Then the body on 
the larger ellipse moves slower than the other, so, after enough time 
elapses, the body moving on the inner orbit will get well ahead of the 
body moving on the outer orbit, thus proving that Liapunov stability 
cannot take place. Therefore the most we usually hope for in celestial 
mechanics from this point of view is orbital stability. However, most 
of the time we are happy to find linearly stable solutions. 

In October 2010, I asked Carles Simo whether he would like to 
investigate numerically the stability of the Lagrangian solutions of 
the curved 3-body problem. He accepted the challenge in January 
2011, and came up with much more than some numerical insight. To- 
gether with Regina Martinez, he studied the case of 3 equal masses, 



94 



Florin Diacu 



mi = m 2 = m 3 = 1, for the sphere § 2 of radius 1. In the end, they 
found an analytic proof of the following result. 

Theorem 4. (Martmez-Simo) Consider the Lagrangian solution, 
with masses nil = m 2 = m 3 = 1 of the 2- dimensional curved 3-body 



problem given by system (44) with n = 3 and k = 1, i.e. on the sphere 
§ 2 of radius 1, embedded in IR 3 , with a system of coordinates x,y,z. 
Let the motion take place on a non- geodesic circle of radius r in the 
plane z = constant, with < z = \J\ — r 2 < 1. Then the orbits are 
linearly stable (or totally elliptic) for r e (r 1; r 2 ) U (r 3 , 1) and linearly 
unstable for r e (0, ri) U (r 2 , r 3 ), where, with a 35-digit approximation, 

n = 0.55778526844099498188467226566148375, 

r 2 = 0.68145469725865414807206661241888645, 

r 3 = 0.92893280143637470996280353121615412. 

In the unstable domains, the local behaviour around the orbits consists 
of two elliptic planes and a complex saddle. At the values ri,r 2 , and 
r 3 , there occur Hamiltonian-Hopf bifurcations. 

In the proof of the above theorem, the only numerical computation 
is that of ri,r 2 , and r 3 , as roots of a certain characteristic polynomial. 
But the existence of these roots is easy to demonstrate, so the proof is 
entirely analytic, [BT] . 

As expected, when r is close to 0, i.e. when approaching the Eu- 
clidean plane, the Lagrangian orbit is unstable, as it happens with the 
classical Lagrangian solution of equal masses. In the Euclidean case, 
the Lagrangian orbits are stable only when one mass is very small. It 
is therefore surprising to discover two zones of stability, which occur in 
the intervals (ri, r 2 ) and (r 3 , 1), in the curved problem. This important 
result implies that the curvature of the space has decisive influence over 
the stability of orbits, a fact that was not previously known. 



22. Future perspectives 

More than 170 years after Bolyai and Lobachevsky suggested the 
study of gravitational motions in hyperbolic space, the recent emer- 
gence of the n-body problem in spaces of constant curvature, given 
by a system of differential equations set in a context that unifies the 
positive and the negative case, provides a new and widely open area 
of research. The results obtained in this paper raise several interesting 
questions, which I will outline below. 
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22.1. Relative equilibria. As shown here, the curved n-body prob- 
lem has 5 natural types of relative equilibria: 2 in and 3 in H^. Cri- 
teria 1 through 7 provide the necessary conditions for the existence of 
these relative equilibria. Let us consider the masses mi, m 2 , . . . , m n > 
0, n > 3. We have seen in the above examples that if the geomet- 
ric configuration of a relative equilibrium has a certain size, and the 
orbit is not generated from a fixed-point configuration, then we can 
find 2 frequencies (equal in size and of opposite signs, a and —a) in 
the case of orbits with a single rotation. If such a relative equilibrium 
is generated from a fixed-point configuration, any non-zero frequency 
yields a solution. When a double rotation takes place, we can find two 
distinct frequencies, a and (3. Let us factorize the set of these rela- 
tive equilibria by geometric size and by the value of the frequencies. 
In other words, all we care about is the shape of the relative equilib- 
rium. While in the Euclidean case there is only one way to look at the 
shape of a configuration, the concept is a bit more complicated here 
since we discovered orbits, such as those of the curved 6-body problem 
when 3 bodies move along a great circle of a great sphere of §^ and 
the other 3 bodies move along a complementary great circle of another 
great sphere, where there is not only one shape to take into consider- 
ation. Nevertheless, a decomposition of the orbit into a finite number 
of shapes solves the problem. 

We can now ask how many classes of relative equilibria exist in vari- 
ous contexts, such as for each of the 5 types of relative equilibria proved 
to exist or for n = 3, 4, 5, . . . If we restrict to simply asking whether the 
number of such classes is finite or infinite, and if it is infinite, whether 
the set is discrete or contains a continuum, then the question extends 
Smale's Problem 6, discussed at the beginning of Part 5, to spaces of 
constant curvature. But even finding new classes of relative equilibria 
for various values of n and mi, rri2, ■ ■ ■ , m n > is a problem worth 
pursuing. 



22.2. Stability of periodic and quasiperiodic orbits. We referred 
in Section 21 to the stability of Lagrangian solutions in the case of 
curvature k — 1. Since Theorem [4] is proved only for that particular 
case, it would be interesting to solve the general problem in order to 
understand how the stability of these orbits changes for various values 
of the curvature k, using a method similar to the one Martinez and 
Simo developed in [61]. The same method can be further applied to all 
the periodic orbits we constructed in Part 5 of this paper. Moreover, 
Martinez and Simo analyzed in [UT], for k = 1, the stability of a La- 
grangian homographic orbit, which turns out not to be periodic, but 
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quasiperiodic, as shown in [24J. Again their method could be further 
applied to other quasiperiodic orbits, such as the relative equilibria 
constructed in Subsections 17.4 and 17.5[ towards getting a better un- 



derstanding of how curvature affects the stability of orbits. In the mean 
time, it has been used to determine the stability of tetrahedral orbits 
in § 2 , where one body, of mass M, is fixed at the north pole, while the 
other bodies, all of mass m, lie at the vertices of an equilateral triangle 
that rotates in a plane parallel with the plane of the equator, [23] . This 
paper also opens new avenues, such as the study of the orbit's stability 
in § 3 . 

Acknowledgment. The research presented above was supported in 
part by a Discovery Grant from NSERC of Canada. I am indebted to 
the anonymous referee for his/her thorough report, which helped me 
improve the content of this paper. 



Relative equilibria in the 3-dimensional curved n-body problem 



97 



References 

[1] R. Abraham and J.E. Marsden. Foundations of Mechanics (2nd ed.) Benjamin- 

Cummings, Reading MA, 1978. 
[2] A. Albouy and A. Chenciner, Le probleme des N corps et les distances 

mutuelles, Inventiones math. 131 (1998), 151-184. 
[3] A. Albouy and V. Kaloshin, Finiteness of central configurations for five bodies 

in the plane, Ann. of Math, (to appear). 
[4] A. Baker, Matrix Groups: An Introduction to Lie Group Theory, Springer 

Verlag, London, 2002. 
[5] J. Bertrand, Theoreme relatif au mouvement d'un point attire vers un centre 

fixe, C. R. Acad. Sci. 77 (1873), 849-853. 
[6] W. Bolyai and J. Bolyai, Geometrische Untersuchungen, Hrsg. P. Stackel, 

Teubner, Leipzig-Berlin, 1913. 
[7] H. Bruns, Uber die Integrale des Vielkorper-Problems Acta Math. 11 (1887), 

25-96. 

[8] H. Cabral and F. Diacu, Classical and Celestial Mechanics: The Recife Lecture, 

Princeton University Press, Princeton, N.J., 2002. 
[9] J.F. Carihena, M.F. Ranada, and M. Santander, Central potentials on spaces 

of constant curvature: The Kepler problem on the two-dimensional sphere S 2 

and the hyperbolic plane H 2 , J. Math. Phys. 46 (2005), 052702. 
[10] J.F. Carihena, M.F. Rahada, and M. Santander, The Kepler problem and 

the Laplace-Runge-Lenz vector on spaces of constant curvature and arbitrary 

signature, Qual. Theory Dyn. Syst. 7 (2008), 87-99. 
[11] J. Chazy, La theorie de la relativite et la mecanique celeste, Gauthier-Villars, 

Paris, 1930. 

[12] A. Chenciner and R. Montgomery, A remarkable periodic solution of the three- 
body problem in the case of equal masses, Ann. of Math. 152 (2000), 881-901. 

[13] A. Chenciner and H.J. Perez, Angular momentum and Horn's problem, 
larXiv:1110.5030 v2. 

[14] H.S.M. Coxeter, Regular poly topes, Methuen & Co., London, 1948. 

[15] H.S.M. Coxeter, Introduction to Geometry, 2nd edition, John Wiley & Sons, 
New York, 1969. 

[16] T. Damour, M. Soffel, and C. Xu, General-relativistic celestial mechanics. I. 
Method and definition of reference systems, Phys. Rev. D 43, 10 (1991), 3273- 
3307. 

[17] T. Damour, M. Soffel, and C. Xu, General-relativistic celestial mechanics. II. 

Translational equations of motion, Phys. Rev. D 45, 4 (1992), 1017-1044. 
[18] T. Damour, M. Soffel, and C. Xu, General-relativistic celestial mechanics. III. 

Rotational equations of motion, Phys. Rev. D 47, 8 (1993), 3124-3135. 
[19] T. Damour, M. Soffel, and C. Xu, General-relativistic celestial mechanics. IV. 

Theory of satellite motion, Phys. Rev. D 49, 2 (1994), 618-635. 
[20] F. Diacu, Near-collision dynamics for particle systems with quasihomogeneous 

potentials, J. Differential Equations 128 (1996), 58-77. 
[21] F. Diacu, On the singularities of the curved n-body problem, Trans. Amer. 

Math. Soc. 363, 4 (2011), 2249-2264. 
[22] F. Diacu, Polygonal homographic orbits of the curved 3-body problem, Trans. 

Amer. Math. Soc. 364, 5 (2012), 2783-2802. 



98 Florin Diacu 

[23] F. Diacu, R. Martinez, E. Perez-Chavela, and C. Simo, On the stability 

of tetrahedral relative equilibria in the positively curved 4-body problem, 

|arXiv:1204.5729| 32 pages. 
[24] F. Diacu and E. Perez-Chavela, Homographic solutions of the curved 3-body 

problem, J. Differential Equations 250 (2011), 340-366. 
[25] F. Diacu, E. Perez-Chavela, and J.G. Reyes Victoria, An intrinsic approach 

in the curved n-body problem. The negative curvature case, J. Differential 

Equations 252 (2012), 4529-4562. 
[26] F. Diacu, E. Perez-Chavela, and M. Santoprete, The n-body problem in spaces 

of constant curvature, arXiv:0807.1747, 54 p. 
[27] F. Diacu, E. Perez-Chavela, and M. Santoprete, The n-body problem in spaces 

of constant curvature. Part I: Relative equilibria, Journal of Nonlinear Science 

22, 2 (2012), 247-266. 
[28] F. Diacu, E. Perez-Chavela, and M. Santoprete, The n-body problem in spaces 

of constant curvature. Part II: Singularities, Journal of Nonlinear Science 22, 

2 (2012), 267-275. 

[29] F. Diacu, E. Perez-Chavela, and M. Santoprete, Saari's conjecture for the 
collinear n-body problem, Trans. Araer. Math. Soc. 357, 10 (2005), 4215-4223. 

[30] F. Diacu, T. Fujiwara, E. Perez-Chavela, and M. Santoprete, Saari's homo- 
graphic conjecture of the 3-body problem, Trans. Amer. Math. Soc. 360, 12 
(2008), 6447-6473. 

[31] F. Dillcn and W. Kiihnel, Ruled Weingarten surfaces in Minkowski 3-space, 

Manuscripta Math. 98 (1999), 307-320. 
[32] C. Doran and A. Lasenby, Geometric Algebra for Physicists, Cambridge Univ. 

Press, 2003. 

[33] A. Eddington and G.L. Clark, The problem of n bodies in general relativity 

theory, Proc. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 166 (1938), 465-475. 
[34] A. Einstein, L. Infeld, and B. Hoffmann, The gravitational equations and the 

problem of motion, Ann. of Math. 39, 1 (1938), 65-100. 
[35] B. Elmabsout, Sur l'existence de certain configurations d'equilibre relatif dans 

le probleme des N corps, Celest. Mech. Dyn. Astr. 41 (1988), 131-151. 
[36] L. Euler, Institutiones calculi differentials cum eius usu in analysi finitorum 

ac doctrina serierum, published as a book in 1755, Opera Omnia, Ser. 1, vol. 

10, (ed. G. Kowalewski), Birkhauser, 1980. 
[37] R. L. Faber, Foundations of Euclidean and Non-Euclidean Geometry, Marcel 

Dekker, New York, 1983. 
[38] V. A. Fock, Sur le mouvement des masses finie d'apres la theorie de gravitation 

einsteinienne, J. Phys. Acad. Set. USSR 1 (1939), 81-116. 
[39] I.M. Gelfand and S.V. Fomin, Calculus of Variations, Prentice-Hall, Engle- 

wood Cliffs, New Jersey, 1963. 
[40] G. Goe, Comments on Miller's "The myth of Gauss's experiment on the Eu- 
clidean nature of physical space," Isis 65, 1 (1974), 83-87. 
[41] L. Garcia Gutierrez and M. Santander, Levi-Civita regularization and geodesic 

flows for the 'curved' Kepler problem, | [a7Xiv:0707.38l"0> 2, 2007. 
[42] H. Hopf, Uber die Abbildungen der dreidimensionalen Sphare auf die Kugel- 

flache, Math. Ann. 104 (1931), 637-665. 
[43] L. Infeld, On a new treatment of some eigenvalue problems, Phys. Rev. 59 

(1941), 737-747. 



Relative equilibria in the 3-dimensional curved n-body problem 



99 



L. Infeld, Quest — The evolution of a scientist, Readers Union Ltd., London, 
1942. 

L. Infeld and A. Schild, A note on the Kepler problem in a space of constant 
negative curvature, Phys. Rev. 67 (1945), 121-122. 

W. Killing, Die Rechnung in den nichtcuklidischcn Raumformen, J. Reine 
Angew. Math. 89 (1880), 265-287. 

W. Killing, Die Mechanik in den nichtcuklidischen Raumformen, J. Reine 
Angew. Math. 98 (1885), 1-48. 

W. Killing, Die Nicht-Eukildischen Raumformen in Analytischer Behandlung, 
Teubner, Leipzig, 1885. 

V. V. Kozlov and A. O. Harin, Kepler's problem in constant curvature spaces, 
Celestial Mech. Dynam. Astronom 54 (1992), 393-399. 

T. Levi-Civita, The relativistic problem of several bodies, Amer. J. Math. 59, 
1 (1937), 9-22. 

T. Levi-Civita, he probleme des n corps en relativite generale, Gauthier-Villars, 
Paris, 1950; or the English translation: The n-body problem in general relativ- 
ity, D. Reidel, Dordrecht, 1964. 

D. Lewis and T. Ratiu, Rotating n-gon/fcn-gon vortex configurations, J. Non- 
linear Sci. 6 (1996), 285-414. 

H. Liebmann, Die Kcgclschnitte und die Planctcnbewcgung im nichteuklidis- 
chen Raum, Berichte Konigl. Sachsischen Gesell. Wiss., Math. Phys. Klasse 
54 (1902), 393-423. 

H. Liebmann, Uber die Zentralbewegung in dcr nichtcuklidische Geometric, 
Berichte Konigl. Sachsischen Gesell. Wiss., Math. Phys. Klasse 55 (1903), 
146-153. 

H. Liebmann, Nichtcuklidische Geometric, G. J. Goschen, Leipzig, 1905; 2nd 

ed. 1912; 3rd ed. Walter de Gruyter, Berlin, Leipzig, 1923. 

R. Lipschitz, Untersuchung eines Problems der Variationrechnung, in welchem 

das Problem der Mechanik enthalten ist, J. Reine Angew. Math. 74 (1872), 

116-149. 

R. Lipschitz, Extension of the planet-problem to a space of n dimensions and 
constant integral curvature, Quart. J. Pure Appl. Math. 12 (1873), 349-370. 
N. I. Lobachevsky, The new foundations of geometry with full theory of paral- 
lels [in Russian], 1835-1838, In Collected Works, V. 2, GITTL, Moscow, 1949, 
p. 159. 

D.W. Lyons, An elementary introduction to the Hopf fibration, Math. Mag. 
76, 2 (2003), 87-98. 

J.E. Marsden and T.S. Ratiu, Introduction to Mechanics and Symmetry, (2nd 
ed.) Springer Verlag, 1999. 

R. Martinez and C. Simo, On the stability of the Lagrangian homographic 
solutions in a curved three-body problem on § 2 , Discrete Contin. Dyn. Syst. 
Ser. B (to appear). 

R. McGehee, Triple collision in the collinear three-body problem, Inventiones 
math. 27 (1974), 191-227. 

A. I. Miller, The myth of Gauss's experiment on the Euclidean nature of phys- 
ical space, Isis 63, 3 (1972), 345-348. 

R. Mosseri and R. DandolofF, Geometry of entangled states, Bloch spheres and 
Hopf fibrations, J. Phys. A 34 (2001), 10243-10252. 



100 Florin Diacu 

[65] M. Nakahara, Geometry, Topology and Physics, Institute of Physics Publish- 
ing, Philadelphia, 1990. 

[66] J. Hano and K. Nomizu, On isometric immersions of the hyperbolic plane into 
the Lorentz-Minkowski space and the Monge- Ampere equation of certain type, 
Math. Ann. 262 (1983), 245-253. 

[67] E. Perez-Chavela and J.G. Reyes Victoria, An intrinsic approach in the curved 
n-body problem. The positive curvature case, Trans. Amer. Math. Soc. 364, 
7 (2012), 3805-3827. 

[68] L. Perko and E. Walter, Regular polygon solutions of the iV-body problem, 

Proc. Amer. Math. Soc. 94 (1985), 301-309. 
[69] W.F. Reynolds, Hyperbolic geometry on a hyperboloid, Amer. Math. Monthly 

100, 5 (1993), 442-455. 
[70] D. Saari, Collisions, Rings, and Other Newtonian N -Body Problems, American 

Mathematical Society, Regional Conference Series in Mathematics, No. 104, 

Providence, RI, 2005. 
[71] E. Schering, Die Schwerkraft im Gaussischen Raume, Nachr. Konigl. Gesell. 

Wiss. Gottingen 13 July, 15 (1870), 311-321. 
[72] E. Schering, Die Schwerkraft in mehrfach ausgedchnten Gaussischen und Riem- 

manschen Raumen, Nachr. Konigl. Gesell. Wiss. Gottingen 26 Febr., 6 (1873), 

149-159. 

[73] E. Schrodinger, A method for determining quantum-mechanical eigenvalues 

and eigenfunctions, Proc. Royal Irish Acad., Sect. A 46 (1940), 9-16. 
[74] P.J. Serret, Theorie nouvelle geometrique et mecanique des lignes a double 

courbure, Librave de Mallet-Bachelier, Paris, 1860. 
[75] S. Sgries, Erinnerung and "dunkle Seiten," Gedenkkolloquium fur den 

jiidischen Mathematiker Heinrich Liebmann, Heidelberger Nachrichten, June 

19, 2008. 

[76] A.V. Shchepetilov, Comment on "Central potentials on spaces of constant 
curvature: The Kepler problem on the two-dimensional sphere S 2 and the 
hyperbolic plane H 2 ," [J. Math. Phys. 46 (2005), 052702], J. Math. Phys. 46 
(2005), 114101. 

[77] A.V. Shchepetilov, Reduction of the two-body problem with central interac- 
tion on simply connected spaces of constant sectional curvature, J. Phys. A: 
Math. Gen. 31 (1998), 6279-6291. 

[78] A.V. Shchepetilov, Nonintegrability of the two-body problem in constant cur- 
vature sp aces, J. Phys. A: Math. Gen. V. 39 (2006), 5787-5806; corrected 
version at |math.DS/0601382| 

[79] A.V. Shchepetilov, Calculus and mechanics on two-point homogenous Rie- 
mannian spaces, Lecture notes in physics, vol. 707, Springer Verlag, 2006. 

[80] S. Smale, Mathematical problems for the next century, Mathematics: frontiers 
and perspectives, Amer. Math. Society, Providence, RI, pp. 271-294. 

[81] A. F. Stevenson, Note on the Kepler problem in a spherical space, and the 
factorization method for solving eigenvalue problems, Phys. Rev. 59 (1941), 
842-843. 

[82] A. Wintner, The Analytical Foundations of Celestial Mechanics, Princeton 
University Press, Princeton, N.J., 1941. 



